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UNDERWATER  ELECTROMAGNETIC  TRANSMISSION  AND  ITS 
APPLICATION  TO  SHIP  POSITION  CONTROL 

by  E.G.C.  Burt  and  L.  Rigby 
Imperial  College,  London 


ABSTRACT 

The  paper  discusses  the  possible  use  of  electromagnetic  methods  for  the 
position-detecting  element  in  a  ship  control  system,  where  the  need  is  to  keep  the 
ship  over  a  selected  point  on  the  ocean  floor  in  deep  water  -  e.g.  for  sea  bed 
drilling  operations.  Such  systems  currently  employ  acoustic  methods  for  the  deter¬ 
mination  of  position;  receivers  attached  to  the  ship’s  hull  detect  the  arrival  times 
of  pulses  from  an  ultrasonic  transmitter  on  the  ocean  floor,  the  time  differences 
giving  a  measure  of  the  lateral  positional  error.  The  method  however  suffers  from 
the  effects  of  anomalous  acoustic  propagation,  which  can  lead  to  gross  errors  in  the 
estimate  of  position. 

Electromagnetic  waves  in  a  conducting  medium  such  as  salt  water  are  attenuated 
to  a  far  greater  extent  than  acoustic  waves  at  all  but  the  lowest  frequencies;  on 
the  other  hand,  it  is  to  be  expected  that  at  these  long  wavelengths  the  effects  of 
small-scale  irregularities  in  the  medium  will  be  far  less  pronounced.  The  trans¬ 
mission  of  both  transient  and  s inusoidal  E-M  waves  is  investigated,  and  the  paper  also 
discusses  the  effects  of  the  sea-air  interface  on  wave  propagation  and  diffusion. 

The  results  of  some  preliminary  experiments  are  given,  and  their  relevance  to  the 
ship  control  problem  is  briefly  examined. 

TRANSIENT  FIELDS  IN  A  CONDUCTING  MEDIUM 

Linear  Source 

We  consider  first  the  simplest  possible  E-M  source:  a  current-carrying  linear 
element  of  length  dJ l.  Expressions  have  been  derived  for  the  electric  and  magnetic 
field  produced  by  such  a  source  when  the  current  is  a  6-function;  for  example,  the 
magnetic  field  is  given  by 
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with  similar  expressions  for  E_  and  E.  (H  “  H„  -  E.  -  0).  The  element  dP  is  at  the 
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origin  of  coordinates  (r,  0,  <f>)  and  directed  along  the  polar  axis;  the  medium  has 
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conductivity  O,  permeability  y  and  permittivity  e,  with  3 
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and  the  argument  of  the  Bessel  functions  I 
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3(t  -a  )  .  These  formulae  are  derived  on  the  assumption  that  the  medium  is  infinite 

in  extent;  so  they  only  apply  if  both  source  and  field  points  are  clear  of  bounding 
surfaces  between  different  media  -  e.g.  the  ocean  floor  or  the  sea  surface  (but  see 
below) . 

It  will  be  seen  that  the  response  at  a  distance  r  from  the  source  is  zero  until 

t  *  —  ,  at  which  time  there  is  a  <5  impulse  of  H, ,  whose  coefficient  is  however 
v  <p 

enormously  attenuated  with  distance  from  the  source  (the  term 
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2~)r  ,  which  for  sea  water  becomes  e  *  when  r  is  a  few 


centimetres).  The  main  response  is  due  to  the  term  involving  I  and  1^  ,  propagating 
at  a  speed  which  decreases  rapidly  with  distance. 

This  behaviour  may  be  seen  in  Fig.  1,  which  shows  the  magnetic  field  due  to  a 
unit  step  of  current,  obtained  by  integrating  equation  (1): 
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where  we  have  used  normalised  time  t  and  range  r,  given  by 
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It  will  be  observed  that,  at  short  ranges,  the  response  approximates  to  a 
delayed  step,  but  that  as  the  range  increases  the  sharp  rise  disappears,  the  time 
rate  of  change  diminishing  rapidly  with  distance. 

The  magnitude  of  the  initial  step  at  r  =  t  (i.e.  at  t  =  ~  )  is 
e‘r  (1  +  r  * 

which  is  almost  zero  for  r  =  10.  In  sea  water,  for  which  we  may  take 

0*4  ohms  *  m  \  V  *  4  xlO  ^  Hm  ^ ,  and  G  38  —  x  10  ^  Fm  ^  ,  we  have  r  =  1.2  cm 

*  36  *o 

and  t  *  .35  ns;  so  the  step  vanishes  at  ranges  greater  than  about  12  cm  from  the 
source . 

We  shall  be  concerned  with  ranges  and  times  which  are  several  orders  greater 

than  r  and  t  .  Under  these  conditions  (1)  and  (2)  reduce  to 
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for  a  6-function  and  a  step  of  current 
derived  by  putting  e  =  0,  but  they  are 
0t  >>  1.  In  fact,  a  finite  value  of  t 
everywhere  finite,  whereas  (3)  and  (4) 


respectively.  These  results  could  have  been 
valid  for  any  e  provided  that  t  »  a  and 
ensures  that  the  propagation  velocity  is 
imply  an  infinite  speed  near  the  source. 


Ring  Source 


The  field  due  to  an  impulsive  current  in  a  ring  of  radius  "a"  can  be  deduced 
from  the  dipole  results  given  above.  If  the  range  considered  is  much  greater  than 
the  size  of  the  ring  (which  will  certainly  be  so  in  the  application  we  are  contem¬ 
plating)  it  is  found  that 
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where  the  coordinates  are  now  cylindrical  (p,  z) ,  and  the  centre  of  the  current 
ring  is  at  the  origin,  the  plane  of  the  ring  being  normal  to  the  z  axis. 

,  _  \io  ,  2  2.  . 

s  ~  ^  \P  z  )  is  proportional  to  the  square  of  the  distance  from  the  source. 

The  field  components  have  been  written  as  WR  ,  etc,  because  they  are  in  fact 

z 

weighting  functions;  and,  since  the  E-M  equations  are  linear,  the  fields  due  to  any 
current  f\t)  which  is  zero  for  t<0  can  be  found  immediately  as 


Hz(p»z»t)  ■  |  WH  (p,z,x)  f(t-x)  dx 
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Equations  (5)  show  that  the  form  of  the  transient  field  components  as  a 
function  of  time  depends  strongly  on  p  and  z.  As  these  become  larger  the  transient 
persists  for  a  longer  time,  as  well  as  diminishing  in  amplitude.  The  outward 
velocity  of  the  transient  also  falls  with  range,  and  is  different  in  different 
directions.  For  example,  the  time  rate  of  change  of  the  axial  component  H  (p=0) 

reaches  a  peak  value  at  z  when  t  =  .047uaz2  ,  and  the  velocity  of  propagation  of 

the  peak  is  l/.094UOz.  For  p  -  4nxU)'7  henries  m"1  ,  o  =  4  ohms'1  m'1  (typical  of 
sea-water)  and  z  -  100  m,  the  peak  arrives  at  2.3  ms  and  has  a  speed  of  21  km/s  at 
this  point.  At  1  km  the  peak  occurs  at  .23s,  with  a  speed  of  only  2  km/s.  This 
clearly  shows  the  diffusion-like  character  of  the  propagation,  which  is  dominated  by 
the  conductivity  term. 

The  transfer  functions  for  the  field  components  -  i.e.  the  Laplace  transforms 
of  the  weighting  functions  -  can  be  obtained  from  (5)  as 
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and  patting  s  *  iu>  gives  the  steady-state  results  for  a  sinusoidal  source  current. 

THE  USE  OF  E-M  FIELDS  IN  SHIP  POSITION  CONTROL 

The  depths  of  interest  for  sea-bed  drilling  operations  are  from  100m  -  300m. 

If  at  these  ranges  one  or  more  of  the  field  components  can  be  detected  for  reasonable 
values  of  power  supplied  to  the  source,  we  have  in  principle  a  direct  replacement 
for  t*  e  acoustic  system.  The  pulsed  E-M  method  gives  arrival  time  differences  of 
the  same  order  (ms)  as  for  the  acoustic  case,  and  could  therefore  be  us  id  in 
exactly  the  same  way  to  derive  error  information.  Either  a  linear  element  or  a 
solenoid  as  a  detector  would  give  a  voltage  proportional  to  the  rate  of  change  of 
the  magnetic  field;  altemativ*  ly,  a  linear  array  could  be  employed  to  detect  the 
electric  field. 

In  addition  to  this  method  based  on  pulse  arrival  times,  there  is  the  possi¬ 
bility  of  using  amplitude  and  phase  characteristics  for  a  sinusoidal  excitation  of 
the  source.  From  (7),  the  amplitudes  of  Hz  and  H  for  a  current  ring  are 
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where  we  have  written 

x  -  U&a)1  -  (~^)S  (p2+z2)i 

If  the  transmitting  coil  or  solenoid  is  at  a  depth  z  below  the  surface  with 
its  axis  vertical, the  magnitude  of  the  voltages  induced  in  vertical  and  horizontal 
coils  located  at  any  point  on  a  circle  of  radius  p  centred  on  a  point  vertically 
above  the  transmitter  will  be  proportional  to 

w|h  |  and  <*>| H  |  respecti vely.  In  particular,  for  p  =  0  we  have  ] H  |  =  0  and 
z  p  p 

«|Hx|  =  (1  +  2x  +  2x2)^e  X 

2z 


so  for  any  z  there  is  an  optimum  frequency  which  maximises  the  received  voltage, 
given  by 


yoz2 


If  account  is  taken  of  the  "Q"  of  a  tuned  receiving  coil,  the  voltage  is  more 

2  2 
nearly  proportional  to  (0  |H  ] ,in  which  case  the  optimum  frequency  is  u  =  48/pQz  . 


For  z  =  100m,  this  is  152  Hz,  and  for  z  =  300m  it  is  17  Hz.  These  very  low 
frequencies  result  of  course  from  the  attenuation  factor 


exp 
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Fig.  2  exhibits  | H ^ [  and  | f  as  functions  of  p/z.  On  the  sea  surface  the 
contours  of  constant  [ {  and  |h^|  are  circles  centred  on  the  point  vertically 

above  the  transmitter,  the  former  decreasing  and  the  latter  increasing  with 
distance.  These  contours  could  be  utilized  as  shown  in  Fig.  4,  where  the  ship  is 
equipped  with  four  receivers:  a  vertical  coil  at  the  bow  measuring  | | ^ ,  say, 

and  a  similar  coil  at  the  stern  giving  ,  with  two  horizontal  mutually 

perpendicular  coils  on  both  the  port  and  starboard  sides,  measuring  [H^l^  and  |  |  ^ 

respectively.  The  fore  and  aft  receivers  give  voltages  proportional  to  and 

(Fig.  4)  which  in  turn  give  the  error  in  ship  axes  both  in  magnitude  and  direction, 
apart  from  a  180°  ambiguity.  This  is  resolved  by  the  lateral  receivers:  the 
desired  position  lies  to  starboard  if  the  starboard  receiver  gives  the  smaller 
reading  and  vice  versa. 

The  derivation  of  the  error  vector  in  ship  axes  avoids  the  need  for  axis 
transf ormation  and  the  measurement  of  ship's  heading  with  respect  to  a  space  datum; 
the  effects  of  yaw  are  automatically  accounted  for  in  the  measurement  system.  For 
example,  when  the  ship  is  in  the  desired  position  the  measured  error  remains  zero 
for  all  yaw  angles  ard  no  power  need  be  wasted  in  maintaining  an  adverse  heading 
with  respect  to  wind,  wave  or  tide. 

We  have  so  far  taken  no  account  of  changing  z  (due  to  surge)  or  of  pitch  and 
roll.  However,  there  is  sufficient  redundancy  in  the  system  outlined  above  to 
eliminate  their  effects  -  in  fact,  the  received  signals  furnish  an  estimate  of  z  as 
well  as  the  error  D.  There  is  also  further  information  available  in  the  phase  of  the 
signals:  these  are  shown  in  Figure  3. 

SEA  SURFACE  EFFECTS 

The  foregoing  results  are  derived  on  the  assumption  that  the  conducting  medium 
is  infinite  in  extent,  so  both  source  and  field  points  are  taken  to  be  remote  from 
boundaries  with  different  media.  This  however  is  by  no  means  the  case  for  the 
application  which  we  have  in  mind;  the  transmitter  will  be  located  near  the  sea 
floor,  while  the  receiver  will  be  just  below,  or  above,  the  sea  surface. 

The  presence  of  such  boundaries  leads  to  considerable  complication,  since  the 
field  in  each  medium  depends  on  the  properties  of  both  media.  Moreover,  there  are 
four  different  cases  to  consider:  magnetic  or  electric  dipole  transmitters,  with 
the  dipole  axis  either  normal  to  or  parallel  with  the  surface  separating  the  two 
media  (the  field  from  any  other  source  can  be  obtained  as  a  linear  combination  of 
the  fields  due  to  these  four  sources). 

For  the  general  case  in  which  both  media  have  different  conductivities, 
emissivities  and  permeabilities,  the  cylindrical  components  of  the  fields  can  be 
formulated  as  integrals  involving  Bessel  functions.  Tf  the  two  media  are  salt 
water  and  air  the  conductivity  of  the  former  is  the  dominant  factor,  and  the 
emissivities  of  both  media  can  be  neglected.  In  this  situation  analytical  solutions 
have  been  found  for  the  surface  field  components  due  to  any  of  the  four  sources 
mentioned  above  when  they  are  excited  harmonically.  Comparison  with  computer 
solutions  of  the  full  equations  have  confirmed  that  the  approximate  solutions  are 
very  accurate  at  the  low  frequencies  of  interest. 

Tn  general  the  total  field  may  be  ascribed  to  the  direct  field  from  the  source 
(as  if  there  were  no  surface),  a  similar  field  from  an  imago  source ,  and  a  third 
source  due  to  electrical  currents  distributed  over  the  surface.  Since  we  have 
assumed  the  same  magnetic  permeability  for  each  medium,  all  components  except  F. 
are  continuous  at  the  surface.  To  illustrate  the  effect  of  the  sea-air 
surface  we  take  the  particular  case  of  a  vertical  magnetic  dipole  as  the  source. 


AP  1-9 


ON  IDENTIFICATION  OF  NONLINEAR 
SPEED  EQUATION  FROM  FULL  SCALE  TRIALS 

by 

M.  Blanke 

Technical  University  of  Denmark 
Lyngby,  Denmark 


ABSTRACT 

The  paper  describes  a  novel  approach  on  identification  of  the 
nonlinear  dynamic  speed  equation  for  a  surface  ship  from  full  scale 
measurements.  Different  methods  of  analysis  are  discussed,  and  the 
application  of  frequency  analysis  techniques  is  shown  to  be  attract¬ 
ive  regarding  separation  of  the  effects  of  different  nonlinear  terms 
and  concerning  measurement  of  small  amplitude  loss  terms. 

The  findings  reported  are  the  relative  large  influence  on  pro¬ 
pulsion  from  steering  activity  due  to  wake  fluctuations  and  the 
practical  demonstration  of  hydrodynamic  memory  effects  being  import¬ 
ant. 


The  dynamical  speed  equation  is  finally  discussed  in  the  paper 
concerning  autopilot  performance. 


INTRODUCTION 

Speed  loss  of  vessels  due  to  steering  activity  is  achieving  an 
increasing  interest  because  significant  reductions  in  operating 
costs  can  be  obtained  by  even  minor  improvements  in  propulsion  ef¬ 
ficiency  during  the  long  periods  of  automatic  course  keeping  at  sea. 

Loss  of  speed  is  due  to  the  movements  of  the  hull  and  deflect¬ 
ions  of  the  rudder  but  also  the  influence  on  propeller  performance 
from  sway  and  yaw  is  contributing.  In  undisturbed  conditions,  the 
average  speed  loss  due  to  rudder  activity  and  hull  motions  can  ac¬ 
curately  be  expressed  using  existing  mathematical  models  ([1]  and 
[2]),  but  under  dynamic  conditions  the  effects  of  varying  inflow 
velocity  to  the  propeller  and  its  interaction  with  the  engine  must 
be  considered.  Such  dynamic  conditions  are  met  e.g.  when  waves  and 
wind  influence  the  coursekeeping. 

The  accurate  modelling  of  the  combined  dynamic  performance  of 
ship,  propeller  and  engine  is  essential  when  optimum  performing  au¬ 
topilots  are  to  be  designed,  in  particular  when  the  governor  control 
loop  is  to  be  attached  to  the  course  keeping  loop  as  to  achieve  mi¬ 
nimum  propulsion  losses. 

The  paper  is  concerning  verification  of  such  a  mathematical  mo¬ 
del  by  means  of  steering  experiments  on  full  scale  ships. 


medium  is  to  introduce  a  small  phase  change  in  gradually  increasing  with 

increasing  p,  and  to  smooth  out  the  sudden  phase  reversal  of  H  ,  the  degree  of 
smoothing  depending  on  the  value  of  0. 

For  these  experiments  the  measured  conductivity  was  2.5  reciprocal  ohms  per 
metre,  and  Figures  12(b)  -  15(b)  show  the  theoretical  curves  of  phase  as  a  function 
of  P for  this  value  of  <3,  together  with  the  experimentally  measured  values  at  3  ft 
intervals.  The  transmitter  and  receiver  angles  are  those  given  above  for  the 
amplitude  plots. 

It  will  be  seen  that  when  the  receiver  is  fully  submerged  the  measured  and 
theoretical  values  for  the  phase  of  Hz  are  in  close  agreement  (Fig.  13(b)).  As  the 

distance  from  the  source  increases  (Figs.  12(b)  and  14(b))  the  agreement  becomes 
less  good:  however,  the  amplitudes  for  large  p  are  then  very  small,  so  it  is 
difficult  to  obtain  a  good  measure  of  phase  (the  ratio  of  two  small  quantities) 
with  the  analogue  equipment  used. 

The  expected  change  in  the  phase  of  (Fig.  15(b))  is  less  than  1  degree, 

smaller  than  the  resolution  of  the  signal  processing  unit.  Nevertheless,  the 
increasing  trend  in  the  experimental  points  is  clearly  matched  by  the  theoretical 
curve . 

CONCLUDING  REMARKS 

The  investigations  outlined  in  this  paper  suggest  that  underwater  electro¬ 
magnetic  transmission  provides  a  feasible  means  of  error  detection  for  a  ship 
position  control  system,  and  might  well  offer  advantages  in  terms  of  accuracy, 
simplicity  and  reliability.  The  transmitter  depth  required  for  an  operational 
system  is  of  course  greatly  in  excess  of  that  available  for  the  experimental  work 
described  in  the  paper,  although  the  frequency  used  (23  Hz)  is  appropriate  for 
depths  of  100  -  300  m.  Further  work  is  planned  to  gather  data  under  conditions 
closer  to  those  in  which  the  apparatus  would  actually  be  used  as  part  of  a  control 
system;  in  particular,  the  measurement  of  the  signal  to  noise  ratio  for  greater 
depths  and  for  various  frequencies. 

It  is  also  proposed  to  examine  other  components  of  the  electromagnetic  field 
for  both  magnetic  and  electric  dipole  sources.  Because  of  the  complexity  of  the 
ship  control  problem  it  may  be  desirable  to  take  advantage  of  the  additional  infor¬ 
mation  available  (electric  and  magnetic  components,  horizontal  and  vertical,  and 
their  phases  relative  to  each  other  and  to  the  source)  by  detecting,  correlating 
and  filtering  the  signals  in  an  optimum  way,  taking  account  of  the  different  noise 
structure  of  each  component. 
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Figure  15.  Receiver  near-horizontal  at  water  surface:  Corrected  Resul ts 
Depth  of  transmitter  24.13  ft. 

Transmitter  axis  0.045  radians  from  vertical 
Receiver  axis  0.015  radians  from  horizontal 
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Figure  14.  Receiver  fully  out  of  water:  Corrected  Results 
Depth  of  transmitter  24.13  ft 
Transmitter  axis  0.1  radians  from  vertical 
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Figure  13.  Receiver  just  submerged:  Corrected  Results 
Depth  of  transmitter  24.13  ft. 

Transmitter  axis  0.1  radians  from  vertical 
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Figure  12.  Receiver  Halfway  in  water:  Corrected  Results 
Depth  of  transmitter  24.13  ft 
Transmitter  axis  0.1  radians  from  vertical 
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The  transmitter  was  suspended  vertically  by  its  cable  from  a  small  dinghy,  the 
lower  end  of  the  transmitter  being  attached  to  a  weight  resting  on  the  bed  of  the 
lake  in  order  to  minimise  transmitter  motion  due  to  currents.  To  the  same  end  an 
attempt  was  made  (not  wholly  successfully)  to  prevent  dinghy  motion  by  means  of 
lines  from  dinghy  to  shore.  The  receiver  was  on  a  raft  which  could  be  moved  in  a 
fixed  direction  and  at  known  distances  from  the  point  of  suspension  of  the  trans¬ 
mitter.  To  avoid  saturating  the  receiver  at  the  short  distances  involved 
(maximum  42  ft)  it  was  necessary  to  reduce  the  transmitter  supply  voltage  from 
150V  to  50V,  so  the  power  consumption  of  the  transmitter  was  rather  less  than  1  watt. 

Readings  of  H  and  H  (or  rather  their  time  rates  of  change)  were  taken  at 
intervals  of  3  ft,  measured  along  the  surface  from  the  point  vertically  above  the 
transmitter.  At  each  such  value  of  p  the  vertical  component  was  measured  at 

three  points:  with  the  receiver  (suspended  from  its  cable)  just  submerged,  half¬ 
way  out  of  the  water  and  fully  out  of  the  water.  The  horizontal  measurement  H 
was  taken  with  the  receiver  just  below  the  surface,  and  held  horizontal  (as 
near  as  could  be  judged)  by  means  of  ropes  lashed  to  receiver  and  raft. 

Amplitude  Measurements 

The  results  so  obtained  are  shown  in  Figures  8-11,  where  the  full  lines  are 
the  theoretical  curves.  The  differences  in  the  three  sets  of  readings  for  the 
amplitude  of  H  (Figures  8,  9  and  10)  are  not  of  course  due  to  surface  effects: 
they  are  merely  the  result  of  changes  in  the  z  coordinate. 

While  agreement  between  theory  and  experiment  is  fair,  the  results  are  less 
accurate  than  expected,  even  allowing  for  the  somewhat  rudimentary  nature  of  the 
experimental  arrangement.  Later  investigation  however  showed  that  the  distance 
measurements  were  in  error,  because  the  datum  points  used  on  transmitter  and 
receiver  housings  were  displaced  from  the  centres  of  the  coils.  These  errors 
(4"  for  the  transmitter  and  6"  for  the  receiver)  would  of  course  be  quite  insig¬ 
nificant  at  ranges  of  100  metres  or  so,  but  they  become  important  at  the  short 
distances  involved  in  these  experiments. 

Figures  12  -  15  show  the  effect  of  correcting  the  distance  measurements,  and 
also  of  introducing  small  angles  to  account  for  the  possibility  that  the  trans-* 
mitter  axis  was  not  exactly  vertical  and  that  the  receiver  axis  was  not  exactly 
horizontal  for  the  H  measurements.  The  best  fit  to  the  H  data  was  found  to  be 
P  2 

with  a  deviation  of  the  transmitter  from  the  vertical  of  0.1  radians.  For  the 
data,  which  was  taken  the  following  day,  the  theoretical  curves  shown  are  for  a 
transmitter  angle  of  0.045  radians  from  the  vertical  and  a  deviation  of  the 
receiver  axis  of  0.015  radians  from  the  horizontal. 

Angular  errors  of  this  magnitude  could  certainly  have  existed.  With  the  above 
adjustments  the  differences  between  the  experimental  points  and  the  theoretical 
values  (the  full  lines  in  Figures  12(a)  -  15(a))  are  very  small,  and  are  within 
the  limits  to  be  expected  from  the  known  inaccuracies  inherent  in  the  analogue 
equipment  (multipliers,  operational  amplifiers  etc.)  used  in  the  signal  processing 
unit. 

Phase  Measurements 

As  mentioned  earlier,  the  conditions  for  the  experiments  described  above  are 
such  that  the  amplitude  of  and  are  almost  exactly  as  they  would  be  in  air. 

This  however  is  not  true  of  the  phases:  these  begin  to  show  differences  before  any 
distinction  between  amplitudes  becomes  apparent.  In  air,  at  the  very  long  wave¬ 
lengths  involved,  the  phase  of  wit^  respect  to  the  transmitter  does  not  change 

with  p,  while  simply  reverses  sign  at  p  *  /lz.  The  effect  of  the  conducting 
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Figure  10.  Receiver  fully  out  of  water 

Depth  of  transmitter  23.75  ft. 
Transmitter  axis  assumed  vertical 
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Figure  11.  Receiver  horizontal  near  water  surface 
Depth  of  transmitter  23.75  ft 
Transmitter  axis  assumed  vertical 
Receiver  axis  assumed  horizontal 
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Figure  8.  Receiver  Halfway  in  Water 

Depth  of  transmitter  23.75  ft 
Transmitter  axis  assumed  vertical 
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Figure  9.  Receiver  just  submerged 

Depth  of  transmitter  23.75  ft. 
Transmitter  axis  assumed  vertical 
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mitter,  but  the  coil  has  20,000  turns,  giving  an  inductance  of  136  henries  and 
therefore  a  tuning  capacity  of  about  .  35yF.  The  receiver  unit  also  contains  a 
remote  tuning  device  and  an  active  filter  network,  thus  ensuring  that  the  signal 
transmitted  through  the  cable  to  the  signal-processing  unit  is  at  a  reasonable 
voltage  level. 

The  tuning  is  controlled  from  the  far  end  of  the  cable  by  a  DC  voltage  which 
alters  the  effective  resistance  of  a  field-effect  transistor  incorporated  in  the 
feedback  circuit  of  an  operational  amplifier.  There  is  also  positive  capacitative 
feedback,  so  the  effective  capacitance  across  the  receiver  coil  terminals  is 
adjustable  via  the  DC  voltage.  This  arrangement  has  the  advantage  that  the  cable 
wiring  is  not  included  in  the  tuned  circuit,  thus  obviating  a  possible  source  of 
noise . 


the 


Further  operational  amplifiers  within  the  receiver  tube  provide  a  filter  with 
transfer  function: 


1  +  k2sT 


kisT  ) 

k  sT  +  s2T2  / 


2 


where  k^t  k^  are  chosen  to  give  a  high  gain  at  w  “  2tt .23  *  and  a  band  width  of  a 
few  Hertz  about  this  frequency. 

Signal  Processing  Unit 


As  mentioned  above,  the  transmitter  power  is  supplied  remotely  from  a  DC  source. 
The  firing  of  the  transistor  gives  rise  to  a  small  current  pulse  down  the  cable,  and 
from  this  (after  filtering)  a  sine  wave  can  be  derived,  its  phase  being  linked  to 
that  of  the  transmitter  coil  current,  and  thus  to  the  phase  of  the  generated  mag¬ 
netic  field. 

If  Asinwt  is  the  reference  signal  so  obtained,  and  the  received  signal  is 
Vsin (u)t+4>)  +  N(t),  where  N  is  noise,  cross-correlation  of  the  two  gives 


(Vsin(u)t+<J>)  +  N(t))Asincot  =  AV  cos  <p 

provided  signal  and  noise  are  uncorrelated. 

The  multiplication  is  effected  in  the  signal  processing  unit  by  an  analogue 
multiplier,  and  the  averaging  by  a  simple  CR  network  with  an  adjustable  time 
constant  (2-10s).  A  second  multiplier  is  arranged  to  give  V  sin  <J>  (by  using  the 
reference  signal  shifted  by  90  ),  so  smoothed  values  of  V  and  $  can  be  obtained  as 
DC  voltages. 

EXPERIMENTAL  RESULTS 

In  order  to  test  the  validity  of  the  theoretical  work,  the  apparatus  described 
above  has  been  used  in  a  sheltered  salt-water  lake  (Horsea  Lake,  Portsmouth),  and 
the  results  so  obtained  are  presented  in  Figures  8-15. 

The  experimental  equipment  was  designed  to  operate  at  depths  of  up  to  100  metres 
unfortunately  the  greatest  depth  available  at  the  salt-water  lake  used  for  the 
trials  was  only  about  ten  metres.  For  h  *  10m  and  f  *  23Hz,  the  parameter  h  -  0.15, 
for  which  the  amplitudes  of  Hp  and  are  almost  exactly  as  they  would  be  in  air, 

the  presence  of  the  water  and  of  the  surface  having  a  negligible  effect  (Figure  5a). 
(There  is  however  a  detectable  difference  in  phase:  this  is  discussed  below). 
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and  consider  the  amplitudes  of  the  two  components  (H^  and  H^)  of  the  magnetic  field 

at  the  surface.  The  fields  can  be_conveniently  expressed  as  functions  of  two  non- 
dimensional  parameters  h  *  kh  and  p  *  kp  ,  where  h  is  the  depth  of  the  source  below 

t  ^ 

the  surface  and  k  =  (^p-) 5 .  The  factor  K  *  4n/Mk  where  M  is  the  magnetic  moment 
of  the  source. 


Figures  5,  6  and  7  show  the  amplitudes  of  and  H z  at  the  surface  for 

h  -  0.2,  0.8  and  2.0  respectively;  in  each  case  curves  are  also  given  for  the 
response  in  air  (i.e.,  ignoring  the  conductivity  term)  and  for  that  in  an  infinite 


medium  (ignoring  the  presence  of  the  surface).  It  will  be  observed  that,  for 
h  =  0.2,  the  three  curves  are  identical,  so  the  field  is  virtually  independent  of 


the  medium.  For  h  *  0.8  the  curves  begin  to  diverge,  while_at  h  *  2 .0  they  are 


distinctly  different.  It  is  interesting_to  note  that,  for  h  *  2.0,  the  presence  of 


the  surface  reduces  the  value  of  Hp  for  p  less  than  about  3.5,  but  for  H ^  the 
reverse  is  true  -  the  surface  enhances  the  amplitude  of  this  component. 


If  we  take  p  =  4  ohm  m  and  f  *  25Hz,  then  k  is  about  0.02,  and  the  curves 
are  for  h  =  10,  40  and  100  metres;  so  with  these  parameters  the  surface  effects 
become  noticeable  at  source  depths  greater  than  about  40  metres.  At  a  given  source 
depth  the  effects  become  more  apparent  as  the  frequency  (or  the  conductivity)  is 
increased  -  e.g. ,  at  f  =  100  Hz  the  figures  represent  the  situation  at  depths  of 
5,  20  and  50  metres. 


EXPERIMENTAL  APPARATUS 


As  mentioned  earlier,  the  optimum  frequency  for  the  depths  of  interest  for  ship 
control  (100  -  300m)  ranges  from  about  150Hz  down  to  15Hz.  In  order  to  investigate 
the  feasibility  of  the  method  an  experimental  system  has  been  constructed  to 
operate  on  a  frequency  of  23Hz:  this  should  give  detectable  signals  at  ranges  of 
100m  or  more,  depending  on  the  ambient  noise  level. 

Both  transmitter  and  receiver  units  are  enclosed  in  sealed  plastic  tubes,  to 
which  are  attached  cables  for  the  supply  of  power  to  the  transmitter  and  the 
extraction  of  signals  from  the  receiver.  These  cables  (100m  for  the  transmitter  and 
25m  for  the  receiver)  also  act  as  supports,  and  are  arranged  in  such  a  way  that 
theie  is  no  strain  on  the  cables  at  the  points  where  they  enter  the  tubes.  The 
complete  unit  can  operate  on  dry  batteries  where  mains  power  is  not  available. 

Transmitter  Unit 

The  transmitter  coil  consists  of  635  turns  of  wire  wound  on  a  laminated  iron 
core  having  a  circular  cross-section  3  cm  in  diameter  and  a  length  of  91  cm.  The 
coil  inductance  is  0.4  henries;  a  capacity  of  I20pf  is  therefore  needed  to  give  a 
natural  frequency  of  23Hz.  The  twelve  10yF  capacitors  are  housed  alongside  the  coil 
in  the  plastic  tube,  which  also  contains  a  drive  circuit  to  maintain  oscillations  at 
a  constant  level.  The  circuit  is  fed  via  the  100m  cable  from  a  150  V  DC  supply, 
and  a  current  of  50  mA  (i.e.  a  power  consumption  of  7.5W)  is  sufficient  to  give  a 

B  of  about  10  kilogauss  (near  saturation  of  the  iron).  The  external  magnetic 
max 

field  due  to  the  transmitter  is  proportional  to  the  dipole  moment  Bv,  where  v  is 
the  volume  of  the  iron;  so  the  field  can  be  increased  by  increasing  v,  and  at  the 
same  time  providing  more  power  to  accommodate  the  increased  iron  loss  (also  propor¬ 
tional  to  v  for  a  given  B).  The  dimensions  chosen  here  are  a  compromise  between 
the  desire  for  a  high  field  strength  and  the  need  to  minimise  the  weight  for 
reasonable  portability.  (The  transmitter  unit  weighs  17  kg). 

Receiver  Unit 

The  receiver,  weighing  11  kg,  has  an  iron  core  identical  to  that  of  the  trans- 
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Figure  7.  Amplitude  of  Magnetic  Field  for  h  =  2.0 
A  In  air 

O  In  an  infinite  conducting  medium 
0  In  a  conducting  half-plane  with  an  air  interface 
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Figure  5.  Amplitude  of  Magnetic  Field  for  h  *  0.2 
A  In  air 

O  In  an  infinite  conducting  medium 

ID  In  a  conducting  half-plane  with  an  air  interface 
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A  MATHEMATICAL  MODEL 


The  dynamics  of  ship  manoeuvring  has  been  treated  in  various 
papers.  Mostly,  the  mathematical  models  are  developed  for  predict¬ 
ing  ship  steering  qualities  and  for  the  programming  of  ship  hand¬ 
ling  simulators.  Those  models  are  able  to  be  adapted  to  a  variety 
of  ship  types  and  applicable  over  a  wide  range  of  ship  speeds  and 
engine  revolutions.  They  give  precise  representations  of  manoeuvring 
responses  provided  a  proper  set  of  coefficients  is  obtainable.  The 
equations  appear  as  quite  complex  as  numerous  derivatives,  coeffi¬ 
cients  and  characteristics  are  employed  which  require  sophisticated 
techniques  for  being  determined. 

When  autopilot  design  is  considered,  a  dynamical  model  must  be 
available,  which  exactly  describes  the  structure  of  the  system  to  be 
controlled. Small  uncertainties  in  model  parameters  generally  do  not 
influence  closed  loop  performance,  whereas  lack  of  information  on 
couplings  between  parts  of  the  system  may  prevent  the  proper  funct¬ 
ion  of  a  controller.  A  precise  model  having  too  many  parameters 
which  are  hardly  obtainable  is  therefore  less  feasible  than  one 
which  gives  approximate  responses  only  but  has  the  correct  structure 
and  is  described  by  obtainable  parameters.  Such  simplified  mathema¬ 
tical  models  are  derived  from  the  more  advanced  descriptions  by  e.g. 
considering  only  a  limited  range  of  operational  conditions.  For 
course  keeping  autopilot  design,  the  simplified  yaw  equation  is  an 
example  of  such  a  model  which  has  been  used  during  more  than  a  de¬ 
cade. 


Models  of  reduced  complexity  are  also  required  when  free  sail¬ 
ing  full  scale  experiments  are  used  for  model  validation,  because 
the  parameters  in  the  complete  descriptions  may  not  be  able  to  be 
identified  separately.  Instead,  the  composite  coefficients  enter¬ 
ing  in  descriptions  on  state  variable  form  or  in  the  transfer  funct¬ 
ions  can  be  determined.  Which  parts  of  the  equations  may  actually 
be  identified  depend  strictly  on  the  actuators  and  measurements 
available  for  the  experiments. 

Most  control  systems  design  procedures  are  based  on  models  be¬ 
ing  on  state  variable  form  or  given  as  transfer  functions  from  in¬ 
puts  to  outputs.  When  a  model  is  derived  on  this  form,  it  is  thus 
feasible  for  being  used  both  as  a  basis  of  autopilot  design  and  for 
model  validation  by  input-output  experiments. 

Aiming  design  of  autopilots  which  minimize  propulsion  losses 
a  model  of  the  ship  steering  dynamics  which  include  interaction  with 
screw  and  engine  is  outlined  below.  The  model  is  derived  from  more 
complete  descriptions  by  assuming  a  constrained  set  of  motions  of 
the  hull  in  order  to  relax  the  equations  from  certain  extraneous  non¬ 
linear  terms.  The  form  of  the  model  is  also  chosen  to  be  adequate 
for  identification  of  the  propulsion  dynamics  by  free  sailing  experi¬ 
ments  . 
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Ship  Dynamics 

For  completely  describing  the  motions  of  a  ship  six  degrees  of 
freedom  are  required,  translations  in  three  directions  and  rotations 
around  three  axes.  When  employed  for  describing  the  propulsion  loss 
mechanisms,  however,  it  is  sufficient  to  consider  sway,  yaw  and 
surge  motions  only,  as  the  pitch,  heave  and  roll  dynamics  form  a  set 
of  equations  which  are  only  insignificantly  coupling  to  the  remain¬ 
der  of  the  motional  equations. 

In  closed  loop  autopilot  operation,  coupling  between  roll  and 
yaw  can  be  induced  through  rudder  activity  in  certain  conditions. 

For  full  form  merchant  ships  the  coupling  is  rather  small  and  is 
mainly  perceptible  around  the  natural  roll  eigenfrequency .  Generally 
this  is  much  higher  than  the  band  of  frequencies  where  propulsive 
loss  effects  are  preponderant,  so  by  proper  autopilot  design  the  ef¬ 
fect  can  be  suppressed.  The  roll  equation  is  therefore  further  on 
not  considered. 

The  formulation  of  the  equations  is  well  known  from  the  tho¬ 
rough  description  by  Abkowitz  M  ]  and  Norrbin  [2]  .  It  is  not  in¬ 
tended  to  repeat  details  of  those  works  here,  but  it  is  feasible  to 
use  fragments  of  their  developments  as  a  base  for  discussing  the 
simplifications  possible  when  the  restricted  model  is  derived. 

To  describe  the  equations  of  motion  the  coordinate  system  fixed 
to  the  ship  is  shown  in  fig.  1. 

The  foundation  of  the  dynamical  behaviour  of  a  ship  is  the  con¬ 
servation  of  linear  and  angular  momentum.  Mathematically  this  is  ex¬ 
pressed  by  the  laws  of  Newton  in  eq.  1  where  only  sway,  surge  and 
yaw  motions  are  considered. 


ax 

X 

Y 

X 

X 

= 

+ 

+ 

+ 

a  y 
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abs 

hydrodyn. 

L  -1 

'prop 

L  - 

rudder 

L  J 

(D 

external 


Iz  t^abs  Njhydrodyn.  +  ^prop  +  ^rudder  +  ^ external 


The  right  hand  side  force  and  torque  components  X,  Y  and  N  are 
separated  according  to  their  origin:  Hydrodynamical  effects,  pro¬ 
peller  excitation, rudder  activity  and  external  disturbances  from 
waves  and  wind.  Ship  mass  is  m  and  inertia  around  the  vertical  axis 
I  .  Force  and  torque  components  refer  to  the  directions  of  ship  fix¬ 
ed  coordinates  whereas  the  accelerations  (a  ,  a  ,  A  )  are  relative 
to  an  inertial  frame. 

The  accelerations  are  more  conveniently  expressed  in  terms  of 
the  velocities  relative  to  ship  axes.  If  xG  is  the  x  coordinate  of 
the  center  of  mass,  the  momentum  vector  becomes 
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Space  frame 


Fig  1 .  Definition  of  ship 
fixed  coordinates  and  of 
componentes  of  force  and 
torque . 


5  rudder  angle 
vp  sway  velocity  at  stern 
u  surge  velocity 
v  sway  velocity 
O  angular  velocity  vector  in 
inertial  space. 


Symbols : 

D  screw  diameter 
h  height  of  rudder 
a  is  D/h 

Yp  inflow  velocity  at  screw 
U*  flow  velocity  downstream 
due  to  screw. 

V.+k-U.  flow  velocity  past 
R  rudder 


Fig  2.  Geometry  of  screw  and  rudder.  Flow  velocities  are 
relative  to  ship. 
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(2) 


max 

p 

m(il-rv-xGr  ) 

may 

= 

m(^+ru+xGf ) 

JzP. 

Izf'+mxG('>+ru) 

abs 


Below  the  force  and  torque  components  in  the  motional  equation  are 
treated  separately  according  to  their  origin. 

Hydrodynamic  Forces.  The  accelerations  are  important  in  the 
determination  of  transient  hydrodynamic  performance  of  the  hull.  Be¬ 
cause  exchange  of  energy  between  the  hull  and  the  streaming  water 
is  involved  in  changes  of  ship  velocity  vector  (it,  $,  t)  effects  of 
virtual  added  mass  and  added  inertia  occur  in  the  equations .  Terms 
equivalent  to  the  acceleration  components,  eq.  2  are  in  fact  present, 
but  the  magnitude  of  each  coefficient  has  to  be  separately  deter¬ 
mined. 

Steady  state  velocity  also  gives  rise  to  resistance  in  a  vis¬ 
cous  fluid.  The  reaction  on  the  hull  from  the  steady  forward  advance 
is  in  essence  the  resistance  of  the  wet  surface.  This  is  proportio¬ 
nal  to  the  square  of  the  forward  speed.  Wavemaking  and  other  effects 
further  increase  the  resistance  at  high  speeds. 

Reactions  from  steady  sway  and  yaw  velocity  can  be  approximated 
by  the  readily  calculated  forces  and  torques  on  a  slender  wing  of 
equivalent  dimensions  as  the  hull.  Those  reactions  appear  to  vary  bi- 
linearly  with  speed  times  sway  or  rate  of  turn.  By  large  excursions 
in  sway  and  yaw  rate,  however,  nonlinear  terms  are  predominant  due 
to  cross  flow  effects  and  at  high  frequencies  of  excitation,  the  re¬ 
actions  depend  also  on  the  tiraehistory  of  motion. 

Here  the  objective  is  to  utilize  the  equations  for  operative 
conditions  as  course  keeping  only.  The  inclusion  of  the  simple  bi¬ 
linear  terms  therefore  suffices  for  the  sway  and  yaw  rate  equation. 

The  hydrodynamic  forces  employed  in  the  equation  of  motion  thus 

are 
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In  the  equation  the  derivative  jN/su^r  has 
Nur  etc- 


been  abbreviated 
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Screw  Behaviour.  The  Propeller  action  in  a  field  of  steady  in¬ 
flow  velocity  is  described  by  relations  between  developed  thrust  T, 
propeller  shaft  torque  Q|  and  velocity  components  according  to  fig. 2. 
By  using  momentum  theory  and  considering  the  lift  developed  by  the 
screw  blades,  thrust  and  torque  can  be  shown  to  depend  on  bilinear 
and  quadratic  terms  in  screw  angular  velocity  n  and  inflow  velocity 
at  the  propeller  VA .  During  the  steady  course  keeping,  fluctuations 
in  propeller  load  are  not  excessive  compared  to  the  average  load. 

From  practical  experience  then  propeller  thrust  and  torque  can  be 
readily  assessed. 

Ship  screw  performance  is  commonly  presented  in  nondimension 
form  by  the  torque  and  thrust  coefficient  KQ  and  KT  varying  over  the 
range  of  propeller  loads  as  reflected  in  the  advance  number  J.  For 
the  small  load  variations  considered,  KQ  and  KT  can  be  expressed  sa¬ 
tisfactory  accurate  merely  using  a  linear  approximation 


Kt  -  ICjq  +  Kjj  J 
KQ  =  KQ0  +  KQJ  J 

where 

%,  =  T  /  (p4n2n2D4) 

Kq  =  Q1/  (p4n2n2D5) 

J  =  VA/  (27inD) 

D  is  propeller  diameter 

n  is  propeller  angular  velocity  in  rad/sec. 


(4) 


In  dimensional 

form  we 

get 

T  =  T 

2 

n  +  T  „ 

nV, 

nn 

nv 

A 

Q,  =  Q  „ 

n2  +  Q 

nV. 

1  nn 

nv 

A 

(5) 


The  inflow  velocity  at  the  propeller  VA  is  related  to  the  ship  for¬ 
ward  speed  by  the  wake  coefficient  w  as 


VA  =  (!-w)  u 


(6) 


where  u  is  the  ship  forward  speed. 

The  propeller  action  is  less  predictable  when  the  inflow  veloci¬ 
ty  is  varying.  Very  few  surveys  have  ever  been  made  of  the  screw  be¬ 
haviour  when  subject  to  varying  sway  velocity,  and  due  to  the  highly 
complex  boundary  layer  interaction  with  propeller  and  stern,  no  simp¬ 
le  analytical  expression  describes  the  phenomena  concisely.  Etapiri- 
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cally,  the  effect  of  local  sway  velocity  at  the  propeller  can  be 
expressed  by  its  influence  on  the  average  wake  coefficient 


w  =  wQ  +  W(r,v) 


(7) 


For  a  single  screw  vessel,  typically  a  plot  of  W(r,v)  like 
fig.  3  will  result  [7].  The  wake  is  seen  to  vary  as  an  odd  function 
of  local  sway  velocity  in  the  range  normally  exploited  at  steady 
course.  The  odd  dependancy  of  wake  coefficient  from  local  sway  is  un¬ 
derstood  to  be  the  build  up  of  antisymmetric  surge  and  pressure 
fields  due  to  the  one  way  screw  rotation. 

The  sway  at  the  propeller  is  a  composite  of  ship  sway  v  and 
yaw  rate  r,  vp  =  v  -  l^-r.  The  coefficent  1PG  is  the  geometric 
distance  from  centre  of  gravity  to  the  propeller  position. 


Linearization  of  the  wake  function  is  obtained  by  expanding  in 
first  order  components  in  ship  sway  and  rate  of  turn 


w 


w0  +  wvv  -  wvlpGr 


(8) 


The  inflow  velocity  to  the  propeller  thus  becomes 


VA  =  (l"wo)  u  "  u 


(9) 


Equations  (5)  and  (9)  then  constitute  the  description  of  pro¬ 
peller  performance. 


Rudder  forces.  A  rudder  may  be  considered  a  small  aspect  ratio 
wing.  'The  lift  force  then  vary  with  the  square  of  advance  velocity 
at  the  rudder  and  for  small  rudder  angles  it  is  proportional  to 
rudder  deflection  6 .  Accordingly  rudder  drag  is  proportional  to  the 
square  of  rudder  deflection. 

Transformation  of  rudder  force  to  ship  fixed  coordinates  is  re¬ 
quired  for  fitting  in  the  equations  of  motion  for  the  ship.  This 
transformation  and  approximation  of  sine  and  cosine  terms  by  assum¬ 
ing  small  rudder  deflections  result  in  the  rudder  forces 


XR  ~  Xccss  °2  &2 


YR  =  Ycc6  °2  a 


(10) 


The  torque  exposed  to  the  hull  from  the  rudder  is  expressed 
through  the  distance  ljy,  from  the  rudder  to  centre  of  gravity 

Nr  =  1RGYR  =  XRG  Ycc$  c2a  ^ 1 1 


Here  c  is  the  mean  flow  velocity  at  the  rudder.  The  velocity  c 
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accounts  for  the  effects  of  propeller  loading  when  expressed  as  fol¬ 
lows  [8] 


c2  =  a(VA  +  kRUA)2  +  (1-a)  V* 

where 


(12) 


a  is  the  ratio  between  screw  diameter  and  height  of  rudder 
kRUA  is  the  propeller  induced  velocity  at  the  rudder 


The  term  UA  is  the  theoretical  screw  induced  velocity  infinite¬ 
ly  downstream  and  is  expressed  from  momentum  theory  by 

UA  +  f  -  VA  03) 

Model  tests  reveal  the  coefficient  kR  to  equal  unity  for  prac 
ical  purposes,  meaning  the  propeller  slipstream  has  nearly  reached 
its  maximum  value  of  TJA+VA  at  the  rudder.  Thus  the  mean  flow  £,from 
eq.  ( 1 2 )  and  (13),  simplifies  as 

c2  =  V2  +  -^2  T  (14) 

71  pD 

-  v2  +  c2  T 
VA  +  °T  1 


The  combined  performance  of  screw,  rudder  and  engine  is  com¬ 
pleted  when  the  dynamics  of  the  engine  is  included. 


Engine  Equation.  From  eq.  (5)  it  is  apparent  that  the  engine 
model  required  must  generate  the  screw  angular  velocity  according 
to  some  function  of  command  inputs  and  in  turn  generate  a  torque  to 
compensate  for  the  propeller  load.  The  model  has  to  be  accurate  on¬ 
ly  to  the  extent  it  copies  dynamic  response  when  subject  to  load- 
variations  and  command  settings.  Formulation  in  terms  of  energy  ba¬ 
lance  then  is  adequate  for  adaption  to  the  different  types  of  en¬ 
gines  and  associated  control  mechanisms. 

Power  inlet  to  the  motor  Pm  is  proportional  to  fuel  flow,  and 
power  outlet  equals  developed  torque  Qm  times  screw  angular  velocity 
n.  Further  any  difference  between  developed  torque  Qm  and  load  torque 
will  change  the  timerate  of  screw  revolutions.  The  basic  equations 
thus  are 

V1  =  ^Pm  0  5) 

Imn  =  nQm  *  nQl 


where  Im  is  the  total  inertia  of  rotating  parts  plus  added  inertia 
from  the  water  as  referred  to  the  shaft.  n  is  the  efficiency  of  the 
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motor. 

The  build  up  of  developed  torque  from  power  inlet  reflects  the 
dynamics  of  the  motor.  In  eq.  15  this  is  taken  as  instantaneous  be¬ 
cause  the  engine  response  is  very  fast  compared  to  the  remainder  dy¬ 
namics  of  the  ship. 

The  governor  control  loop  contributes  by  its  low  frequency  per¬ 
formance  which  determines  the  stiffness  of  RPM  control  to  load  vari¬ 
ations.  In  fig.  4  a  block  diagram  of  the  general  engine  model  is 
shown. 


The  Speed  Loss  Equation.  From  the  preceeding  sections  the  dyna- 
mical  speed  loss  equation  is  formulated  for  a  vessel  without  exter¬ 
nal  disturbances  acting  upon  it. 


(““V*  =  (VmxG)i,+  (Yur-m)ru+Yuvuv+Ycc6C 


-2 


=  (N«--mxG^+(Nur-'nxG^ru+Nuvuv+Ncc6C^5 
m,i  =  xuuu2+(xvr+m)vr+(xrr+mxG)r2+xcc66E2tl2+(l-t)T 

V  =  VQ1 


T  =  T  n' +T  nV, 
nn  nv  A 


(16) 


Q,  =  Q  n  +Q  nV, 
1  nn  nv  A 


VA  =  (l-wo)u-wv(v-lpGr)u 

=  VI  +  4  T 


Hn  =  H 


The  ir  and  ft  relations  show  the  well  known  dynamical  equations 
of  motion  and  those  of  0  and  ft  constitute  the  propulsion  dynamics. 
The  remaining  expressions  form  the  algebraic  bounds  brought  up  by 
screw  and  wake  behaviour. 

The  upper  two  sway  and  yaw  rate  equations  in  (16)  have  the  rud¬ 
der  angle  as  driving  element.  Speed  and  torque  dependancy  is  affect¬ 
ing  the  coefficients  only.  Therefore,  this  part  describing  the  hull 
dynamics  is  decoupled  from  the  speed  loss  equations  for  the  small 
excursions  from  steady  advance  considered.  The  sway  velocity 

and  the  rate  of  turn  are  then  literally  inputs  to  the  speed  equat¬ 
ion  although  they  are  generated  by  antecedent  rudder  angle  movements 
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External  Disturbances.  Excitations  on  a  vessel  from  waves  and 
wind  are  complicated  functions  of  the  shape  of  hull  and  superstruct¬ 
ure  and  of  the  angle  of  attack.  When  forces  in  the  motional  equation 
vary  with  the  ships  heading,  the  course  becomes  necessarily  a  state  in 
the  equations,  and  the  external  disturbances  then  establish  feed¬ 
back  paths  from  the  course  to  yaw  and  sway  accelerations.  These  tim§ 
varying  coefficients  have  not  been  included  in  this  exposition  be¬ 
cause  they  do  not  affect  the  speed  loss  equation.  Further  the  full 
scale  trials  were  all  carried  out  in  nice  weather. 


When  autopilot  performance  is  considered  a  model  is  adapted, 
which  utilizes  a  set  of  transf erfunctions  relating  disturbance  po¬ 
wer  density  with  power  density  of  the  resulting  movement. 


Fig. 5  Wave  spectrum  as  seen  from  ship 
at  different  angles  of  attack. 


The  excitation  then 
follows  from  the 
spectrum  of  the  actu¬ 
al  disturbance  pas¬ 
sing  the  complex  trans- 
ferfunction.  The 
spectrum  of  wind 
pressure  is  approx¬ 
imately  white  whereas 
the  wave  spectrum  has 
bandpass  nature.  When 
the  ship  meets  the 
disturbances  at  steady 
advance  velocity,  the 
wavespectrum  is  trans¬ 
formed  by  a  nonlinear 
convolution  due  to  the 
natural  dispersion  of 
ocean  waves.  In  fol¬ 
lowing  sea  this  ef¬ 
fect  can  concentrate 
major  wave  energy  in 
a  tiny  bandwidth  at  a 
fairly  low  frequency 
of  encounter.  This 
situation  calls  for 
specific  attention 
when  basic  autopilot 
performance  is  inve¬ 
stigated.  A  plot  of 
the  encounter  spect¬ 
rum  is  shown  in  fig. 5. 


Summary  of  Ship  Dynamics. 

In  this  chapter  a  model  of  a  vessel  has  been  outlined  which  in¬ 
cludes  the  effects  of  propeller  interaction  with  hull  movements 
through  dynamical  variations  of  the  wake.  A  simple  input  -  output 
model  of  the  engine  completed  the  description  as  to  form  the  pre¬ 
sent  model  of  the  speed  loss  dynamics. 
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TECHNIQUES  OF  SYSTEM  IDENTIFICATION. 

Application  of  modern  identification  techniques  has  been  a 
fast  growing  field  of  practical  exercise  ever  since  efficient  mathe¬ 
matical  methods  and  proven  computer  algorithms  were  developed  some 
years  ago.  Methods  of  least  squares  fitting  to  curves  in  the  phase 
plane  and  different  methods  of  frequency  and  time  response  analysis 
has  a  tradition  of  long  practical  experience .when  deterministic  sig¬ 
nals  are  involved.  The  novel  methods  in  contrast  offer  the  possibi¬ 
lity  of  quantifying  stochastic  disturbances  also  in  a  mathematical 
framework  which  is  well  posed  for  applications  of  modern  control 
theory. 

At  the  time  of  planning  of  the  experiments,  however,  maximum 
likelihood  methods  had  not  proven  efficient  when  nonlinear  systems 
were  involved  and  were  still  at  a  state  of  active  research.  Actually, 
the  first  applications  to  nonlinear  ship  steering  dynamics  were  re¬ 
cently  reported  [11]. 

Instead  frequency  analysis  techniques  were  considered.  They  re¬ 
quire  the  system  to  be  excitated  by  a  signal  having  its  major  energy 
at  one  or  more  distinct  frequencies.  In  the  case  of  the  nonlinear 
speed  equation, excitation  by  a  single  sinewave  in  the  rudder  angle 
or  rate  of  turn  holds  a  unique  achievement.  A  linear  term  in  the 
equation  will  respond  at  the  excitation  frequency  whereas  a  square 
term  will  produce  a  response  containing  double  the  excitation  fre¬ 
quency.  Thereby  different  nonlinearities  appear  as  transparent 
effects.  Furthermore,  methods  of  spectral  analysis  are  well  fitted 
to  depict  the  very  small  amplitude  loss  terms  of  interest. 

The  frequency  analysis  method  does  not  allow  for  the  sophistic¬ 
ated  automatic  validation  of  model  structure  as  does  the  system 
identification  approach.  Instead  the  physical  comprehension  as  re¬ 
flected  in  eq.  (16)  has  to  be  relied  on  and  a  more  laborious  effort 
is  required  to  validate  the  model . 


Reformulation  of  Speed  Loss  Equation.  The  structure  of  the 
speed  equation  must  be  changed  from  eq.  (T6)  to  an  explicit  input 
output  form  when  used  for  identification  purposes. 

The  sway  and  rate  of  turn  equations  combine  to  form  the  two 
transf erf unctions  below,  both  having  rudder  angle  as  input. 


1+str 

cr  (1+st1 )(1+st2) 


and 


(17) 


Cv  ( 1+ST1 ) ( 1 +ST2 ) 


(18) 


The  timeconstants  and  gain  factors  are  easily  derived  from  eq. 
(16)  and  has  been  given  in  the  literature. 

The  speed  loss  equation  is  mostly  truncated  to  include  the  rv 
and  the  62  loss  terms  only  because  RPM  control  is  assumed  to  be  ideal 
In  actual  vessels,  the  governor  controller  has  a  certain  stiffnes  to 


A 
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load  variations  and  variations  in  RPM  are  induced.  Combined 

with  the  antisymmetric  wake  variation  with  sway  velocity  at  the 
stern,  this  can  contribute  significantly  to  propulsion  losses.  This 
is  seen  by  expanding  the  speed  equation  from  steady  state. 

By  elaborating  eq.  (16)  using  the  notation  U0 ,  n0,Wg  and  V0 
for  the  steady  state  values  of  ship  speed,  screw  revolutions,  wake 
and  inflow  velocity  at  the  propeller,  the  set  of  small  signal  re¬ 
lations  eq.  (19)  are  obtained.  The  symbols  u  and  n  now  denote  the 
deviations  of  speed  and  screw  rotation  from  the  steady  state  va¬ 
lues.  aQ„  is  differential  torque  and  G(n(t))  denotes  governor  respon- 
S0  •  ^ 

=  2XuuV+<Xvr+m>vr+Xcc6  6  d-w0)2U^2+(l-t)Tnnn2 

+ ( 1 -t ) ( 2Tnnn0n+TnvnQ ( 1 -wQ ) u+TnvU0wvn( v+lpQr ) )  (19) 

Xmfl  =  ^Sn-2Qnnnon+QnvnoUowv(v-1PGr) 

-Qnvno(l-wo)u-QnvU0(1'wo)n 
AO^  =  G(n(t) ) 


A  block  diagram  showing  the  structure  of  the  small  signal  speed 
loss  equation  is  shown  in  fig.  6  .  Propeller  revolution  is  seen  to  be 
excitated  by  the  stern  sway  velocity  and  the  resulting  variation  in 
RPM  gives  a  net  drag  force  proportional  to  n(v-lPGr).  The  magnitude 
of  this  term  is  discussed  below. 

By  neglecting  the  second  order  influence  on  speed  and  RPM  from 
the  propulsion  loss,  the  transferfunctions  listed  in  eq.(20)  result. 
Expressions  for  timeconstants  and  gain  factors  in  eq.  (20)  are  listed 
in  table  1 . 


(v+lpQr)  “  VSn  (st^+1  )(sts+1  )-kmkgaia2+(sTs+1  )G(s)km  ^20) 


eq. 


u  _  kmG(s)4-(sTm-H) 

^Toss  '  s  (sT“H')(STa+l)-\kBaia2+rsT8+l)fl(s)^; 


Total  drag 
21  . 

T  = 

loss 


force  from  the  loss  terms  is  finally  summarized  in 


(xvr+m)vr+xCC55(l-wo)2U262  + 

(1-t)(Tnnn2+TnvUoWvn(v+lpGr)) 


(21  ) 


In  the  transferfunctions  eq.  (20)  G(s)  is  the  transf erfunction  of 
the  governor,  from  RPM  deviation  from  its  setpoint  to  developed  mo¬ 
tor  torque.  This  controller  normally  includes  a  pure  integration, and 
the  two  transferfunctions  are  generally  of  third  order. 
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Term 

Expression 

Magnitude 

Unit 

kv 

vVQnvnoUo 

-4- 1 05 

Nns/m 

km 

’/t2fln nWo^1-*)) 

3 • 1 0~6 

rad/Nns 

k 

s 

5-10-6 

m/Ns 

ai 

-Qnvno^-wo^ 

1  .3-105 

Nns/m 

a2 

(l-t).2Tnvn0(1-w0) 

-2.3* 105 

Ns /rad 

Tm 

V<2Bm.Wo^o)) 

1  -5 

sec 

Ts 

-n’/(Vo+(1-t)Tnvno(,-vo)) 

100 

sec 

Table  1 .  Expressions  for  timeconstants  and  gain  factors  in 
the  speed  equation. 


Term 

Magnitude 

Unit 

Term 

Magnitude 

Unit 

VA 

9.0 

m/s 

Tnnno 

2.5'106 

N 

no 

10.0 

rad/s 

TnvnoVA 

-1  .9-106 

N 

T 

nn 

2 . 1 • 1 04 

Ns2/rad2 

T 

0.6-106 

N 

Tnv 

-1 . 8 • 1 04 

p 

Ns  /m  -  rad 

O  n2 
nn  o 

2 . 5 ' 1 06 

Nn 

Qnn 

2 . 1 • 1 04 

Nns2/rad2 

QnvnoVA 

-1  .6-106 

Nn 

Qnv 

-1 . 6 • 1 04 

Nns  /rad  m 

Q 

0.9-106 

Nn 

VUo 

0. 4-3.0 

xuu/ra 

-0.2'10"3 

m-1 

wo 

0.20 

xvr/m 

0.5 

(1-t) 

0.8 

xcc65Va/T 

-1  .2 

m/s^  Tad 

Table  2.  Typical  set  of  data  for  a  cargo  ship. 
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Identif icab ility  and  Techniques  of  Analysis.  Provided  measure¬ 
ments  of  sway  velocity  and  rate  of  turn  are  available  the  transfer- 
function  n  /(  v-lpQ  r)  can  be  derived  from  measurements  in  terns  of  val¬ 
ues  of  poles  zeros  and  static  gain  factor.  These  are  composite  ex¬ 
pressions  of  the  hydrodynamic  coefficients  and  a  simple  analysis 
shows  that  they  cannot  individually  be  obtained.  The  slope  of  the 
wake  variation  is  not  obtainable  either.  By  combining  excitations 
xrom  engine  torque  with  disturbances  in  v  and  r  from  rudder  movements 
however,  the  effective  coupling  kv  from  sway  velocity  at  stern  to 
load  torque  variation  can  be  determined. 

The  drag  force  from  losses  must  be  computed  from  the  different 
squareterms  and  crossproducts  62  ,  n2,  vr  and  n(v-lBrr).  When  hull 
movements  are  brought  about  solely  by  rudder  fluctuations,  v  and  r 
relate  to  6  through  eqs.  (17)  and  (18).  As  further  n  relates  to  v 
and,X  throuSh  (20)  it  is  impossible  to  identify  the  loss  term 
coefficients  from  eqs.  (21)  and  (20)  directly.  By  using  uncorrelat¬ 
ed  excitations  from  rudder  and  engine  torque  instead,  the  coeffici¬ 
ents  are  in  theory  identificable  when  the  dynamic  relations  between 
the  variables  are  also  taken  into  account. 


In  the  actual  experiments  only  rudder  movements  were  used  to 
create  hull  motions  because  of  hardware  constraints.  Numerical  valu¬ 
es  of  the  individual  loss  term  coefficients  are  therefore  not  at¬ 
tained  but  the  magnitude  of  speed  loss  generation  can  be  verified. 


When  the  excitation  signal  is  a  single  frequency  sinewave,  drag 
iorce  xrom  the  loss  appears  as  second  harmonic  of  the  input  frequen¬ 
cy  in  the  speed  recording.  Instantaneously  thrus  :  oscillations  at 
the  input  frequency  occur  from  the  variations  in  wake.  Because  of 
the  tight  coupling  of  the  system,  according  to  fig. 6  ,  a  calibration 
or  the  speed  reading  is  thus  provided. 


Fig. 6.  Small  signal  speed  loss  equation. 
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ABSTRACT 

This  paper  discusses  some  of  the  approaches  that  have  been  considered  for  use  in 
the  control  of  large  SES,  the  general  arrangement  of  the  ship  control  system,  the  Ride 
Control  System  (RCS),  and  the  Maneuvering  System. 

The  algorithms  considered  for  the  more  critical  functions  of  maneuvering  control 
and  heave  control  are  discussed  and  the  current  treatment  of  the  control  allocation 
law  is  presented. 

INTRODUCTION 

The  majority  of  advanced  ship  designs  place  at  least  some  emphasize  on  high  speed 
operation  and  reduction  of  manning.  This  trend  has  resulted  in  technology  programs  to 
provide  a  greater  level  of  automation  in  all  aspects  of  ship  control.  In  addition, 
the  marked  increases  in  availability  and  reliability  of  data  processing  equipment  have 
made  the  application  of  this  equipment  highly  competitive  with  more  traditional 
techniques.  The  design  of  a  large  SES  thus  provides  a  timely  opportunity  for  ap¬ 
plication  of  some  of  this  new  technology. 


Figure  1.  The  general  configuration  of  a  large  SES. 
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CONCLUSION 


From  the  mathematical  description  of  ship  forward  speed  dynamics 
it  was  found  that  variations  in  wake  fraction  due  to  varying  direct¬ 
ion  of  inflow  to  the  propeller  cause  fluctuations  in  RPM,  and  also 
gives  rise  to  a  drag  force  component 

TnvUowvn(v+1PGr) 

From  the  full  scale  experiments  it  is  demonstrated  that  the  con¬ 
tribution  from  this  effect  can  be  significant  and  certainly  when  au¬ 
topilots  are  considered,  which  possess  very  low  steering  generated 
drag  force.  Coordination  of  RPM  control  with  course  control  there¬ 
fore  seems  advantageous. 

The  identification  of  the  nonlinear  speed  equation  by  means  of 
frequency  analysis  techniques  showed  successful  in  validating  the 
structure  of  the  mathematical  model  and  in  determining  a  set  of  para¬ 
meters  in  the  input-output  description  of  the  system.  The  values 
found  compare  reasonably  to  those  predicted.  Explicit  determination 
of  the  individual  loss  terms  was  not  possible  from  the  experimental 
data  for  two  reasons.  Only  excitations  from  the  rudder  were  used  in 
the  experiment, and  the  exploration  of  the  dynamics  of  one  decade  in 
frequency  is  hardly  enough  for  the  identification. 

Finally  hydrodynamic  effects  appeared  as  an  essential  element 
of  the  dynamic  wake  variation. 
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AUTOPILOT  PERFORMANCE 

The  performance  of  course  keeping  autopilots  regarding  the 
speed  loss  from  steering  activity  is  determined  through  its  response 
to  disturbances  from  waves  and  wind. 

The  excitation  of  the  hull  from  waves  can  be  calculated  from 
strip  theory  methods,  and  computer  programs  excist  that  calculate 
amplitude  and  relative  phase  of  the  disturbed  motions  expressed  in 
sway  velocity  and  yaw  rate.  They  add  to  v  and  r  generated  by  the  rud¬ 
der,  to  form  the  resulting  hull  motions.  The  frequency  spectra  of 
the  disturbances  have  a  form  similar  to  the  wave  spectra  in  fig.  5- 
In  following  sea  the  excitations  gather  at  a  relatively  low  frequency, 
which  is  normally  somewhat  higher  than  the  closed  loop  bandwidth  of 
the  system.  As  rudder  movements  have  little  effect  in  that  frequency 
range,  the  effect  of  wave  interaction  is  mainly  to  cause  rudder 
generated  speed  loss.  The  vr  generated  speed  loss  is  not  necessarily 
excessive  when  large  excursions  in  v  and  r  are  brought  about  by 
waves,  because  the  phasing  of  v  and  r  in  some  situations  of  follow¬ 
ing  sea  are  in  quadrature.  The  rudder  induced  drag  in  this  situation 
will  be  proportional  to  the  square  of  the  high  frequency  gain  of  the 
autopilot  and  considerable  savings  are  thus  gained  by  decreasing  the 
system  bandwidth. 

The  contribution  to  speed  loss  from  the  wake  generated  variat¬ 
ions  in  RPM  is  best  judged  from  an  example. 

Example .  Assume  that  a  0.2  deg/sec  rate  of  turn  and  0.5  m/s  sway 
velocity  is  generated  in  quadrature  by  waves.  Local  sway  velocity  at 
stern  then  is  approximately  of  0.6  m  amplitude  and  30°  phase.  From 
the  estimated  value  of  k  ,  this  gives  an  amplitude  in  RPM  of  0.3 
rad/sec.  Note  that  from  the  n/r  plot  the  governor  will  not  suppress 
the  variations  in  the  range  of  frequencies  relevant  to  following  sea 
excitations.  From  table  2  and  3,  the  drag  force  is  then  found  to  be 

Tnv'Vwv'n(v+1PGr)  =  8*1o3n 
as  total  thrust  is  0.6*10^  N,  this  is  1.3$. 

The  wake  induced  loss  is  thus  significant  when  optimal  perform¬ 
ing  autopilots  are  considered.  Coupling  of  the  RPM  controller  to  the 
course  keeping  autopilot  therefore  seems  to  be  advantageous  as  to 
attempt  proper  phasing  of  torque  from  the  engine  and  steering  activi¬ 
ty  in  order  to  minimize  propulsion  losses. 
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Instead,  from  the  r/6  function,  the  magnitude  of  x,  was  obtain¬ 
ed,  and  the  remaining  parameters  in  the  n/r  and  u/r  functions  were 
then  estimated,  assuming  the  value  of  x,  as  fixed.  This  did  give  con¬ 
sistent  results  that  show  good  agreement  with  the  theoretical  pre¬ 
dictions.  The  estimated  parameters  are  listed  in  table  4. 

The  accuracy  of  the  log  reading  is  of  special  concern.  By  using 
the  relation  between  u/r  and  n/r,  one  gets  a  prediction  of  the  gain 
in  the  u/r  function  of  102  m/rad,  based  on  the  measured  n/r  gain 
estimate.  The  estimation  of  u/r  gives  128  m/rad.  The  logreading  is 
therefore  relied  on  although  some  spread  in  the  measurements  are 
apparent.  The  outlier  in  the  u/r  plot  originates  from  an  experiment 
having  a  very  large  amplitude  in  yaw  rate.  Influence  from  the  boun¬ 
dary  layer  is  presumably  causing  the  misreading  in  this  case. 

On  the  plots  of  u/r  and  n/r  the  measured  phase  is  seen  to  have 
different  shape  than  that  expected  from  the  transf erfunctions  in 
table  3.  The  measured  phase  has  an  ever  decreasing  asymptote  instead 
of  approaching  a  limited  value  as  expected.  This  is  explained  by  the 
presence  of  a  timedelay  in  the  build  up  of  wake  from  rate  of  turn. 

A  computation  of  the  delay  gives  consistent  values  for  the  transfer- 
functions  u/r  and  n/r,  and  also  for  u/s2 .  In  the  speed  equation  the¬ 
refore  kymust  be  replaced  by  kv*e-ST.  This  result  is  not  quite  as 
significant  from  the  n/r  measurements  as  from  u/r  due  to  spread  of 
the  high  frequency  values.  It  is  accepted  due  to  the  structure  of  the 
system,  which  is  verified  by  the  measurements  anyway. 

Hydrodynamic  memory  effects  are  described  in  the  literature 
l3],  and  are  known  to  give  rise  to  timedelays.  The  magnitude  was  sur¬ 
prising,  however,  and  even  more  as  the  r/&  measurements  bear  no  in¬ 
dications  of  memory  effects. 

Validation  of  the  loss  term  equation  is  difficult  because  of  its 
complexity.  At  low  frequencies  the  rv,  nr  and  n2  terms  will  dominate, 
and  at  high  frequency  rudder  induced  loss  is  preponderant.  The  fre¬ 
quency  spacing  of  the  double  pole  in  1 /x,  to  the  pair  of  zeros,  which 
depend  on  the  magnitudes  of  x„  ,  *44  ’  etc.,  determine  the  relative 
contributions  of  the  losses.  In  fig.  12  an  characteristic  is  drawn, 
which  approximates  the  measured  values.  It  has  the  poles  according  to 
table  3  and  two  zeros  0.15  Hz.  The  position  of  the  zeros  indicates 
the  ratio  of  the  x„  and  x^  loss  terms  to  be  of  magnitude  30.  This  is 
not  too  far  from  values  obtained  from  model  tests,  but  the  result  is 
not  significant.  The  DC  value  of  3  knots/degree  of  rudder  is  also 
much  larger  than  expected  from  practical  experience. 


Results  of  measurements. 

It  has  been  demonstrated  that  the  use  of  frequency  analysis  is 
feasible  to  verify  the  structure  of  the  speed  loss  equation,  and  that 
values  of  comparable  magnitudes  to  those  anticipated  are  obtained. 

The  measurements  demonstrate  that  the  effects  of  dynamic  wake  vari¬ 
ations  on  the  full  scale  vessel  are  perceptible  to  RPM  and  shaft 
torque.  Further  the  presence  of  memory  effects  in  the  build  up  of 
wake  from  hull  motions  is  rendered. 

The  experiments  fail  to  reveal  magnitudes  of  the  different  loss 
terms,  but  values  of  entire  speed  reduction  from  drag  forces  are 
obtained.  The  speed  loss  due  to  wake  variation  is  assessed  by  using 
the  measured  value  of  kv  together  with  thrust  coefficients  of  the 
screw. 
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Data  Analysis. 


From  the  recordings,  various  powerspectra  and  crossspectra  have 
been  computed  to  obtain  the  transf erf unctions  of  interest.  The  spect¬ 
ral  analysis  was  carried  out  using  fast  fourier  transform  techniques 
together  with  application  of  data  smoothing  in  order  to  obtain  mini¬ 
mal  resulting  variance. 

For  the  problem  at  hand  simply  averaging  transforms  of  conse¬ 
cutive  timeslices  showed  advantageous. 

In  fig.  8  some  characteristic  spectra  are  plotted.  In  the  part¬ 
icular  run,  excitation  frequency  is  0.02  Hz.  The  spectra  of  rudder 
angle  and  rate  of  turn  exhibit  excellent  signal  to  noise  ratios  as 
also  apparent  from  the  plots  of  timehistory.  In  the  spectrum  of  the 
rudder  angle  some  3.  order  and  higher  harmoni  s  are  present.  This 
is  due  to  slow  rate  saturation  of  the  rudder  which  can  also  be  seen 
from  the  nearby  triangular  shape  on  the  recording  of  rudder  angle 
over  time . 

The  most  interesting  spectrum  is  that  of  the  speed.  It  clearly 
contains  significant  power  at  both  the  excitation  frequency  from 
the  wake  variation  and  at  double  the  excitation  frequency  due  to 
the  loss  terms.  Signal  to  noise  ratio  of  the  speed  measurement  is 
not  impressive  but  it  is  sufficient  to  obtain  significant  results 
when  computing  transf erf unctions . 

The  transf erfunction  from  rudder  angle  squared  to  second  har¬ 
monic  of  measured  speed  constitutes  the  input-output  model  of  the 
speed  loss, when  only  excitation  from  the  rudder  is  considered. 

Values  of  this  transf erfunction  are  obtained  from  the  cross-spectra 
between  6s  and  u  and  the  powerspectra  of  the  square  of  s.  The 
values  can  then  be  compared  with  the  theoretical  expectation  obtained 
from  the  previous  equations . 

By  direct  substitution  into  eq.  (20)  and  (21)  high  order  expres¬ 
sions  result  which  are  irrelevant  to  the  range  of  frequencies  of 
interest  .  Some  simplifications  are  possible,  however,  to  achieve 
reduced  order  relations.  They  are 

-  The  engine  timeconstant  x_<<  thus  x_  ~  0.0 

^  m  s  m 

-  The  governor  is  approximated  by  a  pure  integration  G(s)  ~  k  /s 

O 

-  The  timeconstants  x,  and  x2  are  considered  to  cancel  in  the 
r/6  transf erfunction  when  u/&2  is  computed. 

The  set  of  reduced  order  transf erfunctions  listed  in  table  3  are 
those  used  for  the  model  validation.  Note  that  the  n/r  expression 
contain  a  pure  differentiation  due  to  the  integrating  part  of  the 
governor. 

The  results  obtained  are  plotted  in  fig.  9  to  12.  Each  computed 
value  is  circled  by  an  0  and  the  0.95  confidence  limits  of  each  is 
marked  by  bars . 

It  was  first  attempted  to  estimate  the  transf erfunctions  from 
least  squares  fitting  to  each  set  of  data.  This  was  not  successful. 
The  complexity  of  the  transferfunctions  and  the  spread  of  measure¬ 
ments  gave  inconsistent  values  of  poles  and  zeros. 
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Fig.  7.  Recordings  from  oscillating  experiment. 
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THE  CLIFFORD  MAERSK  EXPERIMENT 

The  frequency  analysis  technique  described  in  the  previous  sect¬ 
ions  has  been  used  on  a  container  ship  Clifford  Masrsk  owned  by  the 
Danish  shipping  company  A.  P.  Mailer.  A  number  of  experiments  were 
conducted  where  the  vessel  was  oscillated  about  a  straight  course  by 
means  of  the  rudder.  The  wavy  course  was  brought  about  by  a  sinewave 
signal  added  to  the  course  error  voltage  inside  the  autopilot. 


Experimental  setup. 

The  Clifford  Maersk  is  equipped  with  a  diesel  engine  having  RPM 
controlled  by  a  governor.  For  achieving  information  on  values  of  the 
control  inputs  and  output  signals  required  for  the  identification, 
the  following  instrumentation  was  used: 

course  from  gyro  compass 
r  rate  of  turn  from  rate  gyro 

6  rudder  angle  from  potentiometer  at  the  rudder  stock 
n  screw  revolutions  from  tachometer 
u  from  pitot  tube  log  located  near  the  stern 

Q  not  directly  measured  but  an  index  setting  sensor  was 

mounted  and  the  shaft  torque  was  obtained  through  cali¬ 
bration. 

The  experiments  were  all  carried  out  in  calm  water  in  order  not 
to  complicate  the  data  analysis  from  external  disturbancies.  Measure¬ 
ments  were  made  at  various  frequencies  between  0.005  Hz  and  0.04  Hz, 
the  range  which  is  most  important  to  the  manoeuvring  dynamics  of  the 
particular  vessel.  The  range  was  chosen  from  practical  reasons 
also,  however.  At  lower  frequencies  run  time  of  each  experiment  ne¬ 
cessarily  becomes  excessive, and  above  0.04  Hz  the  rudder  will  not 
respond  reasonably  to  large  commands  due  to  slew  rate  saturation. 

A  sample  of  the  timehistory  of  one  experiment  is  plotted  in 
fig.  7. 


A  measurement  of  sway  velocity  v  was  not  available.  Instead  the 
similarity  between  r,  the  rate  of  turn  and  v  must  be  used  to  elimi¬ 
nate  v  from  the  transf erf unctions,  eq.  20.  When  the  hull  motion  is 
due  to  rudder  activity  eqs.  (17)  and  (18)  combine  to  express  v  from 
r  as 


1+STv  r 

1+STr  'r 


(24) 


The  sway 

v-lpGr  = 


velocity  at  the  stern  is  then  related  to  r  through 

(25) 


-^PG-  e? 


cyN  h +s ' <  <  1PCTr-cvTr/cr  > 1P0“cv/cr  J } 


1+ST, 


The  constants  in  eq.  25  are  all  positive  because  c  is  negative 
due  to  the  sign  convention.  In  the  next  section  the  gain  in  eq.  (25) 
is  referred  to  as  1*  and  the  numerator  timeconstant  is  abbreviated 

1  *  . 


A 
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Estimation  of  Transf erfunctions .  Digital  time  series  analysis 
by  means  of  frequency  domain  techniques  is  well  established,  and  va¬ 
rious  methods  excist  for  estimating  transf erfunctions  for  a  system 
from  recordings  of  inputs  and  outputs. 

If  the  powerspectrum  of  the  input  x  is  Gx(f),  where  f  is  the 
frequency,  and  the  transf erfunction  of  the  system  is  H(f),  then  the 
powerspectrum  of  the  output  y  will  be  Gy(f). 

Gy(f )  =  | H(f ) I  2  Gx(f)  (22) 

The  spectrum  of  the  cross  correlation  function  between  x(t)  and 
y(t)  is  Gxy  (f) . 

Gxy(f)  =  H(f)-Gx(f)  (23) 

The  powerspectra  are  real  valued  functions  and  contain  inform¬ 
ation  on  amplitude  only.  The  cross  spectrum  is  complex  valued,  and 
from  (23)  GXY (f)  is  seen  to  have  its  imaginary  value  equal  to  that 
of  the  system.  By  means  of  GXY (f)  and  Gx(f)  it  is  thus  possible  to 
obtain  both  amplitude  and  phase  of  the  transf erfunction  H(f).  The 
redundant  information  on  the  amplitude  |H(f)|  when  Gy(f)  is  also  a- 
vailable  is  used  to  establish  the  confidence  of  the  measurement.  In 
the  plots  of  transf erfunctions  in  the  next  chapter,  the  intervals  of 
confidence  level  0.95  are  indicated.  The  likelihood  is  95%  that  each 
value  is  within  the  intervals . 

Here  fast  fourier  transform  techniques  were  used  to  compute  pow¬ 
er  and  cross  spectral  densities , and  data  smoothing  methods  were  care¬ 
fully  tested  because  very  optimistic  confidence  intervals  may  result 
due  to  bias  effects. 


Summary  of  Identification  Methods. 

In  this  chapter  frequency  analysis  methods  have  been  discussed 
and  found  attractive  regarding  the  separation  of  the  linear  from  the 
nonlinear  contributions  to  speed  variation.  The  speed  equation  was 
expanded  from  a  steady  state  assuming  small  signal  deviations  and 
thereby  a  set  of  input-output  relations  was  obtained.  The  consecut¬ 
ive  discussion  of  identificability  concluded  that  the  full  set  of  co¬ 
efficients  can  not  be  determined  when  only  excitations  from  the  rud¬ 
der  are  used. 


Figure  2.  Relationship  of  the  main  reverser  tube 
assembly  to  the  inboard  and  outboard 
waterjet  steering  sleeves. 

A  large  SES,  originally  planned  for  2000  tons  but  currently  fixed  at  3000  tons, 
has  been  under  study  for  some  time.  A  variety  of  designs  have  been  considered  but 
most  versions  have  had  the  general  configuration  indicated  in  Figure  1.  The  ship  is 
supported  on  a  cushion  of  air  confined  within  the  plenum  defined  by  the  ship's  two 
side  walls,  its  bow  and  stern  seals,  and  the  sea  surface.  The  air  mass  for  the  plenum 
and  for  the  bow  and  stern  seals  is  supplied  from  gas  turbine  driven  centrifugal  fans 
and  is  distributed  via  air  ducts  positioned  across  the  width  of  the  ship.  Propulsive 
drive  is  supplied  by  four  stern  mounted  waterjets  whose  inlets  are  located  on  the 
lowermost  portion  of  the  respective  sidewalls.  The  waterjets  are  vectored  (Figure  2) 
by  steering  sleeves  mounted  at  the  waterjet  exits.  Reversing  forces  are  generated  by 
bringing  "U"  shaped  tubes  called  reversing  tubes  into  place  in  front  of  the  respective 
outboard  waterjets. 

The  SES  design  includes  provisions  for  both  local  and  remote  control  of  the  bow 
and  stern  seals,  the  lift  system  engines,  the  "reversers, "  and  the  propulsion  system, 
including  vectoring  of  the  waterjet  nozzles.  The  interconnection  network  for  the 
remote  control  is  a  distributed  data  system  of  the  form  indicated  in  Figure  3.  The 
system  consists  of  local  data  terminals  which  relay  data  between  two  control  stations, 
the  Ship  Control  Console  (SCC)  located  in  the  pilot  house  and  the  Central  Control 
Station  (CCS)  located  in  close  proximity  to  the  SCC  but  one  deck  below. 

The  functional  split  between  the  SCC  and  CCS  is  shown  schematically  in  Figure  4. 
The  various  operational  and  ship  plant  monitoring  signals  are  picked  up  by  the  closest 
data  terminal  (with  the  exception  that  certain  critical  signals  are  placed  on  separate 
data  terminals  to  preclude  total  loss  of  ship  capability  through  loss  of  a  single 
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SHIP  CONTROL  CONSOLE 


CENTRAL  CONTROL  STATION  - 


ELECTRONIC  EQUIPMENT 
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CENTRAL  CONTROL 
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PROPULSION/ LI  FT 
CONTROL  CONSOLE 


DAMAGE  CONTROL,  AUXILIARIES 
&  ELECTRIC  CONSOLE 


Figure  3.  Location  of  the  Ship  Control  Console,  Central  Control  Station,  and 
Data  Terminals. 


Figure  U.  Functional  block  diagram  of  major  ship  control  system  elements 
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data  terminal)  and  conditioned  and  multiplexed  prior  to  transmital  to  the  Control 
Electronics  rack  where  the  various  signals  are  relayed  to  the  SCC  and  CCS.  Although 

not  shown,  "hard  wired"  paths  are  provided  for  certain  critical  functions. 

The  SCC  provides  two  virtually  identical  stations  (Figures  5  and  6) .  The  console 
is  manned  by  a  "Ship  Control  Officer"  and  an  "Assistant  Ship  Control  Officer."  The 
"helm"  normally  resides  with  a  single  specific  location;  however,  each  station  has 

full  capability.  Only  the  Propulsion  and  Lift  Control  Panels  and  the  Ride  Control  and 

Autopilot  Panels  are  shared.  Although  the  navigation  and  collision  avoidance. system 
(NAVCAS)  is  treated  as  a  separate  functional  system,  multi-purpose  displays'  for 
collision  avoidance  data,  surface  search  radar,  and  "electronic  charts"  are  included 
at  each  station.  The  Commanding  Officer's  station  is  located  just  aft  of  the  Pro¬ 
pulsion  Control  Panel  so  that  all  controls  settings  and  displays  are  readily  visible. 


Figure  5.  The  Ship  Control  Console  in  plan  view. 


MANEUVERING 

The  3KSES  design,  of  which  the  above  is  probably  most  representative,  will  be 
unusual  in  that  it  is  neither  a  prototype  nor  a  pure  technology  ship;  rather,  it  is 
intended  to  prove  out  the  military  utility  of  the  SES  concept.  In  addition,  a  goal  of 
the  SES  design  is  to  maintain  manning  at  the  minimum  level  consistent  with  the  safety 
and  operational  requirements  of  an  ocean  going  ship.  The  above  design  attempts  to 
satisfy  both  these  goals  by  normally  operating  with  only  two  men  at  the  CCS  and  with 
the  previouly  mentioned  "Ship  Control  Officer"  and  "Assistant  Ship  Control  Officer"  at 
the  SCC.  As  is  apparent,  it  will  be  possible  for  the  functional  loading  of  the  latter 
two  operators  to  become  quite  high  under  certain  conditions.  The  proposed  solution 
to  this  problem  is  to  supply  a  simple  maneuver  oriented  control  system  for  control  of 
the  ship  in  manual  and  emergency  modes.  This  is  intended  to  allow  two  operators  to 
maintain  control  of  the  ship  even  in  high  loading  situations.  This  same  maneuver 
oriented  system  also  includes  a  computer  assisted  primary  mode  whose  function  is  to 
provide  for  automatic  sequencing  of  effectors  and  yaw  stabilization.  This  is  intended 
to  minimize  operator  loading  during  normal  operation  and  fatigue  over  a  normal  watch 
period.  The  functional  relationship  between  the  elements  of  the  Maneuvering  Control 
System  is  indicated  in  Figure  7.  Both  a  direct  link  and  a  computer  assisted  link  is 
supplied  to  all  ship  effectors  including  the  four  engines,  four  waterjet  steering 
sleeves,  and  the  port  and  starboard  thrust  reversers. 
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Automatic  sequencing  of  the  effectors  and  the  control  logic  is  supplied  by  the 
Control  Allocation  Law  which  resides  in  the  main  control  processor.  The  potential 
effectiveness  of  the  Control  Allocation  Law  is  based  upon  explicitly  incorporating  all 
ship  control  effectors  into  the  maneuvering  control  law.  This  is  possible  since  a 
multiplicity  of  effectors  can  be  used  simultaneously  to  control  the  SES.  This  feature 
can  be  used  to  advantage  in  a  number  of  ways  of  which  the  obvious  choice  is  to  use 
this  redundancy  to  increase  the  "availability*'  of  the  ship.  Since  it  is  possible  to 
maintain  performance  near  the  maximum  practical  capability  of  the  ship  by  automat¬ 
ically  sequencing  out  "failed  effectors,"  this  approach  provides  for  a  graceful 
degradation  of  ship  performance  without  causing  an  excessive  increase  in  complexity. 

In  addition,  the  gas  turbines  which  drive  the  water jet  thrusters  can  be  very  wasteful 
of  fuel  during  low  speed  maneuvering.  The  control  algorithm  can  be  used  to  control 
fuel  consumption  during  these  maneuvers  and,  of  equal  or  greater  importance,  to 
minimize  engine  cycling  as  well  as  cycling  of  other  control  effectors.  A  number  of 
tradeoffs  are  required  in  order  to  establish  the  logic  and  algorithms  for  these 
control  processes.  It  is  of  course  necessary  to  assure  ship  saftey.  This  is  acheived 
by  means  of  operating  procedures  and  careful  programing  of  limit  schedules.  The  other 
tradeoff  is  with  regard  to  the  sensitivity  or  priority  of  the  various  sub-algorithms: 
fuel  conservation,  effector  ware,  and  quality  of  control.  Priority  is  established  by 
appropriate  gain  constants  or  schedules  in  the  algorithms. 


Figure  6.  General  arrangement  of  the  Ship  Control  Station. 

MANEUVERING  ALGORITHM 

Active  control  of  the  ship  with  respect  to  propulsion,  steering,  and  reversing 
affects  three  degrees  of  freedom:  surge,  sway,  and  yaw.  When  the  control  axes  are 
close  to  the  principal  inertial  axe  ;,  ship  motion  about  these  axes  is  determined  by 
forces  and  moments  about  the  control  axes  with  minimal  cross-coupling.  From  the 
operators'  point  of  view  this  is  very  desireable  since  the  operator  can  easily 
identify  the  ship's  response  to  his  command.  In  particular,  a  separate  propulsive 
sideforce  control  is  of  great  advantage  for  low  speed  maneuvering  since  it  can  be  used 
to  facilitate  docking  operations,  low  sideslip  turns,  and  general  station  keeping.  At 
high  speed  a  pure  propulsive  sideforce  is  not  nearly  so  helpful  since  a  combination  of 
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propulsive  turning  moment  and  sideforce  exists  which  can  generate  the  desired  turn 
rate  for  small  sideslip  angles.  In  fact,  a  particular  turn  can  be  acheived  by  a 
variety  of  thrust  vectoring  and  thrust  differential  commands.  This  flexibility  in 
implementing  a  desired  maneuver,  as  will  be  discussed  later,  will  be  used  to  acheive 
additional  ship  capability. 

Figures  8  and  9,  which  depict  the  maneuvering  controls  for  the  SES,  divide  the 
problem  into  two  separate  control  processes:  a  high  speed  and  a  low  speed  maneuver. 

The  high  speed  case  includes  all  of  the  elements  of  the  low  speed  case  save  for  the 
sideforce  command. 

In  the  computer  assisted  mode  three  basic  commands  (e.,  j  =  A.  S,  w)  are  avail¬ 
able  to  the  ship  operator  (subscript  A  denotes  forces  and  Commands  associated  with 
forward  and  reverse  thrust,  subscript  S  denotes  factors  associated  with  sidethrust, 
subscript  to  denotes  factors  associated  with  the  turn  command).  These  inputs  from  the 
helm  and  thrust  level  in  turn  must  be  scaled  to  the  force  and  turn  rate  commands  by  a 
factor,  F  (j  *  A,  S,  10),  which  accounts  for  the  efficiency  of  that  element  as  a 
function  if  speed  (e.g..  the  available  thrust  as  a  function  of  speed)  and  by  the 
saturation,  y,  of  a  given  effector  when  presented  with  multi-axis  commands  (e.g.,  if  a 
command  input  lies  outside  the  operational  envelope  of  the  ship,  the  command  signal 
will  be  reduced  by  a  factor  y).  In  general  it  may  not  be  possible  to  schedule  all  y. 
In  this  instance  some  sensing  procedure  will  be  required. 


Figure  7.  The  maneuvering  control  system  and  ten  effectors. 
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Figure  8.  Functional  elements  of  the  low  speed  controller. 


Figure  9.  Functional  elements  of  the  high  speed  controller. 
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The  maneuvering  command  values  will  thus  be  of  the  form: 


Jc  =*  tjFjY»  j  -  A,  S,  to  (1) 

Now,  for  effective  allocation  of  control  effector  commands,  it  is  necessary  that 
an  adequate  mathematical  model  of  the  contribution  of  the  individual  control  effectors 
be  available.  Such  a  relationship  has  been  established,  for  example,  in  Reference 
(2);  however,  these  expressions  are  decidely  non-linear.  Linearization  of  these 
equations  about  a  given  operating  point  is  not  particulary  attractive  since  realistc 
maneuvers  are  not  in  general  characterized  by  small  deviations  from  a  reference  state 
and  since  such  a  procedure  would  imply  a  great  deal  of  scheduling  and  possible  con¬ 
current  storage  demands. 

Fortunately,  if  the  commarided  forces  and  turning  rate  are  written  in  terms  of  the 
control  effector  variables  N,  6,  and  X  (which  describe  engine  speed,  waterjet  nozzle 
angle,  and  thrust  reverser  position): 

jC  *  V*’  T)  (2) 

it  becomes  apparent  that  a  linear  relationship  can  be  established  between  the  rate  of 
change  of  propulsive  forces  (A^,  Sc>  and  moment  (Mc),  and  the  control  variable  rates 
(N,  <S,  X)  simply  by  differentiating  equation  (2). 

This  relationship  is  still  not  in  a  particularly  attractive  form  since  it  implies 
that  the  command  input  must  be  the  time  derivative  of  the  command  forces  and  moments. 
Fortunately  again,  this  result  can  be  acheived  either  by  direct  differentiation  or  by 
feedback  techniques.  The  latter  is  of  course  to  be  preferred  since  the  control  error 
can  be  maintained  small  by  comparing  the  command  input  to  the  desired  output.  This 
will  in  turn  enhance  the  accuracy  of  those  equations  used  in  the  mathematical  model. 
The  control  inputs  for  the  distribution  law  used  in  the  low  speed  case  (Figure  8)  are 
thus : 


jc  «  Gj  (jc  -  jp),  j  =  A,  S,  M  (3) 

where:  =  the  control  gains 

jp  =  predicted  propulsive  forces  and  moment  as  derived  from 
physical  modeling 

=  commanded  propulsive  forces  and  moment. 


As  indicated  in  Figure  9,  in  the  high  speed  case  (above  about  10  knots),  the 

derivative  of  the  sideforce  command,  S  ,  is  not  used.  Instead,  a  derived  expression 
is  used  for  the  total  calculated  sideforce  : 

TP 

STp  =  mVg  (1  +  6)  (4) 


where:  m  =  ship  total  mass 
Vg  =  ship  speed 
u)  =  yaw  rate 

8  =  sideslip  angle  rate. 


It  follows  that  for  small  sideslip  angles,  the  total  radial  turning  force  is 
essentially  the  sum  of  the  hydrodynamic  sideforce  and  the  propulsive  sideforce.  A 
particular  rate  of  turning  can  therefore  be  acheived  by  some  combination  of  these  two 
components,  i.e.,  some  combination  of  thrust  vectoring  and  thrust  differential.  The 
longitudinal  force  command  is  introduced  by  thrust  level  deflection  while  the  lateral 
motion  command,  in  the  form  of  a  turn  rate  command,  is  introduced  by  deflection  of  the 


(2) 

helm.  The  latter  signal  is  the  input  for  the  yaw  stabilization  loop  indicated  in 
Figure  10.  The  loop  provides  yaw  stabilization  by  generating  a  command  which  is 
proportional  to  the  yaw  rate  and  the  deviation  from  the  yaw  rate  command: 

Mc  -  -  K2  Sw  +  (Kj  +  |£)  (u>c  -  u)  (5) 

where:  *  e  F  y. 

L  w  u) 


Figure  10.  The  Yaw  Stabilization  Controller, 

Formulation  of  the  Optimal  Controller 

The  preceding  dicussion  has  developed  the  relationship  between  the  operator  input 
commands,  A^,  u,  and  the  available  control  effectors,  or  more  precisely,  between 

the  rate  commands  and  the  control  effector  command  variables  N^,  6.,  X  .  For  the  low 
speed  case  these  relationships  can  be  summarized  by: 


— 

- , — 

A 

N 

c 

s  1 

=*  s 

6 

,c| 

M 

X 

L  cJ 

_  _ 

where:  S  =  a  3X10  matrix  of  coefficients  which  are  a  function  of 

N^,  6^,  and  X ^ 
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The  values  of  the  command  rates,  j^,  (j  =  A,  S,  M) ,  are  obtained  by  comparing  the 

input  commands,  j  ,  to  predicted  forces  and  moments,  jp,  which  have  been  calculated 
from  the  detailed  ship  model.  Expression  (6)  thus  represents  a  system  of  three  equa¬ 
tions  in  the  ten  effector  commands  N,  6,  A.  Activation  of  either  reverser  requires  the 
centering  of  the  corresponding  water jet  nozzle  so  that  only  eight  independent  commands 
are  available. 


Mathematically  we  have  a  system  of  three  equations  with  eight  degrees  of  freedom. 
This  redundancy  can  be  used  to  great  advantage  to  limit  effector  activity  and  fuel 
consumption,  and  to  increase  the  availability  of  the  ship  control  system.  The  former 
isachievedby  using  a  cost  function''  of  the  form: 


8  =  Kf 


•  5/4 

N.N.  +  K 


i-1 


i  n 
i=l 


+  K . 
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i=l 


6. 2  + 


l  A 

i-M1 


(8) 


The  various  K.  are  weighting  factors  which  determine  the  relative  importance  of 
tha.activity  of  th^  different  classes  of  effectors.  The  first  term  (proportional  to 
N  )  is  a  derived  cost  function  for  the  fuel  consumption  of  a  gas  turbine.  The  next 
three  terms  minimize  activity  of  the  engines,  water jet  nozzles,  and  reversers  by 
penalizing  large  values  of  the  effector  rate  commands.  It  should  be  noted  that  such  a 
scheme  could  be  adapted  to  allow  for  preferred  modes  of  operation.  As  an  example, 
could  be  set  equal  to  zero  for  high  speed  operation  to  assure  minimum  engine  cycling 
during  normal  operation,  but  could  be  set  equal  to  one  for  low  speed  operation  to 
minimize  fuel  consumption  during  docking  or  station  keeping. 


Increased  availability  isachievedby  means  of  a  sequencing  system  which  allows 
the  ship  to  maneuver  at  or  near  peak  performance  despite  the  loss  of  one  or  more 
control  effectors.  One  method  of  realizing  this  goal  is  by  incorporating  a  series  of 
effector-enable-flags  which  are  either  0  or  1  depending  or  whether  the  corresponding 
effector  is  disabled  (due  to  either  failure  or  specific  mode  selection)  or  enabled. 


Formally  if  we  introduce  scaled  variables,  n  ,  and  £ .  and  matrix  H,  x,  c,  and 
p  defined  by: 


n  «  K  1/2  N  , 
ini’ 

«i  -  Kx1/2  XA 

H  =  K 
ij  n  j’ 

m 

CM 

tH 

H 

-  K  -1/2f 

Kd  £j’ 

j  -  5,  6,  7,  8 

.  K  -1/2, 

j  =  9,  10 

(9) 
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(10) 


equation  (6)  becomes: 
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and  the  cost  function  g  becomes: 

g  *  2p^x  +  x  ^  x  (12) 

where  T  denotes  the  transpose  of  the  respective  matrix.  The  controller  thus  minimizes 
the  cost  function,  g,  by  solving  for  the  optimal  set  of  effector  commands,  x,  and 
includes  those  constraints  necessary  to  account  for  failed  effectors. 

Although  somewhat  complex  in  description,  the  Control  Allocation  Law  thus  sup¬ 
plies  a  reasonably  simple  algorithm  for  sequencing  the  multiple  effectors  available 
for  ship  control.  The  major  benefits  of  this  approach  are  decreased  operator  loading 
and  enhanced  available  for  the  ship  maneuvering  system.  An  additional  feature  is  the 
capability  to  adjust  the  activity  of  the  individual  effectors  based  on  field  expe¬ 
rience. 


Figure  11.  The  variable  geometry  fans  are  an  integral  part 
of  the  lift  fans.  The  vent  valves  are  connected 
directly  to  the  main  plenum. 

RIDE  CONTROL  SYSTEM  (RCS) 

The  other  major  task  of  the  ship  control  system  it  to  provide  attenuation  of  the 
heave  motions  of  the  SES.  This  is  accomplished  by  the  RCS  which  controls  pressure  in 
the  SES  plenum  and  which  consequently  controls  the  heave  motions  of  the  SES.  Figure 
1  indicates  the  general  physical  relationship  of  the  lift  system  elements  and  seals; 
while  Figure  11  provides  a  schematic  of  the  lift  system  air  flow. 

Lift  system  fans  are  driven  by  gas  turbines  through  a  series  of  reduction  gears. 
The  fans  are  centrifugal  flow,  double  axial  inlet  fans  with  an  inlet  variable  geometry 
(Vg)  feature.  All  lift  fans  are  identical  but  differ  in  actual  installation  ori¬ 
entation.  The  air  distribution  system  transfers  air  from  the  fans  to  the  seals  and 
main  cushion  by  means  of  independent  duct  systems.  The  forward  fans  discharge  air 
directly  into  the  bow  seal  via  a  short  diffuser  while  the  center  fans  have  a  similar 
discharge  path  to  the  main  cushion.  The  flow  of  the  aft-most  fans  is  transmitted  to 
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the  stern  seal  by  a  slightly  longer  combination  of  round  and  square  ducts  followed  by 
another  set  of  short  diffusers.  The  lift  system  effectors  which  can  be  controlled  by 
the  RCS  are  VG  fan  actuators,  vent  valve  (VV)  actuators,  and  the  rotational  speed  of 
the  VG  fans.  The  RCS  establishes  the  operating  point  of  the  system,  i.e.,  establishes 
the  reference  or  bias  settings  of  the  actuators  and  fan  speed,  and,  provides  feedback 
as  necessary  to  regulate  cushion  pressure  and  to  control  the  crafts  vertical  plane 
motions  (heave,  and  potentially  pitch  and  roll)  based  upon  the  sensing  of  either  heave 
acceleration  or  plenum  pressure. 

The  VG  fans,  as  the  name  implies,  are  fans  whose  general  conf iguration  or  geom¬ 
etry  can  be  varied  in  order  to  extend  the  useful  operating  range  of  the  fans.  In  the 
case  of  radial  fans,  the  angle  of  attack  of  the  blades  is  typically  varied,  while 
in  the  case  of  centrifugal  fans,  the  effective  vane  area  is  varied.  The  variable 
geometry  feature  of  the  centrifugal  fans  consists  of  two  translating  sleeves,  one  at 
each  inlet,  mounted  concentrically  around  the  axis  of  rotation  of  the  fans.  Movement 
of  the  sleeves  along  this  axis  effectively  shrouds  the  inlet  and  varies  the  effective 
inlet  area.  The  sleeves  can  be  positioned  at  a  fixed  point  or  modulated  continuously. 
The  steady  state  feature  is  utilized  to  provide  adequate  stall  margin  for  steady 
state  operation,  while  the  modulation  capability  is  used  during  dynamic  ride  control. 
Total  travel  capability  of  the  sleeves  in  the  latter  mode  is  about  17  inches  with  a 
half  power  point  corresponding  to  about  four  Hertz.  Typical  cushion  pressure  for  a 
3000  ton  SES  will  range  from  200  to  330  PSF  depending  on  operating  conditions. 

Nominal  corresponding  air  flow  will  vary  from  about  32,000  to  38,000  CFS  for  a  A0  to 
90  knot  range  of  ship  speed,  with  peak  flows  to  about  62,000  CFS.  The  W  are  basically 
openings  or  apertures  in  the  deck.  Their  area  is  modulated  by  rotating  vanes  or 
shutters  mounted  in  the  deck  openings.  Typical  cushion  pressure  for  a  2000  ton  SES 
would  range  from  approximately  180  to  290  PSF.  The  nominal  corresponding  vent  valve 
flow  would  range  from  near  zero  to  a  maximum  of  40,000  CFS  at  an  effective  venting 
area  of  144  square  feet. 


Figure  12.  Lift  and  Ride  Controller  functional  block  diagram. 
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The  RCS  can  operate  with  either  VG  fans  or  W  or  both.  The  configuration  now 
considered  to  have  the  best  performance  consists  of  VG  fans  and  W  operating  in  a 
push-pull  mode.  The  VG  actuators  are  biased  partially  "open"  and  serve  to  provide 
the  principal  control.  The  W  actuators  are  biased  "shut”  and  only.ppen  when  the 
venting  requirements  exceed  a  specified  level.  This  configuration^  has  been  shown 
to  be  the  most  efficient  in  terms  of  fuel  consumption.  The  VG  fans  are  basically  more 
fuel  efficient  than  the  W;  but  the  inclusion  of  the  VV  in  this  configuration  provides 
a  wider  range  of  control. 


COMMAND  TO 
ALL  6  VG 
ACTUATORS  (%) 


VG  BIAS 
SETTING  (%) 


Figure  13.  Block  diagram  of  a  Ride  Controller  using  acceleration  feedback. 

Stabilization 

Figure  12  depicts  the  functional  relationship  between  the  elements  of  the  lift 
and  ride  control  system.  A  signal  derived  from  either  the  ship  heave  acceleration  or 
cushion  pressure  is  used  to  provide  feedback  for  the  Ride  Controller.  The  sensed 
signal  is  conditioned  and  fed  to  a  controller  which  consists  of  a  heave  mode  compen¬ 
sation  network  and  a  stabilization  filter  (see  Figure  13  where  for  simplicity  only  the 
VG  fan  control  is  indicated).  A  variety  of  networks  have  been  studied  for  ap¬ 
plication  in  the  stability  filter.  The  basic  choices  are  (1)  filtering  via  vector 
cancellation,  (2)  filtering  via  sensor  location,  (3)  linear  filtering,  and  (4)  non¬ 
linear  filtering.  Filtering  via  vector  cancellation  involves  separate  filtering  of 
each  transducer  signal  before  averaging  or  summing  of  the  total  signal.  This  approach 
foLlows  from  the  fact  that  for  frequencies  at  or  below  the  basic  heave  mode  frequency 
region,  pressure  components  respond  in  phase  in  both  the  bow  and  stern  seals  and 
throughout  the  plenum.  At  higher  frequencies,  phase  differences  develop  between  the 
various  pressure  elements.  It  is  thus  possible  in  principal  to  phase  the  signals  from 
two  or  more  transducers  in  such  a  manner  as  to  cancel  each  other  in  specific  frequency 
regions.  This  method  would  allow  filters  to  be  rolled  off  at  higher  frequencies  than 
would  be  the  case  with  linear  filtering  alone,  thus  allowing  less  phase  lag  and  con¬ 
sequently  providing  better  control.  In  practice  this  approach  has  not  been  widely 
applied. 

From  the  previous  dicussion,  it  is  evident  that  the  phase  and  magnitude  of  the 
transducer  signal  are  location  sensitive  and  thus  provide  the  potential  for  stability 
control.  Again  this  approach  has  had  little  practical  application.  Application  of 
linear  filters  to  the  stabilization  problem  are  reasonably  straight  forward  except  in 
those  instances  when  structural  or  pneumatic  modes  frequencies  fall  in  the  close 
proximity  to  the  basic  heave  resonance  of  the  ship.  When  such  an  instance  occurs  a 
very  elaborate  linear  (higher  order)  filter  may  be*  required  to  separate  the  modes. 

Still  another  approach  is  the^use  of  a  non-linear  filter  of  which  the  Indepen¬ 
dent-Magnitude-Phase  (IMP)  filter'  '  is  one  example.  The  IMP  filter  has  been  con¬ 
sidered  because  it  is  capable  of  providing  both  independent  magnitude  and  phase  as 
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indicated  in  Figure  14.  The  signal  is  first  separated  into  two  channels.  In  the 
first  channel  the  signal  passes  through  a  linear  filter,  F  (s) ,  whose  output  is  then 
passed  through  a  bistable  element.  In  the  second  channel  the  signal  is  passed  through 
another  Linear  filter,  F?(s)  and  the  absolute  value  of  the  filtered  signal  is  then 
taken.  Finally  the  outputs  of  the  separate  channels  are  recombined  by  multiplication. 
As  configured  in  Figure  14,  networks  F.(s)  provides  phase  lead  while  F-Cs)  provides 
signal  attenuation.  The  typical  resultant  transfer  function  is  indicated  in  Figure 
15.  The  filter  obviously  achieves  the  desired  result  of  good  high  frequency  isolation 
with  minimal  lag;  of  course,  careful  attention  to  real  world  considerations  such  as 
noise  background  and  quantization  effects  is  required. 


ABS.  VALUE 


BISTABLE 

FUNCTION 


Figure  14.  Block  diagram  of  an  Independent 

Magnitude  and  Phase  (IMP)  filter. 


Figure 


Heave  Mode  Compensation 


Early  in  the  development  of  the  SES,  a  mathematical  model  for  predicting  the 
loads  and  motions  of  the  SES  was  established.  Motion  data  from  this  modeling  was  in 
turn  used  in  man-machine  simulations''  which  established  the  desirability  of  a  RCS 
for  SES  operation  in  high  sea  spates.  Work  (see  Reference  (9)  for  example)  has 
continued  on  both  the  modeling  and  the  RCS  design  since  that  time  so  that  both  five 
Degree-of-Freedom  (DOF)  frequency  domain  and  six  DOF  time  domain  simulations  a^g. 
currently  available.  The  resulting  equations  of  motion,  even  when  linearized^ 
are  still  rather  formidable.  In  practical  application  the  design  of  the  RCS  has  been 
primarily  restricted  to  an  initial  simulation  of  two  DOF,  pitch  and  heave,  followed  by 
a  more  elaborate  checkout  simulations.  (This  has  been  true  for  two  reasons:  (1) 
current  designs  only  attempt  to  deal  with  heave  motion,  (2)  the  preponderance  of 
motion  is  associated  with  the  heave  and  pitch  DOF.)  Lags  for  control  actuators  and 
air  inertia  are  incorporated  as  are  non-linear  describing  functions  to  account  for 
saturation  of  the  flow  capability  of  the  vent  valves  and  the  VG  fans.  The  various 
gains  are  then  optimized  in  terms  of  this  limited  simulation  and  are  finally  tested 
for  overall  effectiveness  using  the  more  accurate  6  DOF  simulations.  Despite  its 
seeming  unwieldiness,  such  a  paramaterization  allows  a  reasonably  rapid  approach  to 
practical  gain  values. 

For  the  purpose  of  discussion,  it  is  still  possible  to  consider  only  the  two  DOF 
modeling.  Based  on  this  representation  and  including  only  first  order  effects  in  the 
frequency  region  of  the  basic  heave  motion,  the  transfer  function  for  the  acceleration 
of  the  center  of  gravity,  a  ,  relative  to  the  wave  pumping  of  the  cushion  volume, 

Vy ,  is  given  in  terms  of  the  heave  mode  natural  frequegry,  the  heave  mode  lag 

frequency,  and  the  heave  mode  static  gain,  K^,  by 
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cushion  gage  pressure  (PSF) 
cushion  pressure  (PSF) 
total  system  flow  rate  (CFS) 
fan  slope  (CFS/PSF) 

leakage  as  a  function  of  vehicle  and  seal  motion  (CFS/PSF) 
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cushion  area  (ft  ) 

total  vehicle  mass  (slugs) 

specific  heat  of  air 
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Figure  6.  Summary  Block  Diagram  of  Helmsman* s 
Control  Response. 


CONCLUSIONS 

An  examination  of  the  open- loop  steering  describing  function  measurements  in 
the  SES  simulation  has  shown  that,  regardless  of  sea  motion  condition,  both  helms¬ 
men  adopted  average  heading  and  displacement  gains  so  as  to  maintain  the  closed- 
loop  displacement  bandwidth  in  the  neighborhood  of  0.2  rad/sec  with  acceptable 
margins  of  stability  in  phase  and  gain.  This  is  an  adequate  bandwidth  for  course¬ 
keeping  at  SES  cruising  speeds.  Additional  adopted  first-order  lead  equalization 
appears  within  the  heading  loop  between  0.3  and  0.5  rad/sec  to  maintain  a  well- 
damped  characteristic  directional  mode  of  the  ship.  This  lead  equalization  may  be 
generated  by  perceiving  the  turn  rate  indicator. 

Although  no  evidence  for  biodynamic  amplification  appears  in  these  results, 
relatively  lightly  damped  and  amplified  steering  wheel  dynamics  appear  at  about 
1 .25  rad/sec  as  an  intended  artifact  of  the  simulation.  The  relatively  low  damping 
is  due  in  part  to  the  lower- than- recommended  values  of  Coulomb  friction  provided 
in  the  simulated  artificial  feel  for  the  steering  wheel  to  reduce  breakout  torque. 
However,  the  spring  gradient  in  the  steering  wheel  was  also  too  low.  Consequently, 
the  natural  frequency  of  neuromuscular  actuation  was  too  low,  and  the  helmsman 
reduced  the  inherent  steering  wheel  damping  ratio  even  more,  because  a  portion  of 
his  time  delay  exists  within  his  perceived  rudder  angle  feedback  loop. 

That  the  helmsmen  were  scanning  the  rudder  angle  indicator  is  corroborabed 
by  their  commentary.  This  is  presumably  a  part  of  the  transfer  of  training  fcr 
steering  low  speed  displacement- hull  ships.  However,  with  a  stiffer  and  more 
linearly  damped  artificial  feel  system  with  a  lower  breakout  force,  the  proprio¬ 
ceptive  feedback  available  from  the  steering  wheel  will  provide  a  superior 
equivalent  to  visual  rudder  angle  feedback  from  the  panel  instrument .  As  a  result, 
the  neuromuscular  dynamics  will  be  desirably  higher  in  frequency  and  suppressed 
in  amplitude.  This  will  yield  a  better  margin  of  stability  and  improved  course¬ 
following  performance  in  aggravated  seas. 

The  overall  effective  time  delay  of  the  helmsmen  estimated  from  the  open- loop 
describing  function  measurements  is  between  2  and  3  sec.  This  is  about  2  sec  larger 
than  would  be  expected  to  accompany  the  adopted  equalization  when  using  the  elec¬ 
tronic  horizontal  situation  display,  duch  a  relatively  large  time  delay  may  be 
caused  by  scanning  delays  among  the  helmsman’s  instruments  and  visual  field  and  by 
a  division  of  attention  among  his  other  tasks  during  the  measurement  interval.  We 
therefore  recommend  training  in  the  more  effective  use  of  integrated  horizontal 
situation  displays  for  steering. 
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c)  r>w  is  either  0.0ft  (Crew  A)  or  0.15  (Crew  B)  to  fit  the 

describing  function  amplitude  at  the  third  DFA  measure¬ 
ment  frequency.  These  relatively  low  values  of  are 
due  in  part  to  the  lower- than- recommended  values  of 
Coulomb  friction  provided  in  +he  simulated  artificial 
feel  for  steering  Ho;  ever,  they  also  repre¬ 

sent  the  fa-  hat  the  helmsman  will  reduce  the  inherent 
steering  wheex  damping  ratio  anyway,  because  a  portion 
of  his  time  delay  exists  within  his  perceived  rudder 
angle  feedback  loop  or  his  proprioceptive  steering  angle 
feedback  loop. 

d)  is  between  2  and  3  sec  to  fit  the  describing  function 
phase  angle  at  the  third  DFA  measurement  frequency. 

e)  Yy  is  a  pure  gain  and  the  gain  product  YyK^  is  adjusted 
so  that  the  amplitude  of  M/E  matches  that  at  the  lowest 
DFA  measurement  frequency. 

The  quite  satisfactory  results  of  the  fitting  procedure  for  the  helmsman  in  Crew  B 
are  displayed  in  Fig.  h  and  tabulated  in  Table  2  for  both  helmsmen,  and  a  summary 
block  diagram  equivalent  is  shown  in  Fig.  6. 


Table  2 .  Numerical 

Values  Representing  the 

Helmsmen. 

Crew  A 

Crew  B 

VYy  Cn) 

791- 

791. 

V*  (rad/ft) 

0.0000927 

0.00009S1 

V* 

0.0735 

0.077c 

Oijj  (rad/sec } 

1.25 

j  .25 

Tl  (sec) 

3. 

2.5 

"’N 

0.08 

0.15 

xd  (sec) 

2-9 

2. it 

coc  (rad/sec) 

0.18 

0.19 

phase  margin  (deg) 

29. 

28. 

<n  (rad/sec) 

0.5 

0.525 

Gain  margin  (dB) 

7. 

8. 

F2  2-10 


Figure  5.  Bode  Root  Locus  of  the  Zeros  of  [Y^sNsj.  +  Y  N^r]  as 
a  Function  of  the  Pure  Gain  Ratio  Y^/Yy  for  |y »  |y yl 
(Refer  to  the  Numerator  of  Eq.  3) 
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Numerical  estimates  for  the  parameters  in  F^,Yy  were  based  on  the  following  obser¬ 
vations  and  conditions: 


a)  Tl  is  approximately  2  or  3  sec  to  provide  a  phase  cross¬ 
over  frequency  in  the  neighborhood  of  the  second  DFA 
measurement  frequency,  0.5  rad/sec,  with  the  5  or  6  dB 
gain  margin  measured  and  to  maintain  a  well- damped  non- 
oscillatory  characteristic  directional  mode  of  the  ship 
emanating  from  A. 

b)  cjjji  is  slightly  in  excess  of  1  rad/sec  and  approximately 
equal  to  the  third  DFA  measurement  frequency,  1 .25  rad/ 
sec . 
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Substitution  in  Eq.  3  of  the  typical  numerical  values  for  the  ship's  lateral- 
directional  transfer  functions  in  Table  1  will  reveal  the  following  key  points: 


•  M/E  is  of  the  form  K/s2  at  low  frequency  as  we  would  expect. 

•  lY\>l  >:>  I  Yy- 1  in  order  that  sN^r  may  provide  the  requisite 
low-frequency  lead  equalization  to  convert  M/E  to  the  form 
K/s  in  the  region  of  unit  gain  crossover. 

•  The  characteristic  directional  oscillatory  mode  of  the  ship 
is  quite  low  in  frequency  (approximately  0-5  rad/sec)  but 
well  damped. 

•  Y*  (or  F,.,)  should  adopt  lead  equalization  in  the  vicinity 

of  0-5  rad/sec  to  maintain  a  well  damped  closed- loop  directional 
mode . 

•  The  characteristic  rolling  oscillatory  mode  is  higher  in 
frequency  ( approximately  1 . 5  rad/ sec ) ,  lightly  damped,  but 
suppressed  in  amplitude  by  the  zeros  of  the  lateral  and 
directional  numerators. 

We  shall  next  discuss  the  procedure  employed  to  establish  F^,  Y„,  and  Yy 
which  will  provide  plausible  heading  and  lateral  displacement  control  techniques 
to  interpret  the  response  measured  by  the  DFA.  We  shall  illustrate  the  procedure 
by  adopting  pure  gain  equalization  within  Y^.  and  deferring  the  mid- frequency  lead 
equalization  to  F* ,  although  equii  lent  results  can  be  obtained  by  adopting  the 
converse  equalization  technique  because  | Y^,  (  »  |  Yy  f .  In  either  case  pure  gain 
equalization  will  suffice  for  Yy. 

It  was  first  necessary  to  determine  the  behavior  of  the  numerator  zeros  of 
Eq.  5  as  a  function  of  Yijr/Yy  and  to  select  an  appropriately  large  value  for  the 
gain  ratio  Y^/Yy  which  would  fix  the  location  of  the  zeros  for  the  overall  lateral 
displacement  control  so  as  to  provide  the  low-frequency  lead  equalization  apparent 
in  the  measured  response.  The  Bode  root  locus  of  the  numerator  zeros  was  there¬ 
fore  generated  as  depicted  in  Fig.  5.  The  gain  ratio  of  791  was  selected  for 
Y^/Yy  to  provide  M/E  with  a  reasonable  frequency  interval  having  the  form  approach¬ 
ing  K/s  in  the  neighborhood  of  the  unit  gain  crossover  frequency,  0.1 9  rad/sec 
estimated  from  the  DFA. 

With  the  gain  ratio  of  Y^/Yy  thus  established,  M/E  is  now  represented  partially 
in  numerical  form  in  Eq.  U,  by  allowing  Yv  to  covary  with  Yy  in  accord  with  the 
constant  ratio  selected  above. 
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It  now  remains  to  determine  the  product  Fj,Yy  so  that  M/E  in  Eq.  4  will  fit  the 
set  of  three  measurements  in  Fig.  4.  Realizing  that  the  helm  itself  was  a  rela¬ 
tively  lightly  damped  spring-restrained  steering  wheel  with  a  very  low  undamped 
natural  frequency  uojj  on  the  order  of  1  rad/sec,  we  shall  hypothesize  that  the 
important  featui  3s  of  F^  can  be  represented  by: 
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Measurements  represent  the  ensemble  average  and  plus  or  minus  one  standard 
deviation  of  six  runs  by  the  helmsman  in  Crew  B.  The  theoretical  fitted 
describing  function  is  based  on  Eq.  3  and  represents  the  combination  of 
the  SES,  the  helmsman,  displays,  and  steering  control. 

Figure  h.  Example  of  Open- Loop  Describing  Function  Measurements  at 
the  Rudder  Control  Point  for  Manually  Controlled  Course-Keeping 
with  a  Simulated  2000T  SES  . 


loop  lateral  displacement  control  technique  from  rudder  error  to  helmsman's  output 
is  obtained  from  inspection  of  Fig.  1. 
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An  attempt  was  made  to  determine  the  influence,  if  any,  of  performing  the 
subcritical  speed  control  task  on  the  measurements  of  M/E  for  the  steering  tasks. 

No  influence  is  evident  in  the  two  runs  for  the  helmsman  in  Crew  B.  However,  for 
the  helmsman  in  Crew  A  only  one  run  with  both  tasks  is  available;  its  measurements 
appear  more  noisy  with  evidence  of  increased  time  delay  in  the  steering  task,  and 
the  phase  crossover  frequency  apparently  decreases  from  a  value  in  excess  of 
0.5  rad/sec  to  about  0.2  rad/sec  as  a  result  of  the  division  of  attention.  Yet, 
there  is  no  evidence  of  a  corresponding  gain  reduction  by  the  helmsman  in  Crew  A 
to  provide  more  than  2-5  dB  gain  margin  of  stability  when  both  tasks  are  being 
performed.  There  also  was  the  possibility  of  biodynamic  coupling  from  heaving 
motions  to  steering  motions  when  the  helmsman  had  only  his  right  hand  on  the  wheel 
while  performing  the  speed  control  task  with  his  left  hand.  Although  we  shall 
identify  neuromuscular  dynamic  amplification  in  the  results  at  the  third  measure¬ 
ment  frequency,  there  is  no  evidence  for  biodynamic  causation. 

The  mean  values  and  standard  deviation  of  M/E  for  six  selected  runs  by  the 
helmsman  in  Crew  B  are  presented  in  Fig.  k.  Four  are  runs  for  which  the  helmsman 
performed  only  the  steering  task,  and  two  are  runs  for  which  he  performed  both  the 
steering  and  speed  control  tasks  simultaneously.  The  amplitudes  and  phase  angles 
for  M/E  are  presented  only  at  the  three  lowest  measurement  frequencies  correspon¬ 
ding  to  0.188k,  0.501k,  and  1 .25C  rad/sec,  because  M/E  at  the  two  higher  measure¬ 
ment  frequencies  had  very  low  signal- to-noise  ratios  and  was  dominated  by  noise. 

The  average  unit  gain  crossover  frequency  of  the  lateral  displacement  loop  closure 
by  each  helmsman  is  slightly  below  0.2  rad/ sec,  and  the  average  phase  crossover 
frequency  is  in  the  vicinity  of  0.5  rad/sec. 

The  excessively  variable  and  noisy  measurements  obtained  here  are,  in  part, 
the  result  of  deliberately  violating  a  caveat  in  using  the  DFA,  via.,  maintaining 
input  magnitude  sufficient  to  yield  reasonable  displayed  error  deviations(°) .  This 
caveat  was  regretfully  sacrificed  in  favor  of  crew  motivation,  because  the  helmsmen 
complained  that  higher  input  amplitudes  produce.,  abnormally  great  activity  in  turn 
rate  and  rudder  displacement  based  on  their  experience,  even  though  the  input  dis¬ 
turbance  was  applied  to  the  innermost  loop  available  at  the  rudder  control  point 
and  not  directly  displayed  on  the  rudder  angle  indicator.  This  resulted  in  a  K/u>2 
power  spectrum  in  displayed  turn  rate,  whereas  one  should  in  future  tests  employ 
a  sum  of  approximately  equal  amplitude  sinusoids  for  the  input  at  the  rudder.  This 
will  allow  the  SES  controlled  element  to  shape  the  power  spectra  of  the  displayed 
signals  and  will  yield  better  signal- to-noise  ratios  in  the  measurements.  In 
retrospect,  we  placed  too  much  weight  on  crew  experience,  because  at  that  point  in 
time  neither  helmsman  had  had  experience  steering  an  SES  in  aggravated  seas  at 
cruising  speeds  by  reference  to  instruments. 

We  shall  now  describe  the  rationale  for  partitioning  m/e  between  the  (known) 
controlled  element  and  the  helmsman  in  such  a  way  as  to  infer  forms  for  his  de¬ 
scribing  functions.  The  results  of  applying  these  inferred  forms  within  a  theoreti¬ 
cal  model  for  M/E  and  fitting  the  same  to  the  measurements  are  plotted  in  Fig.  k. 

MANUAL  CLOSURE  OF  HEADING  AND  LATERAL  DISPLACEMENT  LOOPS  LEADING 
TO  AN  EXPLANATION  OF  THE  MEASURED  DESCRIBING  FUNCTIONS 

The  topology  of  the  helmsman's  control  technique  represented  in  Fig.  1  is 
founded  on  foreknowledge  of  the  K/s^  form  of  the  SIS  controlled  element  and  includes 
the  necessary  heading  and  lateral  displacement  loops  with  provision  for  lead 
equalization  in  the  heading  loop.  The  helmsman  is  represented  by  partitioned 
describing  functions  Fy ,  Y^,,  and  Yy.  The  describing  function,  m/e,  of  the  open 
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Figure  3.  Exemplary  Time  Histories  of  Surface  Effect  Ship 
Dynamic  Responses  to  a  Step  Change  in  Rudder  Angle. 


The  lags  and  non-minimum  phase  characteristics  combine  to  cause  apparent  ship 
response  delays  in  lateral  displacement  on  the  order  of  several  seconds,  which 
makes  for  a  challenging  dynamic  control  task.  The  various  displays  available  to 
the  helmsman  help  by  allowing  him  to  control  intermediate  states  or  derivatives 
of  the  final  lateral  displacement.  In  Fig.  1  the  Yy,  Y* ,  and  F^,  terms  account  for 
the  helmsman's  operation  on  the  variety  of  displayed  information. 

DESCRIBING  FUNCTION  MEASUREMENTS  WITH  MANUAL  CONTROL 


The  forcing  function  provided  by  the  DFA  is  labeled  "I"  in  Fig.  1,  where  it 
is  injected  into  the  common  path  of  the  control  loops  as  a  disturbance  by  summation 
with  5r,  the  helm's  steering  signal  to  the  rudder.  If  the  disturbance  I  is  viewed 
as  a  "command"  input  and  the  helmsman's  rudder  signal  Bj.  =  -M,  the  negative  motion 
feedback,  the  rudder  error  5e  =  E(=  I  -  M)  acts  to  disturb  heading  and  to  displace 
the  ship  laterally  from  the  desired  course.  Thus,  for  our  purpose,  E  represents 
an  acceptable  measure  of  closed- loop  performance,  the  finite  Fourier  transform  of 
which  is  computed  by  the  DFA,  and  [E/l](j<o)  represents  the  closed-loop  error-to- 
input  describing  function. 

The  open- loop  describing  function  of  the  helmsman’s  lateral  displacement  con¬ 
trol  technique  in  combination  with  the  known  mathematical  model  of  the  SES  can  be 
represented  by  [M/E](j<a),  which  may  be  obtained  from  the  relationship: 


[M/E]  ( jco) 


1  -  [E/lKjca) 
[E/lKjco) 


(2) 


After  identifying  that  portion  of  [M/E](ju>)  which  represents  the  known  .'.EG  dynamics, 
the  remaining  portion  will  represent  the  combined  equi valent  of  the  loops,  gains, 
frequency- dependent  equalization,  and  steering  dynamics  involved  in  the  helmsman' s 
control  technique.  A  digital  computer  program  based  on  Allen'  'J!  is  used  for  i  ienti- 
fying  [M/E](jcu)  and  partitioning  the  result  between  the  machine  'md  the  man. 

Some  of  the  results  of  applying  that  program  to  the  SES  simulation  have  been 
analyzed.  We  have  found  no  significant  differences  in  M/E  which,  correlate  with 
the  different  sea  motion  conditions  tested.  The  only  significant  liffere.ice  be¬ 
tween  helmsmen  is  in  the  average  amplitude  of  m/e  at  the  third  measurement  !><■  ;u eno,v. 
which  is  about  6  dB  higher  for  the  helmsman  in  Crew  A  than  in  Crew  li. 
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Table  1 .  Typical  Lateral-Directional  Open-Loop  Control  Response  Transfer 
Functions  for  a  2000T  SES  at  Cruising  Speed 


Denominator: 

A(s)  =  [s2  +  2(o.898)(o.it69)s  +  (o.W>9)2] 

x  [s2  +  2{0.198)(1  .55)s  +  (1.55)2] 

A(0)  =  O.52S  sec-1* 

Numerator  for  yaw  rate  response  to  rudder  angle: 

N^.(s)  =  0.999(s  +  0.493)[-2  +  2(0.1 91  )(l  .55)3  +  (1 
N^/0)  =  1  .185  sec"^ 

Numerator  for  lateral  acceleration  response  to  rudder  angle: 

N^r(s)  =  -52.o(s  +  o.903)(s  -  1.03) 

x  [s2  +  2(0.191  )(l  .51  )s  +  (1 .31  )2] 

Nfr(0)  =  120  ft/see6-rad 


Path 


Figure  2.  Helm  Control  Tasks. 
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multiloop  response  properties  with  a  single  disturbance  input  as  used  in  these  ex¬ 
periments  are  described  in  Teper(“)  and  McKuer ' ^ '  where  the  principles  are  applied 
to  pilot  control  of  hovering  vehicles  and  driver  control  of  highway  vehicles .  The 
SES  with  its  broad  beam  and  side  wall  keels  exhibits  little  rolling  and  side¬ 
slipping  in  common  with  catamarans  and  highway  vehicles. 

The  helm  itself  was  alt  in.  diameter  marine  steering  wheel.  An  artificial 
feel  spring  was  provided,  and  the  wheel  breakout  torque  and  the  kinetic  Coulomb 
friction  damping  torque  were  adjusted  to  be  less  than  those  values  currently  used 
in  order  to  examine  more  critically  the  possibilities  for  neuromuscular  coupling, 
biodynamic  amplification  thereof,  and  closed-loop  steering  stability  limitations . 
The  steering  wheel  rotation  from  stop  to  stop  was  270  deg  and  the  corresponding 
rudder  travel,  30  deg,  thereby  giving  a  9: 1  steering  ratio  as  in  contemporary 
practice . 

The  helmsman  was  asked  to  perform  the  steering-only  task  with  both  hands  on 
the  wheel  to  minimice  the  possibility  of  biodynamic  coupling.  In  addition  to  the 
primary  steering  control  task,  a  sub-critical*  speed  regulation  task  provided  a 
secondary  surrogate  for  trimming  the  water  speed  of  the  craft.  Speed  was  regulated 
with  the  helmsman's  left  hand  on  a  friction- restrained  quadrant  throttle  lever  with 
60  deg  of  fore-aft  travel,  so  that  when  the  helmsman  was  instructed  to  perform  both 
steering  and  speed  control  tasks,  he  maintained  only  his  right  hand  on  the  steering 
wheel. 

The  ship's  lateral-directional  dynamic  motions  in  response  to  the  helm  and 
disturbance  were  represented  by  linear  constant  coefficient  differential  perturba¬ 
tion  equations  in  body  axes  with  respect  to  a  trimmed  condition  at  constant  speed. 
The  ship's  cab  was  assumed  to  be  located  at  the  center  of  gravity  in  developing  the 
lateral-directional  transfer  functions  given  in  Table  1  for  a  trimmed  cruising 
condition. 

The  differential  equation  representing  the  divergent  (unstable)  controlled 
element  for  speed  regulation,  which  was  independent  of  the  equations  of  lateral- 
directional  craft  motion,  is  given  by  Eq.  1: 


U  =  Xu  +  Kg6t  ( 1 ) 

The  inverse  time  constant  X  was  chosen  as  0.1  rad/sec  for  the  simulation  to  approxi¬ 
mate  a  slow  sub-critical  divergence  which  would  require  consistent  but  not  over¬ 
whelming  attention  to  the  side  task  of  speed  regulation.  Comments  by  the  helmsman 
attested  that  the  speed  regulation  task  and  the  steering  task  together  saturated 
the  helmsman's  workload.  Yet  there  was  only  one  out  of  four  runs  with  both  tasks 
where  performance  on  the  steering  task  appeared  to  degrade. 

The  helmsman's  tasks  are  summarized  in  Fig.  2.  Approximations  to  the  lateral 
dynamics  of  Table  1  are  shown  to  illustrate  the  low-frequency  (i.e.,  long  time 
constant)  nature  of  the  SES  dynamics.  Also,  there  is  a  non-minimum  phase  term  in 
the  lateral  displacement  dynamics  such  that  the  ship  initially  translates  in  a 
direction  opposite  to  its  final  direction  for  a  given  rudder  command  as  illustrated 
in  Fig.  3  by  the  time  history  of  path  response  to  a  step  change  in  rudder  angle. 


* "Sub-critical"  is  used  in  the  manual  control  context  here  and  must  not  be 
confused  with  the  critical  speed  for  the  SES.  A  sub-critical  task  in  the  manual 
control  context  means  that  the  rate  of  divergence  of  the  open- loop  controlled 
element  as  characterized  by  X  is  below  the  critical  level  which  is  the  limit 
of  human  manual  control  capability. 
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tasks.  Participating  helmsmen  were  among  SES  crewmen  with  concurrent  operational 
experience  from  the  Surface  Effect  Ship  Test  Facility,  Naval  Air  Test  Center, 
Patuxent  River,  Maryland. 

Our  purpose  here  will  be  to  report  and  interpret  some  measurements  of  helms¬ 
men's  describing  functions  during  a  precision  straight-course-keeping  task  at  high 
cruising  speeds  in  a  disturbed  sea  and  involving  the  use  of  a  compensatory  elec¬ 
tronic  horizontal  situation  display  of  heading  and  lateral  displacement  errors  not 
specific  to  any  surface  effecs  ship. 

APPARATUS  AND  METHOD 

A  complement  of  3  in.  diameter  rotary  dial  instruments  provided  the  helmsman 
with  compass  heading,  turn  rate,  rudder  angle,  and  water  speed.  A  visual  field 
simulator  provided  a  collimated  external  view  of  the  moving  seascape  and  horizon 
properly  correlated  with  the  ship's  motions  experienced  in  the  cab  on  the  moving 
base.  In  addition,  an  optional  automatic  helm  was  provided  for  relief  of  manual 
steering  duties  while  subjective  rating  analyses  were  being  written  and  while  the 
helmsman  kept  watch  during  four-hour  missions. 

The  random- like  steering  disturbance  was  a  sum  of  five  non-harmonically  re¬ 
lated  and  randomly  phased  sinusoids  whose  relative  amplitudes  were  approximately 
inversely  proportional  to  frequency.  This  disturbance  was  generated  within  the 
NASA-STI  Mark  II  Describing  Function  Analyzer  ( DFA ) ( ^ ( 5 )  and  applied  to  the  mathe¬ 
matical  model  of  the  ship's  rudder  in  linear  combination  with  the  helmsman's  rudder 
command  signal  as  shown  in  Fig.  1.  The  disturbance  did  not  move  the  helmsman's 
wheel  nor  was  it  directly  visible  on  the  helmsman's  rudder  angle  indicator,  but 
its  effects  on  turn  rate,  heading  deviation,  and  lateral  displacement  from  the 
desired  course  were  observable.  The  helmsman  was  instructed  to  minimize  his  lateral 
error  during  each  100  sec  interval  when  the  DFA  was  used  for  cybernetic  performance 
measurement . 
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Figure  1 .  Block  Diagram  of  Heading  and  Lateral  Displacement  Control 
Task  as  Configured  for  the  Describing  Function  Analyzer. 

The  DFA  computes  on-line  the  finite  Fourier  transform,  mean-square,  and  mean 
of  a  signal  in  the  control  loop  —  in  this  case,  the  rudder  error  Se  in  Fig.  1  at 
each  of  the  five  input  frequencies.  The  final  describing  function  and  error  vari¬ 
ance  are  computed  off-line  with  a  digital  computer  program  based  on  some  of  the 
techniques  in  Aliens 5).  The  principles  of  the  measurement  of  helmsman/SES 
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MANUAL  STEERING  OF  A  SIMULATED  SURFACE-EFFECT  SHIP' 


by  Warren  F.  Clement 
and  R.  Wade  Allen 
Systems  Technology,  Inc. 


ABSTRACT 

Crew  performance  in  the  motion  environment  of  a  large  generic  high  speed  sur¬ 
face  effect  ship  has  been  investigated  by  means  of  a  motion  base  simulation.  Some 
of  the  helmsman's  control  tasks  have  been  addressed  with  the  aid  of  a  simulated 
external  forward  visual  field  of  the  seascape  and  navigation  and  steering  displays 
in  the  pilot  house.  In  addition  to  the  primary  steering  control  tasks,  a  sub- 
critical  speed  tracking  task  provided  a  secondary  surrogate  for  trimming  the  water 
speed  of  the  ship.  The  results  of  helmsmen's  steering  describing  function  measure¬ 
ments  are  presented,  and  some  suggestions  for  their  interpretation  are  offered. 

The  likely  steering  loop  closures  comprise  heading  and  lateral  displacement  for 
the  course-keeping  task  investigated.  Regardless  of  the  influence  of  workload, 
steering  technique,  water  speed,  and  sea  state,  the  helmsmen  apparently  adopted  a 
disturbance  regulation  bandwidth  of  about  0.2  rad/sec  for  lateral  displacement. 
Suggestions  for  reducing  the  variability  in  future  helmsmen's  measurements  are 
offered. 

INTRODUCTION 

The  Surface  Effect  Ship  (SES)  is  an  ocean-going  vessel  employing  a  self¬ 
generated  aerostatic  cushion  in  contact  with  the  water  surface  for  vertical  sup¬ 
port  ).  The  SES  has  rigid  shallow-draft  side  walls  with  flexible  fore  and  aft 
skirts  or  seals  to  contain  the  pressurized  air  cushion  while  permitting  the  passage 
of  surface  waves  through  the  cushion  plenum.  The  side  walls  serve  as  keels  to 
provide  lateral  stability  in  the  manner  of  a  catamaran.  Because  the  SES  rides  on 
a  cushion  of  air,  it  is  less  subject  to  the  drag  penalties  which  limit  the  speed 
of  displacement-hull  vessels.  Consequently,  the  SES  is  capable  of  sustained  higher 
speeds,  and  requires  precision  in  course-  and  sea-keeping,  especially  in  aggravated 
sea  conditions. 

The  ability  of  crewmen  to  perform  shipboard  duties  without  undue  fatigue  or 
decreased  proficiency  has  been  the  subject  of  recent  investigations  with  manned 
motion  base  simulation  using  motion  predictions  from  a  2000-ton  SES  mathematical 
raodel(2)(3).  A  simulated  mission  profile  with  assigned  crew  tasks  provided  a 
disciplined  scenario  for  measuring  crew  performance  in  the  simulated  ship  motion 
environment.  Various  tasks  involving  facsimile  shipboard  operations  at  four  duty 
stations  were  performed.  One  of  the  duty  stations  in  the  pilot  house  is  that  of 
helmsman.  In  the  simulation  the  helmsman's  assignment  included  a  division  of  at¬ 
tention  among  steering,  speed  regulation,  obstacle  avoidance,  and  communication 
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tronic  complexity.  Still,  when  advanced  ship  designs  are  considered  or  if  critical 
controls  are  to  be  applied,  a  number  of  potential  advantages  exists.  A  few  of  these 
advantages  have  been  cited  in  the  preceeding  material  and  and  their  implementation 
indicated.  A  particularly  attractive  feature  of  the  discribed  implementation  is  the 
ability  to  adjust  the  algorithms  by  means  of  weighting  functions.  Potential  approaches 
for  automating  the  key  ship  control  functions  of  maneuvering  and  ride  control  were 
also  identified.  It  is  hoped  that  the  future  will  provide  the  opportunity  to  try  some 
of  these  approaches  in  an  ocean  going  SES. 
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Figure  16.  Standard  deviation  of  heave  acceleration 
as  a  function  of  controller  gains. 
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When  a  compensation  network  (Figure  13)  with  gains  K^,  1^,  K.^  and  is  incly^d 

in  this  same  formulation,  the  resulting  equations  for  these  key  parameters  become: 


AC  +  C3- 


(Achc  + 


cl)_1 

UL  -  (Ac(|4)+  CJ  ^Achc  +  ci)'1 

a  Z 


(17) 

(18) 
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-  +  _k  La  - 

-=~  2mv  2AP 

2wt,  H  C 
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where: 

8Q  =  Slope  of  the  VG  fans  (CFS/%) 
3VG 


(19) 


(20) 
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AcKi 

^  *  dVG) 


j=  1,  2,3,4 


Based  on  these  relationships  the  individual  gains  are  seen  to  have  the  following 


ef  fects: 

1. 

As 

Ki 

K2 

increases 

2. 

As 

increases 

3. 

As 

K4 

increases 

4. 

As 

increases 

Such  a  combination  of  gains  should  allow  "fine  tuning"  of  the  RCS  for  a  given 
sea  state  or  to  assure  uniform  response  over  a  band  of  sea  states.  An  example  of 
this  tuning  is  shown  in  Figure  16  which  indicates  the  standard  deviation  of  heave 
acceleration  as  a  function  of  and  with  K.  =  ^  =  0  for  a  ship  with  the  con¬ 

figuration  indicated  in  Table  1.  As  is  evident,  the  RCS  has  the  capability  to 
reduce  heave  acceleration  by  a  factor  of  two. 


Table  1.  Parameters  for  Ship  and  RCS 


PARAMETER 


VALUE 


A 


C 


m 


18,790  ft  2 

166,855  slugs 
254.3  PSF 
35,000  CFS 

632.6  CFS/2! 
5740  CFS/ ft 
349,923  ft3 


CONCLUSION 


A  number  of  tradeoffs  with 
required  before  committing  to  a 
some  of  these  tradeoffs  may  not 


respect  to  cost,  reliability  and  maintainability  are 
given  level  of  automation.  Based  on  past  experience 
reflect  too  favorably  on  the  use  of  increasing  elec- 
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SYMBOL 


a  Naperian  Base,  2.71828... 

E  Error  signal  output  to  DFA  (I  -  M) 

E/l  DFA  error- to- input  describing  function 

F  Helmsman's  heading  feedforward  describing  function  including  neuro¬ 

muscular  steering  actuation  dynamics  l’ dimensionless) 

I  DFA  disturbance  input  to  rudder  steering  axis 

j  vA' 

K0  Throttle  control- to-speed  response  gain  (ft/sec-red) 

K  Gain  equalization  in  (dimensionless) 

M  Closed- loop  motion  output  with  respect  to  DFA  input 
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M/E 

NSr 

*r 


Open- loop  output  to  error  describing  function  as  measured  by  DFA 

Controlled  element  transfer  function  numerator  polynomial  representing 
yaw  rate  response  to  rudder  displacement  (l/sec) 

Controlled  element  transfer  function  numerator  polynomical  representing 
lateral  acceleration  response  to  rudder  displacement  (ft/sec^-rad) 

Laplace  operator,  o  +  ju> 

Time  constant  of  lead  equalization  in  or  F^  (sec) 

Perturbed  longitudinal  water  speed  of  craft  with  respect  to  trimmed 
speed  (ft/see) 

Lateral  displacement  of  craft  (ft) 

Lateral  displacement  command  (ft) 

Lateral  displacement  error  (y  -  y) 

Helmsman's  lateral  displacement  describing  function  (rad/ft) 

Helmsman's  heading  feedback  describing  function  (dimensionless) 

Rudder  angle  error  (rad) 

Perturbed  rudder  deflection  angle  with  respect  to  trimmed  angle  (rad) 
Throttle  displacement  (rad) 

Characteristic  determinant,  transfer  function  denominator 


?jj  Damping  ratio  of  second-order  lag  in  F^  representing  effective  neuro¬ 

muscular  actuation  dynamics  modified  by  proprioceptive  or  visual 
feedback  (dimensionless) 

X  Inverse  time  constant  of  the  first-order  sub-critical  tracking  task 

(rad/sec) 

o  Real  part  of  the  complex  variable  s  (rad/sec) 

Effective  helmsman's  time  delay,  including  transport,  equalization, 
and  scanning  contributions  (sec) 

i(i  Heading  of  craft  (rad) 

t>c  Heading  command  (YyYe)  (rad) 

Ve  Heading  error  (\|/c  -  <r)  (rad) 

o>  Circular  frequency;  imaginary  part  of  the  complex  variable  s  (rad/sec) 

o>c  Unit  gain  crossover  frequency  (rad/sec) 

Undamped  natural  frequency  of  second-order  lag  in  F.  representing 
effective  neuromuscular  actuation  dynamics  modified^by  proprioceptive 
or  visual  feedback  (rad/sec) 

cnu  Unstable  phase  crossover  frequency  (rad/sec) 
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abbreviations 


dB 

deg 

DFA 

ft 

in. 

rad 

sec 

SES 

(•) 


Decibel 

Degree 

Describing  Function  Analyzer 

Foot 

Inch 

Radian 

Second. 

Surface  effect  ship 

(raised  period)  Time  Derivative  Operator  d/dt  (l/sec) 
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SYSTEM  ANALYSIS  TECHNIQUES  FOR  DESIGNING  RIDE  CONTROL 
SYSTEMS  FOR  SES  CRAFT  IN  WAVES 


by  Paul  Kaplan,  Oceanics,  Inc. 
and  Sydney  Davis,  NAVSEA,  PMS-304 


ABSTRACT 

Procedures  for  design  of  a  control  system  (known  as  Ride  Control 
System  or  RCS)  for  reducing  the  vertical  plane  accelerations  of  SES 
craft  by  means  of  altering  the  pressure  variation  within  the  cushion 
are  described.  An  outline  is  given  of  a  method  using  control  system 
analysis  techniques  in  terms  of  the  craft  fundamental  dynamics  in  its 
uncontrolled  mode  (fundamental  frequency  and  damping  characteristics) . 
Particular  feedback  signals  in  terms  of  craft  state  variables  are 
established  in  this  simplified  analysis,  which  is  complicated  by 
saturation  limits  of  the  control  effector,  lags  in  actuator  displace¬ 
ment,  etc.  (an  illustration  of  a  nonlinear  saturation  limit  is  the 
size  of  the  deck  leakage  area  louvers  or  vent  valves  used  as  a  typical 
means  of  SES  control) .  The  analytical  model  is  described  and  initial¬ 
ly  applied  to  louver  systems,  with  and  without  bias  in  the  area  open¬ 
ings.  The  system  analysis  technique  provides  a  measure  of  rms  accel¬ 
erations,  demanded  control  output  and  a  means  of  obtaining  optimal 
responses  as  a  function  of  the  control  gains  in  a  parametric  study. 
Illustrations  are  given  for  control  via  controlled  area  louver  systems; 
axial  fans  with  variable  blade  angles;  and  variable  geometry  centrif¬ 
ugal  fans.  The  basic  system  analysis  technique  illustrated  in  the 
present  paper  have  a  wide  range  of  applicability  to  the  different 
possible  concepts  used  for  establishing  RCS  designs  for  SES  craft, 
and  provide  the  initial  gains  for  establishing  an  optimal  control 
system. 

INTRODUCTION 

The  surface  effect  ship  (SES)  is  a  seagoing  vehicle  that  operates 
at  high  speed  at  the  interface  region  between  air  and  water,  with  its 
major  means  of  weight  support  obtained  from  the  contained  air  pressure 
within  a  cushion  region  under  the  main  deck.  The  concept  underlying 
the  SES  craft,  in  contrast  to  that  of  other  vehicles  such  as  hover¬ 
craft,  involves  the  use  of  rigid  sidewalls  intersecting  the  water 
surface  and  the  presence  of  fore  and  aft  seals  (bow  and  stern)  which 
act  to  retain  the  pressurized  air  within  the  so-called  cushion  region. 
The  resulting  motions  of  this  craft  when  operating  in  waves  are  thus 
primarily  due  to  dynamic  effects  that  arise  from  the  changing  pressures 
associated  with  the  action  of  the  waves  in  modifying  the  cushion  volume. 
A  simplified  description  of  the  basic  dynamics  that  influence  the  .  . 
vertical  plane  motions  of  SES  craft  in  waves  has  been  presented  in'  , 
with  primary  emphasis  on  the  resulting  heave  accelerations  of  the  craft. 

The  high  level  of  vertical  accelerations  experienced  by  SES  craft, 
as  indicated  from  both  experimental  and  computer  simulation  studies, 
has  illustrated  the  necessity  of  incorporating  some  form  of  control 
system  for  reducing  these  acceleration  levels.  The  main  purpose  in 
reducing  these  heave  acceleration  levels  is  to  satisfy  the  requirements 
for  human  operator  habitability  on  board  such  craft. 
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While  the  precise  level  for 

SES  craft  acceleration  spectra  is  being  investigated ' 2 )  for  different 
missions,  the  goals  for  the  action  of  a  control  system  are  primarily 
aimed  at  reducing  the  heave  acceleration  levels  to  the  smallest  pos¬ 
sible  value  while  still  considering  the  possible  penalties  associated 
with  the  use  of  such  a  control  system. 

In  order  to  obtain  an  optimal  control  performance  a  fundamental 
investigation  has  to  be  carried  out  in  order  to  consider  the  effects 
of  different  physical  means  of  control,  as  well  as  to  establish  the 
particular  control  rules  that  should  be  used  in  order  to  obtain  the 
best  performance.  In  addition  to  such  considerations,  recognization 
of  the  realistic  limits  of  actuation  of  different  devices  on  actual 
craft  must  be  made  together  with  consideration  of  lags  and  other 
limiting  effects.  All  of  these  aspects  should  be  investigated  in 
order  to  determine  their  influence  on  the  ultimate  craft  performance. 

In  addition,  in  any  overall  investigation  of  control  systems  for 
reducing  SES  craft  heave  acceleration,  a  detailed  examination  of  the 
power  penalty  associated  with  the  particular  type  of  control  must  be 
made  in  order  to  obtain  a  measure  of  the  overall  efficiencies  of  such 
control  action,  as  the  power  requirements  are  an  obvious  practical 
limitation. 

The  present  paper  will  only  consider  an  initial  analysis  technique 
for  establishing  a  control  system  for  SES  craft,  which  serves  as  a 
basic  approach  to  '.efine  the  nature  of  the  control  system  signal 
processing  associated  with  different  physical  means  of  control,  i.e. 
different  control  effectors  and/or  mechanisms.  The  basic  system 
analysis  technique  relies  upon  utilizing  knowledge  of  the  fundamental 
dynamics  of  the  craft  vertical  plane  motion  presented  in'1'  in  order 
to  provide  the  initial  gains  and  control  rule  for  establishing  an 
optimal  control  system.  The  information  obtained  from  the  method 
described  here  can  then  be  used  for  more  detailed  studies  of  system 
performance  with  a  complete  nonlinear  six  degree  of  freedom  motions 
program  foi  SES,  as  described  in'3',  which  can  then  be  used  to  explore 
a  more  precise  measure  of  motion  performance  for  particular  craft  with 
different  control  concepts.  In  view  of  the  more  simple  approach  based 
on  system  analysis  that  is  presented  in  this  paper,  no  consideration 
of  the  power  penalty  associated  with  different  control  concepts  is 
presented  since  it  is  beyond  the  capabilities  of  that  method  and  is 
more  directly  associated  with  a  detailed  study  using  the  large  computer 
program. 

ASSUMPTIONS  AND  TECHNIQUES  OF  ANALYSIS 

In  carrying  out  the  present  analysis  a  number  of  simplifying 
assumptions  are  made,  most  of  which  have  been  applied  in  the  basic 
dynamics  study  of  (O  which  is  used  as  a  starting  point  in  the  present 
work.  Since  the  major  concern  here  is  the  craft  heave  acceleration, 
the  craft  motion  is  only  represented  by  the  heave  degree  of  freedom 
which  is  coupled  only  with  the  pressure  variable.  Linear  craft  motion 
equations  are  used,  with  the  only  nonlinearity  being  due  to  the  control 
limits  or  other  aspects  of  the  control  action.  There  is  no  additional 
leakage  due  to  motions,  with  the  seals  assumed  to  be  perfect  sealing 
elements  throughout  any  motion  conditions.  The  only  excitation  or 
disturbance  action  that  the  control  system  affects  is  that  due  to  the 
wave  pumping  action,  which  is  evaluated  for  the  head  sea  case  that  has 
the  greatest  expected  heave  acceleration  response.  There  are  no  hydro¬ 
static  and/or  hydrodynamic  forces  acting  on  che  craft,  with  only  the 
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air  pressure  providing  the  support  force  and  disturbance  reactions. 

Thus  the  means  of  control  action  must  therefore  involve  some  method 
of  altering  the  air  pressure  by  affecting  the  air  flow  characteristics 
in  the  cushion. 

The  method  of  analysis  used  is  via  frequency  response  methods, 
together  with  a  quasi-linear  method  of  representing  nonlinear  effects 
which  will  be  commensurate  with  the  methods  used  for  analysis  of  the 
(primary)  linear  elements  in  the  system.  The  SES  craft  will  be 
assumed  to  operate  at  some  forward  speed  in  different  random  sea  states 
represented  by  the  Pierson-Moskowitz  wave  spectrum,  so  that  stochastic 
process  outputs  in  statistical  form  will  be  obtained.  While  the  major 
disturbance  acting  on  the  craft  is  due  to  wave  pumping  action  in  the 
cushion  region,  some  considerations  of  transient  response  properties 
of  the  controlled  system  will  also  be  made.  The  transient  response 
properties  will  be  shown  to  also  influence  the  response  to  the  wave 
distrubances ,  and  in  addition  some  discussion  will  also  be  given  to 
the  responses  of  the  controlled  craft  due  to  leakage  arising  from 
craft  motions  (primarily  pitch  motions) . 


EQUATIONS  OF  MOTION 

The  craft  equations  of  motion  are  given  in  in  the  form 

mz  +  A,  p  p  =  0  (1) 

bo 

and 


K1w  +  K3y  -  aAfez 


-  K2  AAL  "  paVb 


(2) 


where 

m 


K1'K2'K3 


b 

waves 


craft  mass 
cushion  area 

equilibrium  cushion  pressure  (gage) 

p  ~ pQ  ,  nondimensional  pressure  variation 
po 

heave  motion 
density  of  air 

change  in  leakage  area  (controlled) 

system  parameters  from  fan  characteristics,  cushion 
geometry,  equilibrium  pressure 

time  rate  of  change  of  volume  pumping  due  to  waves 


This  equation  system  is  derived  under  the  conditions  of  an 
(effective)  linear  fan  characteristic  curve.  The  fan  curve  is  repre 
sented  by 


Qin 


“  Qo  + 


(H)  po» 


(3) 


where  Q0  is  the  fan  flow  rate  when  p  =  p_  and  ( |^ ) 
linear  fan  slope  (inverse) ,  as  indicated  in  p  ° 


is  the  assumed 
Figure  1. 


A 
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Figure  1.  Linear  Fan  Characteristic  Curve. 


The  parameters  K^,  K 2  and  are  defined  by 


where 

hfa  =  height  of  cushion 

ir  =  ratio  of  specific  heats  for  air 

Pa  =  atmospheric  pressure 

c  =  orifice  coefficient 
n 

The  quantity  Vb  is  defined  for  a  regular  sinusoidal  head 

sea  wave  as  waves 


sin  t\ 9/ \ 


cos  w  t 
e 


where 

a  =  wave  amplitude 
1  =  cushion  length 

i  =  wavelength 

and  we  is  the  (circular)  frequency  of  encounter  given  by 

.  w2u  ,c, 

e  g 

with  oj  =  wave  frequency  and  u  =  craft  forward  speed. 

This  particular  system  of  equations  is  appropriate  to  the  case 
where  additional  controlled  leakage  is  provided  via  louver  area 
variation.  However  other  possible  control  system  implementation 
concepts  can  also  be  represented  using  these  equations,  as  will  be 
shown  later  in  this  paper. 
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CONTROL  SYSTEM  ANALYSIS 


The  design  of  a  useful  system  for  reducing  heave  accelerations 
in  waves,  i.e.  a  ride  control  system  (RCS)  ,  must  consider  a  number  of 
different  items  of  importance.  Low  values  of  rms  heave  acceleration 
must  be  achieved  together  with  good  dynamic  response  to  impulsive 
transient  disturbances,  while  not  demanding  too  large  a  fan  power 
increase.  The  use  of  different  optional  control  techniques,  either 
in  the  frequency  or  state  variable  domain,  would  lead  to  complicated 
compensation  networks  that  are  obtained  from  lengthy  computations.  In 
addition  the  particular  optional  control  would  vary  significantly  for 
different  craft  operating  conditions,  i.e.  speed-sea  state  conditions. 
Also  whenever  large  sea  conditions  are  encountered  the  entire  utility 
of  a  linear  analysis  is  subject  to  questions. 

A  simple  representation  of  the  equations  of  motion  in  state  vari¬ 
able  form  is  given  by 

5  =  Az  +  b  AL  +  cVb  (7) 


■W-  s 

■  [;  3 


and  the  following  definitions  apply: 

PaAb  _  _  K3 


b  =  -SI 


with  this  form,  and  the  nature  of  the  coefficient  matrices,  the 
representation  of  a  simply-structured  control  scheme  can  be  achieved 
in  terms  of  state  variable  feedback  to  achieve  a  desired  form  of 
characteristic  polynomial  equation  describing  the  system.  Since  the 
state  vector  can  be  easily  measured,  the  control  scheme  considered  is 

AAl  =  k^z  +  k2u  (10) 

c 


-  ki  -  k2 

k  —  —  k  =  — 

K1  b  '  2  b 

and  AA  is  the  commanded  change  in  leakage  area. 
Lc 

Control  Analysis  for  Louver  Control 


The  technique  for  control  via  louvers  involves  changing  the 
leakage  area  through  variation  of  an  additional  opening  (in  the  deck 
usually).  This  procedure  allows  controlled  air  leakage  when  the 


F2  4-5 


cushion  pressure  rises,  but  there  is  no  way  to  increase  pressure  by 
means  of  an  additional  leakage  opening.  However  the  most  effective 
means  to  achieve  a  useful  heave  acceleration  reduction  via  louvers 
is  in  allowing  the  area  change  to  occur  about  some  reference  condition, 
so  that  area  variation  in  both  directions,  i.e.  positive  and  negative, 
c;n  be  achieved  (referred  to  as  2-sided  control)  .  The  mean  power 
requirements  increase  significantly  for  this  type  of  control,  depend¬ 
ing  on  the  mean  effective  position  of  the  louver  area  opening,  i.e. 
if  that  value  is  the  same  as  the  base  equilibrium  leakage  area  of  the 
craft  then  the  fan  power  is  doubled.  Thus  it  is  necessary  to  limit 
the  extent  of  the  equilibrium  louver  opening  (and  hence  the  total  area 
variation)  in  order  to  achieve  a  reduced  additional  power  due  to  the 
control,  while  still  obtaining  useful  heave  acceleration  reduction. 

This  is  a  typical  trade-off  illustration,  whose  exact  consideration 
cannot  be  treated  in  this  simplified  analysis  but  only  from  some 
detailed  simulation  studies  (i.e.  with  the  aid  of  the  large  digital 
computer  program) . 

When  applying  the  control  rule  of  Eq.  (10)  to  this  case,  consid¬ 
eration  must  be  given  to  the  limits  imposed  by  the  maximum  leakage 
area  permitted  in  both  directions  (relative  to  the  equilibrium 
reference  louver  area  opening) .  This  is  illustrated  by  the  symmetrical 
saturation  limit  for  the  actual  leakage  area,  as  shown  in  Figure  2. 


controlled  leakage  area  saturation 


Figure  2.  Model  of  Simplified  Heave  Motion  Dynamics  Including 
Controlled  Leakage  Area  Saturation. 
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The  analysis  where  the  wave  inputs  arise,  from  a  random  sea  in¬ 
volves  a  zero-mean,  Gaussian  stochastic  process  where  the  action  of 
the  nonlinear  limit  is  represented  by  the  use  of  a  describing  function 
for  random  inputs.  The  describing  function  for  this  symmetric 
saturation  nonlinearity  is  a  representation  as  a  quasi-linear  gain, 
given  by 


where 


K  (o ) 


L/o/2 


(A  Al  ) 
\  conn./ 


(12) 


(13) 


is  the  mean  square  value  of  the  commanded  leakage  area  change. 


Using  standard  Laplace  transform  notation  for  the  transfer 
functions,  the  following  relations  are  obtained: 


/AA  \  =  [T  ( jw  )  |  T2(w)  S  (u)du 

V  comm . /  *  o  n  n 

J  I2  (w)  S  (u)doj 
^  e  n  n 


-2  ^ 

<V  =  z  =  g 


(14) 

(15) 


where  s  =  ju  is  used  in  the  transfer  functions  T.  ( jui  )  and  T.  (ju  ); 
the  quantity  Tn(w)  represents  the  wave  pumping  vOluml  transfer  e 
function  in  terms  of  wave  frequency  w;  Sn  (w)  is  the  surface  wave 
elevation  spectrum;  and  the  frequencies  are  related  by  Eq .  (6).  The 
transfer  functions  Tj  and  T2  are  defined  by 


T^  ( S) 
T2  (s) 


p  (k2s-k1g> _ 

s2  -  [a2  +  k2K(a)  Is  +  [ajg  +  k  K(u)I 

_  cs2 _ 

s2  -  [a2  +  k2K  (o)]s  +  [a1g  +  k1K(a)I 


(16) 

(17) 


and  the  function  T _^  ( o )  is  given  by 

2gA, 


T  = 


o2l 


SI 


(18) 


It  can  be  seen  from  the  denominators  of  the  transfer  functions  T  and 
T2  where  the  effect  of  the  control  gains  (k,  and  k2)  in  altering1the 
basic  natural  frequency  and  damping  of  the  system  response  is  indicated, 
together  with  the  effect  of  the  nonlinear  limit  operation  in  limiting 
the  influence  of  these  parameters.  The  solution  of  the  above  equations 
requires  iteration  since  the  quantity  a  entering  into  the  definition 
of  K (o)  is  not  known. 


The  solutions  are  obtained  for  different  operating  conditions  of 
sea  state  and  ship  speed.  An  effective  manner  of  presenting  the 
results  is  in  terms  of  dynamic  response  parameters  such  as  natural 
frequency  and  damping,  which  are  obtained  in  the  linear  range  (when 
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K(  )  =1)  with  the  definitions 


0>n  =  /Ta~tlc^Tg  (19) 

C  =  -  (a2  +  k2)/2un  (20) 

with  c  the  damping  ratio  of  the  system. 

Computations  can  be  carried  out  for  particular  cases  where  the 
reference  louver  area  opening  is  some  particular  fraction  of  the  craft 
equilibrium  base  leakage  area  Alq-  T^e  condition  for  symmetric  action 
of  the  louver  area  variation,  as  indicated  in  Figure  2  and  the  present 
analysis,  corresponds  to  an  opening  that  is  1/2  of  the  maximum  possible 
louver  area.  Thus  if  the  saturation  limits  are  ±L,  as  shown  in 
Figure  2,  the  equilibrium  reference  louver  area  is  L  and  the  area  can 
vary  from  0  to  2L.  The  selection  of  the  value  of  L  is  usually  refer¬ 
red  to  the  quantity  Aj,0  for  ease  of  analysis. 

The  computations  are  carried  out  for  different  values  of  k-^,  which 
would  affect  the  linear  natural  frequency,  with  variations  of  k^  to 
produce  desired  values  of  linear  damping  ratio  for  transient  response 
considerations.  All  of  these  computations  are  carried  out  for  each 
particular  speed-sea  state  operating  conditions  in  head  seas.  Results 
of  computations  for  the  louver  control,  as  well  as  other  means  of 
control  to  be  discussed  in  the  following  paragraphs,  are  given  in  a 
later  section  of  this  paper. 

Analysis  of  Fan  Blade  Angle  Control 

Another  possible  method  for  reducing  heave  acceleration  is  by 
direct  action  of  the  fan,  whereby  the  changes  in  certain  characteristics 
of  the  fan  would  result  in  changes  in  the  pressure  acting  in  the  craft 
cushion.  The  particular  fan  characteristic  that  can  be  altered  in  the 
case  of  axial  fans  is  the  pitch  angle  of  the  blades  during  fan  rota¬ 
tion,  which  then  provides  changes  in  the  pressure  and  flow  rate 
properties  of  the  fan  in  accordance  with  the  particular  blade  angle 
value.  Representative  fan  maps  relating  pressure  as  a  function  of  the 
flow  rate  for  different  blade  angles  are  shown  in  Figure  3,  as  obtained 
from  curves  supplied  by  a  manufacturer  of  axial  fans,  which  are  then 
used  to  investigate  this  control  concept. 

For  purposes  of  the  present  analysis,  the  idealization  of  these 
curves  is  an  interpretation  whereby  the  curves  have  the  same  basic 
slope  (in  the  p-Q  fan  maps)  for  the  different  values  of  blade  angle 
when  using  an  assumed  linear  fan  map.  In  the  basic  derivation  of 
Eq.  (2),  as  originally  presented  in  O),  the  flow  characteristic  of 
the  fan  inflow  to  the  cushion  is  represented  by 

°in  =  Qo  +  (%\  W  (21) 


where  this  relationship  is  assumed  to  hold  for  a  particular  fan  blade 
angle,  a  =  aref.,  which  is  the  reference  blade  angle  about  which  the 
changes  in  angle  take  place.  The  control  of  blade  angle  is  assumed 
to  affect  only  the  quantity  Cp  in  Eq.  (21),  so  that  the  representation 
of  the  fan  flow  into  the  cushion  is 


Qin 


('It  40  ) 

1  + 

(|2)  P  u 

(22 

'  “ref J 

\9p  4  ° 
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with  Act  the  change  in  the  fan  blade  angle.  The  value  of  4a  is  limited 
such  that 

1  +  —  >0  (23) 

a  c  — 
ref . 

and  also  with  an  upper  limit  of  —  _<  1,  leading  to 

“ref. 

-1  <  —  <1  (24) 

“ref- 

which  is  the  complete  limit  representation. 

Examining  the  derivation  in('),  the  representation  of  the  cushion 
outflow  relationship  there  is 

p  Q  ^  =  p  Q  +  4  K_A  u  +  K-t  AA7  ,  (25) 

a  out  a  o  2  2  Lq  2  L 


PaQo  =  K2Al  (26) 

o 

With  the  fan  flow  into  the  plenum  given  by  Eq.  (22)  this  leads  to 


P  (Q-  -Q  .)  =  K0At  — 

Ma  out  2  L  a _ * 

o  ref . 


k3p  -  k2aal 


For  the  control  case  where  there  is  no  commanded  leakage  (i.e.  louver 
control) ,  the  quantity  AA^  -  0,  but  it  can  be  seen  from  Eq.  (27)  that 
the  quantity  Al  Aot/ctrej,  “is  essentially  equivalent  to  the  controlled 

leakage  area  chSnge  if  the  sign  of  the  control  gains  is  reversed. 

Thus  any  results  obtained  by  the  use  of  a  two-sided  louver  control, 
with  the  AA  limited  by  +AT  ,  will  be  exactly  the  same  as  from  the 

present  fan  blade  angle  except  for  the  sign  change  of  the  control  gains. 

The  equations  of  motion  for  this  case  with  fan  blade  angle  control 
result  in  Eq.  (2)  being  replaced  by 

K.p  +  K3u  -  PaAbz  =  K2Al  ~  -  PaVb  (28) 

o  “ref.  waves 

with  Eq.  (24)  also  applied,  with  the  control  rule  (command  signal) 
given  by 


(£«.) 

'  '  com . 


-k3z  -k2p 


These  equations  are  then  combined  in  the  same  manner  as  for  the  louver 
control,  using  a  describing  function  representation  for  the  symmetric 
saturation  limit  action  of  Eq.  (24) ,  and  applied  to  the  same  case 
treated  there,  i.e.  the  two-sided  louver  control  that  is  equivalent  to 
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the  present  problem  of  fan  blade  angle  control.  This  control  rule  is 
then  available  as  an  initial  optimal  control  form  to  be  applied  in  the 
system  analysis  of  fan  blade  angle  control,  using  values  of  A^  for  a 

particular  SES  craft  in  the  specified  operating  conditions  (speed-sea 
state  combinations) . 

Analysis  for  Fan  Area  Control 

The  use  of  variable  fan  areas  as  a  means  of  control  is  based  on  a 
recent  development  of  such  a  system  for  centrifugal  fans'I,,  .  The  fan 
inlet  area  is  changed  so  that  the  fan  output  flow  is  then  altered,  at 
almost  constant  pressure,  in  order  to  affect  the  resultant  pressure  in 
the  cushion.  A  general  representation  of  the  fan  map  for  fans  of  this 
type  is  shown  in  Figure  4.  The  three (3)  curves  in  Figure  4  represent 
the  basic  p-Q  variation  for  a  particular  degree  of  fan  inlet  area, 
denoted  as  1/2,  3/4  and  fully  open,  with  all  intermediate  positions 
also  possible.  If  the  condition  corresponding  to  3/4  open  is  con¬ 
sidered  as  a  reference  base  area,  the  possible  changes  in  area  to  1/2 
or  full  open  correspond  to  the  relation 


where  AA  is  the  change  in  fan  inlet  area. 

Since  the  curves  in  Figure  4  indicate  that  the  slope  of  the  fan 
maps  (p-Q  curves)  also  changes  with  the  area  change,  a  fit  of  these 
curves  was  found  to  be  represented  (approximately)  by 


Q  (Q  -AQ-u) (1+1.5  p  )  (31) 

ref. 


where  the  first  parenthesis  term  represents  the  flow  rate  characteristic 
for  the  fan  case  are  of  3/4  open  (the  reference  condition  in  the  present 
case) .  Other  possible  numerical  values  would  be  present  for  different 
fan  curves,  but  the  present  case  is  sufficient  for  illustrative  purposes. 


Since  the  action  of  the  fan  is  to  change  the  input  flow  to  the 
cushion,  and  the  leakage  (louver)  control  is  a  means  to  change  the  out¬ 
flow,  the  control  action  can  be  represented  similarly  with  only  a  change 
in  sign  for  the  command  gains.  The  change  in  the  basic  pressure 
equation  leads  to 


K1u+K3u-oaAbz= 


1.5K2Al 


^  +1 . 5o  U  V, 

Aref.  aV3p/o\Aref  J  al 


(32) 


where  A^  is  equilibrium  leakage  area,  with  the  control  command  signal 

,  o 

given  by 


klZ  -  k2v 


(33) 


subject  to  the  limit  given  by  Eq.  (30). 
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The  problem  here  is  to  find  the  values  of  the  gains  kj  and  k2  to 
provide  an  optional  control  for  this  case.  It  can  be  shown  that  the 
nonlinear  product  zp  that  results  from  Eq.  (32)  has  a  negligible 
contribution  (since  p  •v,  z  the  resulting  product  has  no  contribution 
to  a  describing  function  value)  and  can  be  neglected.  Similarly  the 
p2  product  resulting  from  the  nonlinearity  in  Eq.  (32)  can  be  shown  to 
have  an  influence  of  the  order  of  10%  on  the  equivalent  linear 
representation  of  the  control  input,  when  neglecting  the  effect  of  the 
saturation  limit  action.  This  treatment  is  based  upon  the  assumption 
of  an  equivalent  linearization  of  the  p2  term  by  use  of  the  mean  value 

of  p  that  occurs  for  SES  craft  (a  small  negative  value,  indicating 

deeper  draft  in  waves) ,  as  well  as  avoiding  consideration  of  the 
complexities  arising  from  multiple  nonlinearieites  in  a  describing 
function  analysis.  On  the  basis  of  the  above  reasoning  the  nonlinear¬ 
ity  due  to  the  product  expression  in  Eq .  (32)  is  ignored  (since  a  10% 
change  in  the  gain  factor  is  not  significant  within  the  level  of 
accur  -y  of  the  present  analysis! ,  while  the  nonlinear  limit  action 
indicated  by  Eq.  (30)  is  still  retained.  The  effect  of  such  a  treat¬ 
ment  can  only  be  accurately  assessed  by  means  of  a  more  complete 

simulation  study,  which  also  includes  other  effects  neglected  in  this 
overall  analysis,  so  that  the  results  obtained  here  are  only  to  be 
viewed  as  initial  estimates  of  control  gains  for  extended  simulation 
studies.  ,, 


Based  on  the  analysis  discussed  above,  the  actual  term  entering 
the  equation  system  is 


1. 


5AL  A~ 
o  ref 


1.5Al  (k1z  +  k2p) 
o 


(34) 


and  this  term  is  bounded  by  ±  4al  when  considering  the  action  of  the 
control  as  a  symmetric  variation  about  the  3/4  area  reference. 


When  considering  the  operation  of  the  variable  area  fans  as  two 
separate  groups,  i.e.  a  biased  control  with  half  the  total  fans 
operating  at  1/2  open  and  the  other  half  of  the  fan  system  operating 
at  full  open,  the  equations  are  changed  somewhat.  In  that  case  the 
1/2  open  fan  can  operate  (under  control)  only  to  increase  the  inlet 
area  and  operate  up  to  the  limit  of  the  full  open  condition,  i.e.  to 
allow  an  increase  in  flow  into  the  cushion,  while  the  full  open  fans 
can  operate  under  the  control  to  reduce  the  flow  into  the  cushion,  down 
to  the  1/2  open  condition.  This  operation  involves  knowing  an  equation 
for  the  fan  curve  for  each  operating  condition,  and  then  the  resulting 
equation  for  the  change  of  flow  rate  under  the  actions  of  the  control 
applied  to  each  operating  fan  group.  In  that  case  the  flow  in  can  be 
represented  by 


where  the  1/2  and  1  subscripts  co-respond  to  the  operating  fan  group 
(area  opening  value) ,  with  the  control  action  given  the  limit 
constraints 
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(36) 


These  limit  actions  are  represented  by  asymmetric  saturation  nonline¬ 
arities,  with  the  1/2  area  reference  case  being  a  positive  saturation 
limit  and  the  1  reference  condition  being  a  negative  saturation  limit 


The  resulting  equations  contain  the  values  of  Qoi/2'  Qoi«  etc. 
which  can  be  combined  to  provide  effective  values  of  K3,  a  representa¬ 
tion  of  the  effective  equilibrium  leakage  areas  ALoq/2  an<^  AL0^'  etc- 


These  equations  can  be  simplified  somewhat  by  elimination  of  the  non¬ 
linear  product  terms  pi  and  p2 ,  as  discussed  previously,  and  also 
simplified  by  means  of  a  control  action  that  requires 


-  Cj  (k-^z  +  k2u) 

(37) 


The  values  of  and  C2  are  determined  from  the  actual  fan  curves  in 
the  "fit"  prescribed  by  Eq.  (35)  in  such  a  way  that  the  effective 
control  action  is  then  made  symmetrical,  i.e.  a  symmetric  saturation 
limit  as  shown  in  Figure  2.  For  the  case  considered  previously,  in 
terms  of  the  fan  maps  given  in  Figure  4,  it  can  be  shown  that  the 
effective  representation  of  the  control  input  can  be  given  by 

=  -  A^  (k^z  +  k2w)  (38) 


with  the  limits 


-  1  < 


AA 


-  A 


<  1 


ref 


(39) 


This  particular  control  input  is  then  included  in  the  equations 
in  the  same  manner  as  the  case  for  variable  blade  angle  variation  for 
axial  fans,  with  the  term  Aa/aref  in  Eq .  (28)  replaced  by  AA/Aref  for 
the  present  case  of  variable  area  centrifugal  fans.  Calculated 
results  obtained  in  both  of  these  cases,  assuming  the  same  ship,  ALq 
value,  operating  conditions,  etc.  will  yield  the  same  numerical 
value  for  controlled  responses  when  using  the  same  k^  and  k2  values. 

RESULTS  OF  COMPUTATIONS 


Computations  were  carried  out  in  order  to  determine  the  resulting 
heave  accelerations  for  different  operating  conditions  (i.e.  speed- 
sea  state  combinations)  for  a  representative  2000  ton  SES  design.  The 
computations  covered  a  louver  system  with  the  total  deck  area  opening 
extending  up  to  the  equilibrium  leakage  area  value  Al0  (approximately 
87.5  ft?),  as  well  as  the  two  cases  of  fan  control  via  changing  blade 
angle  (axial  fans)  or  inlet  area  (centrifugal  fans) . 

The  results  for  the  louver  control  can  be  used  to  provide  values 
for  the  responses  and  associated  gains  for  the  case  of  variable  area 
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fans  operating  symmetrically  (about  the  3/4  open  condition,  as  express¬ 
ed  in  Eq.  (30)~(34)),  since  the  range  of  the  limit  given  by  Eq.  (34)  is 
il/2  Al0  -  The  signs  of  the  gains  have  to  be  reversed  for  the  fan  area 
control  case,  as  indicated  previously.  For  the  biased  fan  operating 
conditions,  i.e.  half  the  fans  at  1/2  open  and  the  other  at  full  open 
condition  for  the  reference  conditions,  the  results  are  obtained  sub¬ 
ject  to  the  limit  given  by  Eq.  (28)  and  (29),  which  is  the  same  as  the 
case  for  axial  fan  blade  angle  control  represented  by  Eq.  (28)  and  (24) 

The  computations  were  carried  out  for  an  operating  range  of  a 
representative  2000  ton  craft  in  head  seas,  at  a  speed  of  40  kt.  in 
Sea  State  4;  a  speed  of  50  kt.  in  Sea  States  4  and  5;  and  a  speed  of 
60  kt.  in  Sea  States  4  and  5.  The  results  obtained  are  given  in 
Tables  1-5  for  the  louver  control  case,  as  well  as  for  the  symmetric 
fan  area  control  case,  for  the  range  of  louver  area  opening  up  to 
±1/2  Al0  and  the  fan  area  change  AA/Aref  in  the  range  ±1/2.  The  values 
of  the  gains  kj.  and  k2  given  in  those  tables  apply  to  the  louver  con¬ 
trol,  where 

Al  =  j~  (k^ztkju!  =  52  (k^z+kjW)  (40) 


For  the  symmetric  fan  area  control  the  results  are  represented  by 

~  52  (k, z  +  k,u )  (41) 

°  Aref  1  1 

using  the  same  values  of  kj  and  k2  tabulated  in  Tables  1-5,  with 
essentially  the  same  expected  heave  acceleration  values.  Similar 
considerations  apply  to  the  case  of  the  "biased"  fan  operation,  with 
the  larger  limit  condition  on  fan  area  change,  as  well  as  the  axial 


fan  blade 

area 

control,  with 

the  results  in  those 

cases  given 

Tables  6- 

10. 

Table 

1.  Results  of  Control  Analysis, 

Louver  and 

x, , 

damping 

ratio 

h 

Symmetric 

u  =  40 

(uncontrolled 

a 

Fan  Area  Control 

kt.,  Sea  State  4 
runs  heave  accel,  = 

rms 

heave 

accel.,  g’s 

0 . 20g) 

k2 

0.5 

0 

13.604 

0.0878 

-  7  .796 

1 

17.261 

0.0533 

-10.328 

2 

20.934 

0.0394 

-12.208 

3 

23.225 

0.0317 

-13.773 

1.0 

0 

18.106 

0.0670 

-12.882 

1 

19.450 

0.0422 

-17.947 

2 

22.031 

0.0320 

-21.706 

3 

23.883 

0.0261 

-24.836 

1.5 

0 

21.062 

0.0545 

-17.968 

1 

21.532 

0.0351 

-25.565 

2 

23.310 

0.0270 

-31.204 

3 

24.748 

0.0223 

-35.899 
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Table  2.  Results  of  Control  Analysis,  Louver  and 
Symmetric  Fan  Area  Control 

u  *  50  kt.,  Sea  State  4 
(uncontrolled  rms  heave  accel.  =  0.25g) 


c , 

rms 

damping 

y 

heave 

lr 

ratio 

_i 

a 

accel.,  g's 

2 

0.5 

0 

14.673 

0.1084 

-  7.796 

1 

17.654 

0.0681 

-10.328 

2 

21.106 

0.0512 

-12.208 

3 

23.327 

0.0416 

-13.773 

1.0 

0 

19.115 

0.0814 

-12.882 

1 

20.126 

0.0527 

-17.947 

2 

22.416 

0.0405 

-21.706 

3 

24.132 

0.0334 

-24.836 

1.5 

0 

21.965 

0.0656 

-17.968 

1 

22.335 

0.0432 

-25.565 

2 

23.855 

0.0337 

-31.204 

3 

25.137 

0.0280 

-35.899 

Table 

3.  Results  of 

Control  Analysis, 

Louver  and 

Symmetric 

Fan  Area  Control 

u  =  50 

kt..  Sea  State  5 

(uncontrolled  rms  heave  accel.  = 

0 . 31g) 

e. 

rms 

damping 

heave 

y 

ratio 

_i 

a 

accel.,  q's 

*2 

0.5 

0 

29.825 

0.1725 

-  7.796 

1 

47.006 

0.1211 

-10.328 

2 

75.166 

0.1148 

-12.208 

3 

105.292 

0.1145 

-13.773 

1.0 

0 

45.635 

0.1560 

-12.882 

1 

59.244 

0.1133 

-17.947 

2 

85.744 

0.1068 

-21.706 

3 

114.975 

0.1061 

-24.836 

1.5 

0 

61.147 

0.1484 

-17.968 

1 

75.220 

0.1116 

-25.565 

2 

101.088 

0.1044 

-31.204 

3 

129.854 

0.1027 

-35.899 
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Table  4. 


Results  of  Control  Analysis,  Louver  and 
Symmetric  Fan  Area  Control 


u  =  60  kt . ,  Sea 
(uncontrolled  rms  heave 


State  4 
accel.  = 


0 . 27g) 


?,  rms 


damping 

V 

heave 

ratio 

_i 

a 

accel.,  g's 

2 

0.5 

0 

15.559 

0.1292 

-  7.796 

l 

18.033 

0.0834 

-10.328 

2 

21.277 

0.0636 

-12.208 

3 

23.426 

0.0521 

-13.773 

1.0 

0 

19.918 

0.0958 

-12.882 

1 

20.736 

0.0635 

-17.947 

2 

22.782 

0.0493 

-21.706 

3 

24.374 

0.0410 

-24.836 

1.5 

0 

22.663 

0.0767 

-17.968 

1 

23.021 

0.0513 

-25.565 

2 

24.290 

0.0404 

-31.204 

3 

25.430 

0.0339 

-35.899 

Table  5.  Results  of  Control  Analysis,  Louver  and 
Symmetric  Fan  Area  Control 

u  =  60  kt..  Sea  State  5 
(uncontrolled  rms  heave  accel.  =  0.38g) 


C, 

rms 

damping 

V 

heave 

k 

ratio 

2i 

a 

accel.,  g's 

k2 

0.5 

0 

32.568 

0.2118 

-  7.796 

i 

48.673 

0.1524 

-10.328 

2 

76.533 

0.1453 

-12.208 

3 

106.598 

0.1455 

-13.773 

1.0 

0 

49.887 

0.1924 

-12.882 

1 

63.329 

0.1431 

-17.947 

2 

89.499 

0.1352 

-21.706 

3 

118.576 

0.1344 

-24.836 

1.5 

0 

67.223 

0.1844 

-17.968 

1 

82.074 

0.1416 

-25.565 

2 

107.945 

0.1327 

-31.204 

3 

136.662 

0.1305 

-35.899 
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Table  6.  Results  of  Control  Analysis,  Biased  Fan 
Area  Control  {2  groups)  and  Axial  Fan 
Blade  Angle  Control 

u  =  40  kt..  Sea  State  4 
5,  rms 


damping 

ratio 

^1 

<7 

heave 

accel.,  g's 

k2 

0.5 

0 

13.598 

0.0878 

-  7.796 

1 

17.150 

0.0530 

-10.328 

2 

20.389 

0.0383 

-12.208 

3 

22.153 

0.0302 

-13.773 

1.0 

0 

17.958 

0.0664 

-12.882 

1 

19.144 

0.0415 

-17.947 

2 

21.263 

0.0308 

-21.706 

3 

22.617 

0.0247 

-24.836 

1.5 

0 

20.553 

0.0530 

-17.968 

1 

20.883 

0.0340 

-25.565 

2 

22.216 

0.0257 

-31.204 

3 

23.197 

0.0208 

-35.899 

Table  7.  Results  of  Control  Analysis,  Biased  Fan 
Area  Control  {2  groups)  and  Axial  Fan 
Blade  Angle  Control 

u  =  50  kt..  Sea  State  4 
C ,  rms 


damping 

ratio 

*1 

a 

heave 

accel.,  g's 

k2 

0.5 

0 

14.657 

0.1083 

-  7.796 

1 

17.518 

0.0675 

-10.328 

2 

20.529 

0 . 0497 

-12.208 

3 

22.226 

0.0396 

-13.773 

1.0 

0 

18.874 

0.0802 

-12.882 

1 

19.727 

0.0516 

-17.947 

2 

21.556 

0.0388 

-21.706 

3 

22.787 

0.0315 

-24.836 

1.5 

0 

21.270 

0.0633 

-17.968 

1 

21.506 

0.0415 

-25.565 

2 

22.599 

0.0318 

-31.204 

3 

23.446 

0.0261 

-35.899 
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c 

t 

damping 

ratio 

Table 

h. 

8.  Results  of  Control  Analysis, 
Area  Control  (2  groups)  and 
Blade  Angle  Control 

u  =  50  kt..  Sea  State  5 

rms 

heave 

a  accel . ,  g ' s 

.  Biased  Fan 
Axial  Fan 

k2 

in 

o 

0 

28.487 

0.1623 

-  7.796 

1 

37.435 

0.0943 

-10.328 

2 

46.276 

0.0689 

-12.208 

3 

51.991 

0.0553 

-13.773 

1.0 

0 

38.565 

0.1275 

-12.882 

1 

42.153 

0.0778 

-17.947 

2 

48.684 

0.0582 

-21.706 

3 

53.478 

0.0473 

-23.836 

1.5 

0 

45.531 

0.1056 

-17.968 

1 

46.939 

0.0666 

-25.565 

2 

51.663 

0.0507 

-31.204 

3 

55.541 

0.0417 

-35.899 

?  , 

damping 

ratio 

Table  9.  Results  of  Control  Analysis, 
Area  Control  (2  groups)  and 
Blade  Angle  Control 

u  =  60  kt . ,  Sea  State  4 

rms 

.  heave 

K1  a  accel . ,  q ' s 

,  Biased  Fa 
Axial  Fan 

k2 

0.5 

0 

15.527 

0.1289 

-  7.796 

1 

17.869 

0.0826 

-10.328 

2 

20.666 

0.0617 

-12.208 

3 

22.296 

0.0495 

-13.773 

1.0 

0 

19.577 

0.0940 

-12.882 

1 

20.238 

0.0618 

-17.947 

2 

21.830 

0.0471 

-21.706 

3 

22.950 

0.0385 

-24.836 

1.5 

0 

21.797 

0.0735 

-17.968 

1 

22.018 

0.0491 

-25.565 

2 

22.935 

0.0380 

-31.204 

3 

23.673 

0.0314 

-35.899 
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Table 

10.  Results  of  Control  Analysis, 

,  Biased  Fan 

Area 

Control  (2  groups)  and 

Axial  Fan 

Blade 

Angle  Control 

u 

=  60  kt..  Sea  State  5 

t . 

rms 

damping 

]f 

heave 

Jr 

ratio 

21 

0 

accel.,  q's 

*2 

0.  S 

0 

30.525 

0.1950 

-  7.796 

i 

38.136 

0.1165 

-10. 328 

2 

46.595 

0.0862 

-12.208 

3 

52.191 

0.0697 

-13.773 

1.0 

0 

40.606 

0. 1509 

-12.882 

1 

43.455 

0.0945 

-17.947 

2 

49.434 

0.0716 

-21.706 

3 

53.977 

0.0586 

-24.836 

1.5 

0 

47.477 

0.1241 

-17.968 

1 

48.589 

0.0799 

-25.565 

2 

52.789 

0.0616 

-31.204 

3 

56.362 

0.0510 

-35.899 

Comparing  the  results  obtained  by  these  different  methods,  it  can 
be  seen  that  the  greatest  heave  acceleration  reduction  occurs  when 
using  the  variable  area  fan  in  the  biased  manner,  i.e.  as  two  separate 
groups,  or  by  use  of  the  axial  fans  with  variable  blade  angle  control. 
This  is  due  to  the  larger  range  of  effective  "area"  change  possible  in 
that  case,  since  the  saturation  effect  at  half  that  range  (for  the 
louver  or  symmetric  variable  area  fan)  limits  the  benefits  of  the  other 
techniques.  Of  course  similar  benefits  could  be  obtained  by  using  a 
larger  range  of  louver  area,  with  an  attendant  larger  louver  area  bias 
level,  but  that  would  immediately  demand  greater  lift  system  power. 

The  best  control  system  must  really  be  evaluated  via  more  detailed 
computer  simulation,  where  power  requirements  are  also  determined,  but 
the  present  results  do  indicate  the  expected  region  of  benefit  so  that 
the  required  detailed  computer  simulation  studies  of  such  control 
effects  can  be  limited  to  the  more  promising  types  of  control  concept. 

Examination  of  the  results  given  in  the  tables  shows  that  in 
general  the  rms  heave  accelerations  are  reduced  further  as  the  gains 
k.  and  k2  are  increased.  However,  in  some  results  obtained  in  other 
studies  the  effect  of  the  saturation  is  exhibited,  so  that  an  optimum 
value  for  that  particular  operating  case  is  evident.  The  tables  also 
indicate  the  degree  of  difference  due  to  the  particular  value  of  the 
damping  ratio,  with  the  greater  acceleration  reduction  indicated  for 
the  largest  value  of  damping  ratio.  However  good  dynamic  response  to 
transient  disturbances  is  required  and  too  large  a  value  of  q  would 
give  a  sluggish  transient  response.  Since  the  linearized  damping  ratio 
would  not  be  fully  attained  throughout  the  motion  time  history  of  the 
controlled  craft,  the  most  suitable  value  of  damping  ratio  should  be 
C  =  1.0  for  satisfactory  transient  response. 

While  the  values  tabulated  in  the  tables  indicate  a  significantly 
large  reduction  of  the  rms  acceleration,  it  is  not  expected  that  such 
effects  would  be  realized  in  actual  practice  (or  even  with  the  large 
scale  digital  computer  simulation) .  This  is  due  to  the  presence  of 
other  forces  acting  on  the  craft  that  have  been  neglected,  but  more 
important  is  the  effect  of  leakage  arising  from  the  craft  motions. 

Such  leakage  is  expected  to  occur  primarily  due  to  pitch  motion  that 
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leads  to  leakage  gaps  at  the  bow  due  to  the  relative  motion  with 
respect  to  the  sea  surface  at  that  location.  The  present  analysis 
does  not  consider  these  important  effects,  but  it  is  expected  that  the 
aim  of  designing  for  adequate  transient  response,  and  also  with  an 
effective  increase  in  the  craft  natural  frequency,  would  be  beneficial 
in  reducing  the  responses  due  to  impulsive  disturbances  that  arise  due 
to  transient  leakage.  The  main  use  of  the  data  presented  in  the  tables 
is  therefore  to  establish  a  control  rule,  with  initial  gain  value 
estimates,  for  use  in  the  necessary  large  scale  digital  computer 
simulation  studies  for  arriving  at  an  optimum  RCS  for  SES  craft  in 
waves . 

FULL  SCALE  TESTCRAFT  RESULTS  WITH  RIDE  CONTROL  SYSTEM 

While  the  computational  results  given  in  the  tables  apply  to  a 
representative  2000  ton  SES  design,  such  craft  are  not  presently  built 
and  operating.  However  100  ton  testcraft  have  been  operating  for  some 
time  in  order  to  obtain  data  that  can  be  used  as  a  means  of  predicting 
the  expected  characteristics  of  full  scale  behavior,  as  well  as 
functioning  as  a  test-bed  for  evaluating  various  concepts  that  could 
be  applied  to  multi-thousand  ton  SES  craft  in  the  near  future. 

Some  illustrative  data  that  demonstrate  the  effectiveness  of  the 
RCS  system  using  the  louver-type  control  mechanism  is  shown  in  Figure 
5.  In  that  case  the  time  histories  show  a  specific  illustration  of  the 
general  reduction  in  heave  acceleration  and  cushion  pressure  when  the 
RCS  is  activated.  Further  detailed  information  that  illustrates  the 
benefit  obtained  from  the  use  of  an  RCS  system  is  given  by  Figures  6 
and  7,  where  power  spectra  of  heave  acceleration  and  cushion  pressure 
variations,  both  with  and  without  the  RCS  system  operating,  are  shown. 
These  illustrations  of  actual  full  scale  performance  of  100  ton  test¬ 
craft  with  RCS  are  indicative  of  the  nature  of  the  benefits  that  can 
be  obtained  with  such  a  system.  Continued  analysis  and  prediction 
using  the  initial  methods  described  in  this  paper,  together  with  more 
detailed  large  digital  computer  simulation  studies,  will  allow  develop¬ 
ment  of  an  optimal  RCS  for  use  in  the  presently  developing  multi-thou¬ 
sand  ton  SES  craft. 
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Figure  6.  Mission  111  Power  Spectrum  of  z  (cabin). 
Head  Sea  State  2,  RCS  On  and  Off. 
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ment  of  lightweight  variable  geometry  high 
performance  fans. 

HENRY  J.  DEMATTIA  is  presently  a  program 
manager  in  the  Electromagnetics  and  Acoustics 
Division  (SEA  034)  of  the  Naval  Sea  Systems 
Command.  He  received  a  B.S.  Degree  in 
Mathematics  from  West  Virginia  Wesleyan 
College  and  Graduate  Studies  at  the  George 
Washington  University.  Prior  to  current 
position  with  NAVSEA,  he  was  a  Ship  Design 
Project  Coordinator  (elec,  ronlcs)  for  the 
Naval  Ship  Engineering  Center  (NAVSEC),  and 
also  held  engineering  positions  with  RCA 
Service  Company,  Alexandria,  VA  and  Inter¬ 
national  Electronics  Engineering  Incorporated , 
Washington,  DC. 

JAMES  W.  DONNELLY  Senior  Project  Engineer  was 
born  in  Philadelphia,  PA  on  June  20,  1944  . 

He  received  B.S.  and  M.S.  degrees  in  electri¬ 
cal  engineering  from  Drexel  University  in 
1967  and  1969,  respectively.  From  1969  to 
present,  he  has  been  employed  in  the  Machinery 
Automation  Systems  Department  of  the  Naval 
Ship  Engineering  Center,  Philadelphia  Div. 

He  was  made  co-recipient  of  the  NAVSECPHILA- 
DIV  Technical  Achievement  Award  in  recognition 
of  the  technical  support  that  he  provided  in 
the  development  and  evaluation  of  the  control 
system  for  the  FFG-7  prototype  propulsion 
system. 
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A.M.  DORRIAN  is  a  Chartered  Engineer  and 
Member  of  the  Institution  of  Mechanical 
Engineers.  In  1968  he  obtained  an  H.N.C.  with 
distinction  in  Mechanical  Engineering  and  in 
1970  graduated  from  Strathclyde  University 
with  an  Honours  Degree  in  Mechanical  Engineer¬ 
ing.  He  completed  engineering  apprenticeship 
with  G&J  Weir  Ltd,  Joined  Y-ARD  Ltd  in  1970 
and,  following  work  on  a  variety  of  acoustic 
and  vibration  problems  joined  the  Controls 
Group  where  he  has  been  responsible  for  the 
development  of  the  machinery  controls  philo¬ 
sophy  for  a  variety  of  naval  and  special  pur¬ 
pose  ships.  He  is  currently  Head  of  the 
Surface  Ship  Control  Systems  SectioQ. 

KARE  HARALD  DRAGER  holds  the  M.Sc.  Electrical 
Engineering  (1966)  and  M.Sc.  Industrial 
Engineering  (1973).  Following  studies  of 
Ship  Automation  with  computers  reliability 
analysis  of  bridge  instrumentation  and  systems, 
he  is  currently  Project  Manager  of  a  team 
studying  "Cause  Relationship  of  Collisions  and 
Groundings"  for  Det  norske  Veritas. 


ARTHUR  E.  DRESSER  obtained  the  BSME  in  1965  and 
is  a  Registered  Professional  Engineer  in  the 
State  of  Maine.  He  is  presently  a  Senior 
Mechanical  Engineer  with  Bath  Iron  Works  with 
responsibilities  including  Shipboard  Automation 
Project  Engineering. 


ALBERT  DUBERLEY  served  an  apprenticeship  in 
Marine  Engineering.  He  is  a  chartered  engineer 
and  a  member  of  the  Institution  of  Mechanical 
Engineers.  Joined  Admiralty  Engineering  Lab 
(now  National  Gas  Turbine  Establishment)  in 
195^  where  he  is  now  employed  as  a  Principal 
Scientific  Officer  and  Head  of  the  Machinery 
Control  Section. 


A 1 1-7 


HARUZO  EDA  is  a  Senior  Research  Engineer  at 
Davidson  Laboratory  and  a  Research  Associate 
Professor  for  the  Department  of  Ocean  Engineer¬ 
ing,  Stevens  Institute  of  Technology.  Employed 
at  Davidson  Laboratory  since  1961,  Dr.  Eda  is 
actively  engaged  in  research  on  stability  and 
control  of  ships  and  underwater  bodies.  Dr. 

Eda  serves  as  a  Member  of  the  Panel  H-10  (Ship 
Controllability)  of  the  Society  of  Naval  Archi¬ 
tects  and  Marine  Engineers;  a  Member  of  the 
Maneuvering  Committee  of  International  Towing 
Tank  Conference  and  of  the  American  Towing 
Tank  Conference,  and  also  as  a  Member  of  SNAME 
and  the  Society  of  Sigma  Xi .  Dr.  Eda  holds 
Bachelor  and  Doctoral  Degrees  from  Osaka 
University  in  Japan,  and  a  Master's  Degree 
from  Stevens  Institute  of  Technology. 

EIVIND  KINDEN  ENGEBRESTSEN,  a  1956  graduate 
from  the  University  of  London,  has  been 
associated  since  1965  with  the  control  and 
instrumentation  section  of  the  Ship  Research 
Institute  of  Norway.  Leader  of  the  section 
since  1973.  Mr.  Engebrestsen  is  in  charge  of 
design  specification  and  commissioning  of 
ship  automation  systems. 


R.J.  EVES  graduated  in  Electronic  and  Electri¬ 
cal  Engineering  at  University  of  Birmingham, 
U.K.,  1965.  Has  been  employed  by  Lucas  Aero¬ 
space  since  1969  and  has  been  engaged  in 
systems  design  and  analysis  of  hydromechanical 
and  electronic  control  systems  for  Gas  Turbine 
Engines  for  Aircraft,  Industrial  and  Marine 
applications.  Present  position  -  Performance 
Project  Engineer. 


JAMES  A.  FEIN  graduated  in  1969  from  the 
University  of  Michigan  with  a  B.S.E.  in 
Engineering  Mechanics.  He  has  been  employed 
at  the  Naval  Ship  Research  and  Development 
Center  from  1969  to  the  present.  At  this 
time  he  is  a  naval  architect  in  the  Ship 
Performance  Department  concerned  with  the 
dynamics  of  high  performance  vehicles.  He 
received  an  M.S.  degree  in  Mechanical  Engineer¬ 
ing  in  1972  from  the  University  of  Maryland 
where  he  is  currently  pursuing  a  PhD  in  the 
same  field.  He  is  a  member  of  American 
Institute  of  Aeronautics  and  Astronautics  and 
the  Society  of  Naval  Architects  and  Marine 
Engineers . 
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RAY  FERANCHAK,  after  graduation  from  Youngstown 
University,  started  work  at  Westlnghouse , 
developing  systems  and  controls  for  d-c  drives 
for  towboats,  dredges,  and  Icebreakers.  Since 
then,  he  has  been  involved  with  electrical 
drives  for  wind  tunnels,  missile  pumps,  air¬ 
craft  generators,  and  transit  vehicles.  He 
has  done  graduate  work  in  industrial  control 
systems,  servo-mechanisms,  and  feedback  control 
systems.  He  is  a  member  of  IEEE  and  has 
published  several  papers  on  aircraft  generator 
test  drives. 

ALAN  S.  FIELDS  obtained  B.E.E.  from  the  City 
College  of  New  York  and  a  M.E.E.  from  Drexel 
Institute  of  Technology.  He  has  been  employed 
at  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  since  196*4,  involved  in  the 
design  of  steering,  tracking  and  station¬ 
keeping  control  systems  for  surface  ships. 


CDR  EUGENE  E.  FITZPATRICK,  USN  is  a  1959 
graduate  of  the  U.S.  Naval  Academy  and  holds 
the  MSME  and  degree  of  Naval  Engineer  from 
MIT.  CDR  Fitzpatrick  served  in  the  Naval 
Reactors  program  from  1966  to  1973,  served 
four  years  in  ship  overhaul  management  in 
the  Naval  Sea  Systems  Command  and  is  presently 
Project  Coordinator,  Machinery  Division,  Naval 
Ship  Engineering  Center. 


FLOYD  FOLTZ  has  been  with  Hamilton  Standard 
for  the  past  16  years.  He  is  presently  a 
Project  Manag  -r  responsible  for  the  FT-9 
Engine  Condition  Monitor  and  for  Hamilton 
Test  Systems  TRENDS  systems.  His  major 
programs  have  been  data  acquisition  systems 
and  control  systems,  which  provided  a  back¬ 
ground  of  sensor  applications,  data  acquisi¬ 
tion  techniques,  computer  Interfaces,  and 
field  evaluation.  Mr.  Foltz  was  with  the 
U.S.  Air  Force  for  seven  years. 
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C.  FRENCH  served  an  electrical  apprenticeship 
at  a  Royal  Dockyard  meanwhile  obtaining  a  B.Sc. 
in  Electrical  Engineering.  Following  a  period 
of  Industrial  and  University  research  was 
awarded  a  Ph.D.  in  Engineering.  Undertook 
design  work  on  large  turbo-alternators  with  a 
power  station  plant  manufacturer,  then  re¬ 
entered  the  academic  world  on  the  staff  of  a 
Polytechnic  followed  by  a  University.  Rejoined 
the  MOD  in  the  Ship  Department  six  years  ago, 
spending  the  last  two  as  Head  of  the  Systems 
Development  Group  in  the  Machinery  Controls 
Section. 

IRA  FRIEDMAN  received  his  B.S.  degree  in 
Physics  at  Brooklyn  College.  Over  the  past  25 
years  he  has  held  positions  of  responsibility 
in  navigation,  guidance,  and  control  of  air 
and  surface  platforms,  and  in  Combat  System 
Integration  for  the  3K  SES  (Surface  Effect 
Ship).  He  was  a  member  of  the  NSIA  Committee 
Investigating  new  ASW  Weapon  Concepts  for  high 
speed  ships.  He  Joined  Ingalls  Shipbuilding 
Division,  Litton  Systems,  Inc.  in  1977  where 
he  has  been  involved  in  Propulsion  Control, 

Ship  Data  Multiplexing,  and  Advanced  Combatant 
studies.  He  holds  a  U.S.  Patent  and  is  a 
member  of  the  IEEE. 

HANS  FRIVOLD  graduated  from  the  University  „f 
Oslo  in  1973  with  a  degree  in  Applied 
Mathematics  Hydrodynamics.  From  1973  to  1977 
Mr.  Frivold's  work  was  in  propeller  hydrodyna¬ 
mics,  the  pressure  field  due  to  working  pro¬ 
peller.  This  past  year  his  work  has  been  with 
mathematical  models  for  ship  maneuvering 
simulation. 


DANIEL  H.  GEORGE  received  the  B.S.E.E.  degree 
from  Louisiana  State  University  in  1968.  He 
has  worked  for  the  Tennessee  Valley  Authority 
at  the  Browns  Ferry  Nuclear  Power  Plant  in 
Athens,  Alabama.  Since  1970  he  has  been 
employed  by  the  Panama  Canal  Company  Engineer¬ 
ing  Division  as  a  design  engineer,  and  is 
presently  a  general  engineer  with  wide  respon¬ 
sibility  for  major  studies  and  projects.  Mr. 
George  is  a  member  of  IEEE,  Tau  Beta  Pi,  and 
a  Registered  Professional  Engineer  in  the 
Canal  Zone. 
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F.  WAYNE  GIBSON  attended  the  Royal  Military 
College  of  Canada  in  Kingston,  Ontario,  where 
he  obtained  a  B.Eng.  in  Engineering  Physics 
In  1971.  He  has  been  a  Naval  Officer  in  the 
Canadian  Armed  Forces  since  being  commissioned 
in  1971  and  has  served  as  a  Combat  Systems 
Engineering  Officer  in  HMCS  KOOTENAY  and  a 
Project  Systems  Engineer  at  National  Defence 
Headquarters,  Ottawa,  Ontario. 


CECIL  N.  GOFF  holds  a  BME  from  the  University 
of  Florida  and  the  MSE  from  Pennsylvania  State 
University.  Since  1968  Mr.  Goff  has  been 
employed  at  the  Naval  Training  Equipment  Center 
where  he  has  worked  with  a  Periscope  View 
Simulator,  a  DD963  Class  Propulsion  Trainer 
and  a  1200  PSI  Propulsion  Plant  Trainer. 


FABIO  R.  GOLDSCHMIED  holds  the  degrees  of  BS 
and  MS  in  Mechanical  Engineering  from  Columbia 
University.  For  the  past  ten  years  Mr. 
Goldschmled  has  been  an  Advisory  Scientist  in 
the  Research  and  Development  Center,  Energy 
Systems  Division,  Westinghouse  Electric  Corp. 
He  is  the  author  of  numerous  technical  papers 
for  the  AIAA  Journals  of  Aircraft  and  Hydro- 
nautics,an  Associate  Fellow  of  AIAA  and  a 
member  of  Technical  Committee  on  Marine 
Systems  and  Technologies. 


SIGURD  GORANSSON  graduated  in  Electrical 
Engineering  in  1959  from  Chalmers  University 
of  Technology.  In  196A  he  joined  the  staff 
of  the  Swedish  State  Shipbuilding  Experimental 
Tank,  where  he  is  in  charge  of  the  Hybrid 
Computer  Facilities. 
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LINDSAY  GORRELL  obtained  the  B.A.Sc.  degree 
In  Mechanical  Engineering  from  the  University 
of  Waterloo  in  19 69  and  the  M.Sc.  degree  in 
Ergonomics  from  the  University  of  London  In 
1970.  Mr.  Gorrell  is  presently  employed  as  a 
Defence  Scientist  by  the  Defence  and  Civil 
Institute  of  Environmental  Medicine  in 
Toronto.  Much  of  his  work  involves  human 
engineering  evaluation  and  design  of  electro¬ 
nic  display  systems.  He  is  a  member  of  the 
Human  Factors  Society  and  the  Human  Factors 
Association  of  Canada. 

LYMAN  GRISWOLD  graduated  from  Duke  University 
(BSEE  55)  and  received  a  Master  of  Science  in 
Engineering  from  George  Washington  University 
in  1962.  He  recently  retired  from  Federal 
Service  as  Head  of  the  Sensors  Branch  at  the 
David  W.  Taylor  Naval  Ship  Research  and 
Development  Center.  He  is  a  Captain  in  the 
United  States  Navy  Reserve  and  a  contributor 
to  Bowdltcb  (The  American  Practical  Navigator) 
on  the  subject  of  speed  logs. 


EIGIL  HAALAND  graduated  in  1967  from  the 
Norwegian  Institute  of  Technology.  Following 
six  years  of  research  work  at  the  Institute 
Mr.  Haaland  joined  Norcontrol  where  he  worked 
on  marine  control  systems.  Since  1978,  he 
has  been  an  R&D  section  leader,  instrumenta¬ 
tion  department. 


DONALD  D.  HALL  is  presently  assigned  to  the 
Tactical  Computer  System  Architecture  Division 
at  the  Naval  Ocean  Systems  Center  since  1975, 
where  he  has  served  as  a  microcomputer  net¬ 
work  system  design  engineer.  He  received  his 
BSEE  degree  from  the  University  of  the  Pacific 
and  MSEE  from  the  University  of  California. 
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MAKOTO  HARA  was  born  in  Yamaguchi,  Japan,  in 
1951-  He  graduated  at  Dept,  of  Electronics, 
Hosei  Junior  College  in  1971.  From  1971  to 
1972,  he  was  engaged  in  Kitashiba  Electric  Co. 
Ltd.  He  has  been  with  Dept,  of  Control 
Engineering,  Tokyo  University  of  Mercantile 
Marine  since  1972,  where  he  is  presently  a 
research  assistant.  His  main  interest  is  the 
application  of  microprocessors  to  practical 
marine  control  systems. 


THURMAN  R.  HARPER  is  a  Senior  Simulation 
Analyst  with  Propulsion  Dynamics,  Inc., 
Annapolis,  Maryland.  He  has  participated  in 
and  directed  propulsion  dynamics  and  control 
simulation  studies  of  gas  turbine,  steam, 
diesel  and  combined  ship  propulsion  plants. 

He  has  also  been  employed  at  the  David  W. 
Taylor  Naval  Ship  Research  and  Development 
Center,  Annapolis  engaged  in  ship  automation, 
control  and  mathematical  model  development  of 
gas  turbines  and  ship  systems.  Dr.  Harper 
holds  the  B.S.E.E.  degree  from  Drexel  Univer¬ 
sity,  Philadelphia  and  the  Ph.D.  degree  in 
Electrical  Engineering  from  the  Johns  Hopkins 
University,  Baltimore. 

GRAHAM  JOHN  HINDMARSH  B.Sc.,  M.B.C.S.  obtained 
a  first  class  honours  degree  in  computer 
science  from  the  Council  for  National  Academic 
Awards  while  studying  at  Wolverhampton  Poly¬ 
technic  U.K.  For  the  last  3  years  Mr. 
Hindmarsh  has  been  Manager  of  the  Computing 
and  Simulation  Department  of  the  Products 
Division  of  Vosper  Thornycroft  (UK)  Limited. 

In  this  capacity  he  is  responsible  for  the  use 
of  computers  in  engineering  applications  with 
particular  emphasis  on  marine  machinery 
control,  and  for  the  computer  simulation  of 
ship  hull  and  machinery  performance  leading  to 
the  definition  of  an  optimum  control  specifi¬ 
cation. 

DR.  DAN  HOFFMAN  received  his  undergraduate 
degree  in  Naval  Architecture  from  the  Univer¬ 
sity  of  Strathclyde  in  Glasgow  in  1959.  He 
then  proceeded  with  post  graduate  work  in  the 
Technion,  Israel  Institute  of  Technology  and 
the  Technical  University  of  Norway.  He  has 
been  a  Research  Professor  at  Webb  Institute 
of  Naval  Architecture  for  over  ten  years  and 
President  and  Technical  Director  of  Hoffman 
Maritime  Consultants  Inc.  (HMC)  over  the 
past  five  years.  Dr  Hoffman  has  authored 
numerous  papers  in  the  area  of  vessel  response 
to  the  environment  and  system  approach  to 
operations  in  rough  weather. 
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GEORGE  E.  HOLLAND  received  the  Bachelor  of 
Science  degree  from  the  U.S.  Merchant  Marine 
Academy  and  Bachelor  of  Mechanical  Engineering 
degree  from  the  University  of  Virginia.  Fol¬ 
lowing  active  duty  in  the  Navy  as  Main  Propul¬ 
sion  Assistant  on  USS  FREMONT  (APA-lUl)  he  has 
been  employed  as  a  mechanical  engineer  by  Esso 
Research  &  Engineering  Co.,  Chief  Engineer  by 
Renick  &  Mahoney,  Inc.,  and  Resident  Engineer 
by  Parker  Aircraft  Co.  before  joining  Machinery 
Systems  Division  of  the  Naval  Ship  Engineering 
Center  in  1966.  Since  July  1975  he  has  been 
Head,  Project  Coordination  Office. 

MICHIO  HORIGOME  was  born  in  Miyagi,  Japan  in 
1940.  He  received  the  B.S.  degree  in  Marine 
Engineering  from  Tokyo  University  of  Mercantile 
Marine  in  1963.  Since  1965  he  has  been  with 
the  Dept,  of  Control  Engineering,  Tokyo  Uni¬ 
versity  of  Mercantile  Marine,  where  he  is 
presently  an  Associate  Professor.  He  received 
the  Annual  Award  in  1977  from  the  Marine 
Engineering  Society  in  Japan.  His  research 
interests  are  in  time  series  analysis,  sto¬ 
chastic  process  and  the  application  of  these 
methods  to  Marine  Engine  Systems. 

DR.  THOMAS  M.  HOULIHAN  is  a  native  of  Buffalo, 
NY.  He  received  his  BSME  from  Manhattan 
College  and  his  Ph.D.  from  Syracuse  University. 
Since  1969  he  has  been  teaching  fluid  mechanics 
and  instrumentation  in  the  Mechanical  Engineer¬ 
ing  Department  of  the  Naval  Postgraduate  School 
He  has  received  the  Teetor  Award  of  the  SAE  and 
the  Dow  Award  of  the  ASEE  for  his  teaching 
endeavors . 


M.J.  JOBY  graduated  in  Mechanical  Engineering 
at  University  College,  London  1962.  Has 
since  specialized  in  Control  Systems  ranging 
from  powered  artificial  limbs  to  supersonic 
engines.  Has  worked  for  Lucas  Aerospace  since 
1966.  Present  position  -  Performance  Manager 
(Electronics ) , 
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Naval  Architect  at  Delft  University. 
Following  roil  it  ary  service  lie  has  been  since 
1970,  research  engineer  at  the  Institute  for 
Perception  TNI'  an i  ir.  charge  of  the  ship- 
manoeuvrinr  simulator. 


RAYMOND  C.  SIMAN'WIT':  v ■  ■  ■  ■  ’  v >■  i  fine  decree  of 
3.9  Mathematics,  I  re/..  '.  ’  a  1  v-rsi t.v  and  Mb’ 
Numerical  Science,  .'d. u.'  "  r  eins  University, 
ile  was  employed  at  *.!,•  ■  '  i  W.  Taylor  Naval 

Ship  Research  and  s  it  Center  (DTKSRDC) 

from  19(1-19 , 8  performing  as  a  Computer 
Systems  Analyst.  His  areas  of  expertise 
include  system  design,  software  design  and 
data  base  design. 


JONATHAN  SMALLEY  graduated  RSc  from  Durham 
University  and  MSc  from  Birmingham  University. 
He  joined  EASAKS  Ltd.  as  a  Senior  Systems 
Analyst  in  1976  with  main  interests  in  the 
design  of  Man/Kachine  Interfaces.  He  is 
currently  working  on  the  Man/Machine  Inter¬ 
face  design  study,  which  is  the  subject  of 
this  paper,  with  responsibility  for  human 
factors  and  related  design  aspects. 


BRIAN  SPENCER  graduated  in  1950  with  an 
honours  degree  in  mathematics.  He  became  a 
member  of  the  Scientific  Advisers  group  in  the 
Ship  Department  at  Bath  in  1965  with  partic¬ 
ular  responsibility  for  automatic  control 
systems  and  hydrodynamics.  Has  been  very  much 
involved  with  this  series  of  Ship  Control 
Systems  Symposia,  presented  a  paper  at  the 
First,  organized  the  Third,  a  committee  member 
for  the  Fourth,  co-author  of  a  paper  to  be 
presented  at  the  Fifth,  and  U.K.  coordinator. 
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C.  JOSEPH  B'JBi is  iiv.;;  i«-i'  of  ;  ••  ;  i  :  .•  :  ■. 
Dynamics,  t no .  ,  Annapo'.  is ,  Mary an  i 
Is  engaged  in  simulation  a:id  ana  :  y  t  3  ■. 

of  ship  propulsion  and  v<‘  gat!-  ri  .  .  '  •  ;  .  . 
Previous  to  his  present  j  .:iti  he  un¬ 
scientific  Consultant  and  Hoad  of  the 
Systems  Branch,  David  ft.  Tayl  or  Nava’.  ' ' 
Research  and  Development  center,  Ar;nap"lis 
Laboratoi'y ,  for  over  10  years.  Mr.  Rubis 
received  his  MS  decree  fror.  the  University 
Illinois.  He  was  tlie  technical  co-.-bairr.an 
of  the  Second  Shin  Control  .".'/stems  Symposium 
in  1969- 

LCDR  JAKES  K.  Rl'LANP,  graduated  from 
Tona  College  and  received  his  "aster's  Deere--' 
in  Business  Administration  from  Face  "r.Iver- 
sltv.  He  lias  served  on  the  rss  LACONIC 
(ACC-17),  l.'SS  PERRY  ( DD  ft-,),  ;rAuy  (CO-D 
and  USS  NASHVILLE  ( I 1  (a— 13).  He  also  serve) 
as  officer  in  charge  of  FCr-65  in  Vietmar. 
in  1968.  iie  is  currently  aaslgnei  to  the- 
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Electrical  Engineering  from  the  Illinois 
University  and  the  MS  in  Engineering  from 
U.C.L.A.  Employed  by  the  U.S.  Navy  since 
I960,  Mr.  Dejka  is  presently  a  Senior 
Scientist,  Naval  Ocean  Systems  Center, 
where  he  has  been  developing  long-range 
plans  for  the  Navy  in  the  areas  of  Dis¬ 
tributed  Systems  and  Computer  Science. 

Mr.  Dejka  is  a  senior  member  of  IEEE  and 
Sigma  Xi.  He  is  the  author  of  several 
professional  papers  and  reports. 
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WILLIAM.  M.  ELLSWORTH  is  a  graduate  of 
the  State  University  of  Iowa  where  he  ob¬ 
tained  B.S.  and  M.S.  degrees  in  Engineer¬ 
ing,  majoring  in  fluid  mechanics.  Since 
October  1969  he  has  been  Associate  Tech¬ 
nical  Director  for  Systems  Development, 
and  Head,  Systems  Development  Department, 
Naval  Ship  Research  and  Development  Center. 
He  is  a  licensed  Professional  Engineer  in 
the  State  of  Maryland;  a  member  of  ASNE, 

S NAME,  and  ASME;  and  author  of  a  number 
of  papers  and  reports  in  the  field  of 
naval  engineering.  He  served  as  a  member 
of  the  ASNE  Council  from  1972  to  1974  and 
was  nominated  for  Vice  President  in  1975. 

In  1974  he  received  the  ASNE  gold  medal. 
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CAPTAIN  (N)  E.  J.  (ED)  HEALEY  is  a  gradu¬ 
ate  of  the  Canadian  Military  College 
Royal  Roads  and  the  Royal  Naval  Engineer¬ 
ing  College  (RNEC)  Manado  at  Plymouth, 
England.  He  served  at  sea,  being  engi¬ 
neering  officer  of  three  East  Coast  war¬ 
ships,  interspersed  with  shore  postings 
at  the  Fleet  School  ani  in  the  Dockyard. 
Following  a  year  at  Staff  College,  Captain 
Healey  was  posted  to  NAVSEC  Philadelphia 
were  he  carried  out  the  shore  testing  of 
the  DDH  280  main  machinery.  Moving  to 
Ottawa  he  became  Design  Section  Head  for 
Propulsion  Machinery  and  Project  Officer 
for  the  DDH  280  Propulsion  Machinery. 
Following  successful  sea  trials  of  the 
DDH  280  Class  he  was  posted  to  the  Esqui- 
malt  Naval  Base  as  the  Command  Technical 
Officer.  This  was  followed  by  a  posting 
to  the  Canadian  Forces  Training  Command 
and  a  year  at  the  National  Defence  College. 
He  is  presently  Director  of  Marine  and 
Electrical  Engineering  at  National  Defence 
Headquarters . 


DONALD  H.  KERN  was  educated  at  the  Massa¬ 
chusetts  Institute  of  Technology  where  he 
received  his  Bachelor  of  Science  Degree 
in  Naval  Architecture  and  Marine  Engi¬ 
neering  and  his  Professional  Engineers 
Degree  in  Naval  Engineering.  During  his 
career  in  the  U.S.  Navy,  he  commanded  the 
Portsmouth  Naval  Shipyard  and  served  as 
Project  and  Design  Manager  for  various 
submarine  related  programs.  Currently 
Captain  Kern  is  Vice  President  of  Special¬ 
ized  Systems,  Inc.,  of  Mystic,  CT. 
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REAR  ADMIRAL  JAMES  W.  LISANBY  graduated 
from  the  U.S.  Naval  Academy  and  holds  an 
advanced  degree  in  Naval  Engineering 
(Architecture)  from  the  Massachusetts 
Institute  of  Technology.  Rear  Admiral 
Lisanby  is  an  Engineering  Duty  Officer 
with  wide  and  varied  experience,  both  at 
sea  and  in  shore  assignments.  From  1970 
to  1973,  he  served  as  Supervisor  of  Ship¬ 
building  at  Pascagoula,  Mississippi,  with 
contract  administration  responsibilities 
for  both  the  DD  963  and  the  LHA  1  ship 
acquisitions.  Following  a  brief  tour  as 
Assistant  for  Ship  Design  in  the  Office 
of  the  Chief  of  Naval  Operations,  he 
served  as  Project  Manager  for  the  LHA 
Class  of  Amphibious  Assault  Ships  with 
Headquarters  in  Washington,  DC,  from  197^ 
until  he  assumed  his  present  position  as 
Commander,  Naval  Ship  Engineering  Center, 
in  1978.  Rear  Admiral  Lisanby  holds  the 
World  War  II  Victory  Medal,  the  Korean 
Service  Medal  with  two  battle  stars,  the 
United  States  Medal,  and  the  National 
Defense  Service  Medal.  In  addition  to 
these  campaign  ribbons,  he  has  been  awarded 
the  Meritorious  Service  Medal  for  out¬ 
standing  service  with  the  Commander,  Naval 
Ship  Systems  Command. 

M.  DALE  MARTIN  is  a  BSME  graduate  from 
West  Virginia  University.  He  began  his 
government  career  as  a  project  engineer  in 
the  Bureau  of  Ships,  and  through  a  series 
of  progressively  more  responsible  assign¬ 
ments,  became  head  engineer  for  cargo 
systems  with  the  Fast  Deployment  Logistic 
Ship  Project,  followed  by  head  engineer 
for  systems  performance  for  the  same 
project.  Later  he  became  Director  of  the 
NAVMAT  Maintenance  Technology  Office,  fol¬ 
lowed  by  Director  of  the  Automatic  Test 
Equipment  Management  and  Technology  Office. 
This  office  was  recently  reorganized  into 
the  present  Test  and  Monitoring  Systems 
Program  Office,  where  Mr.  Martin  now  serves 
as  the  Technical  Director. 
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JOHN  TODD  Me LANE  obtained  a  B.S.  degree  in 
Physics  at  Lafayette  College  in  1951.  He 
joined  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center  in  1963  and  has 
studied  human  factors  engineering  at  the 
Catholic  University  of  America.  While  at 
the  Center  he  has  conducted  psychophysical 
experiments  in  night  vision,  displays/con¬ 
trols,  and  color  recognition.  He  has  also 
studied  the  effects  of  ship  motion  on 
sailors  and  the  human  engineering  design 
of  Navy  ship  bridges.  He  is  currently  in¬ 
volved  in  the  man/machine  problems  of 
shipboard  machinery  control  monitoring 
systems.  He  is  the  Center's  coordinator 
for  human  factors  technology. 


NILS  H.  M0RRB1N  was  born  in  Nykoping,  Sweden 
in  1926.  He  holds  the  Civilingenjor  (Naval 
Architecture),  Tekn.  License  and  Tekri.  Dr. 
From  1950  to  1951*  he  was  with  the  Submarine 
Design  and  Preliminary  Ship  Design  Department 
of  the  Royal  Swedish  Naval  Administration. 
Since  1955  he  has  been  with  the  Research 
Department  of  the  Swedish  State  Shipbuilding 
Experimental  Tank  as  Head  of  the  Ship  Dynamic 
Section . 


MICHAEL  0.  PARSONS,  Associated  Professor 
of  Naval  Architecture  and  Marine  Engineer¬ 
ing  at  the  University  of  Michigan  received 
a  Ph.D.  in  Applied  Mechanics  specializing 
in  optimal  control  from  Stanford  University 
in  1972.  He  returned  to  the  University  of 
Michigan  In  1972  as  Assistant  Professor 
with  teaching  assignments  in  marine  engi¬ 
neering,  automatic  control  and  computer- 
aided  ship  design.  His  current  research 
Interests  include  optimal  control,  opti¬ 
mization,  and  computer  applications  in 
naval  architecture  and  marine  engineering. 
Dr.  Parsons  is  a  member  of  Phi  Eta  Sigma, 
Tau  Beta  Pi,  Phi  Kappa  Phi,  Quarterdeck 
Society,  and  Scabbard  and  Blade  Society. 
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ANTON  CHARLES  PIJCKE  entered  the  Royal 
Netherlands  Naval  College  at  Den  Helder 
(branch:  Marine  Engineering)  in  1949  and 
received  his  commission  as  an  officer  in 
1952.  During  his  sea  duty  he  served 
mostly  onboard  destroyers  and  frigates. 

He  has  been  senior  lecturer  in  Marine 
Engineering  at  the  Royal  Netherlands  Naval 
College  during  several  years  and  was  also 
Head  of  the  Naval  Engineering  Department. 

He  left  the  Royal  Netherlands  Navy  as  a 
commander  and  is  now  a  staff  member  of  the 
Netherlands  Maritime  Institute  at  Rotter¬ 
dam.  Mr.  Pijcke  obtained  his  M.Sc.  degree 
at  London  University,  is  a  Fellow  of  the 
Institute  of  Marine  Engineers,  and  is  a 
chartered  engineer.  He  is  author  of  ap¬ 
proximately  twenty  published  papers.  He 
was  General  Chairman  of  the  4th  Ship 
Control  Systems  Symposium  held  at  The 
Hague,  The  Netherlands  in  1975. 

CAPTAIN  FREDERICK  P.  SCHUBERT,  Deputy  Chief, 
Office  of  Marine  Environment  and  Systems,  U.S. 
Coast  Guard  Headquarters,  graduated  from  U.S. 
Coast  Guard  Academy  in  1951  and  received  Mas¬ 
ters  degrees  from  Southern  California  and  Auburn 
Universities.  He  has  served  27  years  in  var¬ 
ious  operational  assignments  including  duties 
in  aviation  and  afloat  as  well  as  several  staff 
assignments  involving  program  management 
respons ib 11 it ies  .  Captain  Schubert  assumed 
fils  present  assignment  in  February  of  1977. 

His  duties  include  broad  program  management 
responsibilities  in  Port  Safety  and  Security, 
Marine  Environmental  Protection,  Aids  to  Navi¬ 
gation  and  Bridge  Administration.  Captain 
Schubert  served  as  a  member  of  the  1977 
Department  of  Transportation  Task  Force  on 
Marine  Oil  Transportation  and  Oil  Pollution  and 
was  Chairman  of  the  Pollution  Response  Working 
Group  on  the  Interagency  Task  Force  on  Tanker 
Safety  which  developed  the  background  for 
President  Carter's  Marine  Oil  Pollution 
Initiatives  of  17  March  1977. 


.CIAM  SPENCER  graduated  in  1950  with  an 
honours  degree  in  mathematics.  He  became  a 
member  of  the  Scientific  Advisers  group  in  the 
Ship  Department  at  Bath  in  1965  with  partic¬ 
ular  responsibility  for  automatic  control 
systems  and  hydrodynamics.  Has  been  very  much 
involved  with  this  series  of  Ship  Control 
Systems  Symposia,  presented  a  paper  at  the 
First,  organized  the  Third,  a  committee  member 
for  the  Fourth,  co-author  of  a  paper  to  be 
presented  at  the  Fifth,  and  U.K.  coordinator. 
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RICHARD  A.  STANKEY  received  the  B.S.  de¬ 
gree  from  the  Illinois  Institute  of  Tech¬ 
nology,  Chicago,  and  a  M.S.  degree  from 
Catholic  University,  Washington,  DC,  in 
Engineering  Sciences  and  Mechanical 
Engineering,  respectively.  He  came  to 
the  Naval  Ship  Engineering  Center  (NAVSEC) 
in  Washington,  DC  in  1966  as  a  hydraulic 
system  engineer.  In  1970  he  transferred 
to  the  Automatic  Controls  Section  of  the 
Ship  Control  Systems  and  Equipment  Branch 
(Code  6165)  of  NAVSEC.  He  has  done  system 
analysis  and  design  evaluation  on  hydrofoil 
flight  control  systems,  submarine  depth 
control  and  hovering  systems,  and  surface 
ship  position  keeping  and  track  keeping 
systems.  He  is  currently  working  on  the 
TRIDENT  and  Mine  Countermeasure  (MCM)  Ship 
Control  Systems. 
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REAR  ADMIRAL  J.  G.  C.  VAN  DE  LINDE  began 
his  higher  education  at  the  University  of 
Delft  in  chemistry  but  was  transported  to 
Berlin  to  a  forced  labor  camp  in  19^3-  He 
escaped  in  19M  and  spent  the  rest  of  the 
war  in  the  Dutch  underground.  In  19^6  he 
went  to  the  Royal  Netherlands  Naval  Col¬ 
lege  and  was  commissioned  as  an  officer  in 
19^8.  He  sailed  in  various  submarines  and 
surface  ships  and  saw  action  in  the  Korean 
war.  He  specialized  in  torpedoes,  weapon 
technology  and  nuclear  reactor  physics. 
After  a  four  year  command  of  the  Royal 
Netherlands  Naval  College  in  Den  Helder 
he  became  Chief  of  Naval  Personnel  in  1976. 


H.  R.  VAN  NAUTA  LEMKE  graduated  in  Electrical 
Engineering  at  the  Delft  University  of  Tech¬ 
nology  in  1950.  From  1950  -  1959  he  worked  in 
industry  (Van  der  Heem,  Philips)  in  the  field 
of  servo  systems,  instrumentation  and  sonar 
equipment.  From  1959  he  is  a  professor  in 
control  engineering  at  the  Delft  University. 
His  interests  are  in  the  area  of  adaptive  and 
optimal  systems  and  the  applications  of  fuzzy 
sets . 
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LCDR  IR.  WILLEM  VERHAGE  was  born  In  Rotterdam 
in  19^0.  After  completing  high  school  he 
enrolled  in  the  school  for  merchant  navy  in 
Amsterdam.  He  joined  the  Royal  Dutch  Navy  in 
1961  and  served  as  an  instructor  at  the  Naval 
Communications  School  in  Amsterdam.  In  19 69 
he  was  granted  his  engineering  degree  in  naval 
architecture  from  Delft  University  of  Tech¬ 
nology.  Since  1970  Ir.  Verhage  serves  at  the 
RNNC  in  Den  Helder  where  he  holds  the  rank  of 
LCDR.  In  addition  to  his  teaching  duties  in 
the  department  of  naval  architecture  he  is  head 
of  a  research  team  which  developed  a  nocturnal 
simulator . 
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